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1993 Lewis Workshop 


I on 

Computational Turbulence Modeling 

1. Objectives 

The purpose of this meeting was to discuss the current status and future de- 
velopment of turbulence modeling in computational fluid dynamics for aerospace 
propulsion systems. Various turbulence models have been developed and applied 
to different turbulent flows over the past several decades and it is becoming more 
and more urgent to assess their performance in various complex situations. In order 
to help users in selecting and implementing appropriate models in their engineering 
calculations, it is important to identify the capabilities as well as the deficiencies of 
these models. This also benefits turbulence modelers by permitting them to further 
improve upon the existing models. 

This workshop was designed for exchanging ideas and enhancing collaboration 
between different groups in the Lewis community who are using turbulence models 
in propulsion related CFD. In this respect this workshop will help the Lewis goal 
of excelling in propulsion related research. 

This meeting had seven sessions for presentations and one panel discussion over a 
period of 2 days. Each presentation session was assigned to one or two branches (or 
groups) to present their turbulence related research work. Each group was asked 
to address at least the following points: current status of turbulence model appli- 
cations and developments in the research; progress and existing problems; requests 
about turbulence modeling. The panel discussion session was designed for organiz- 
ing committee members to answer management and technical questions from the 
audience and to make concluding remarks. 
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Main Issues Raised in the Workshop 

(summarized by workshop committee) 

A total of thirty five papers were presented at the workshop covering a wide 
spectrum of work related to propulsion systems as is evident from the contents of this 
proceedings. Several problems and issues related to the turbulence and transition 
modeling were raised in these presentations as well as in the panel discussion. These 
issues are briefly summarized below. 

1. Transition 

Flow over a turbine blade very often undergoes a transition from laminar to turbu- 
lent and accurate modeling of this phenomenon is very important. The transition in 
propulsion systems (named as bypass transition) has its own features and is caused 
by elevated free stream turbulence, surface roughness, etc. Current work at Lewis 
shows that turbulence models are not being modified for their use in transitional 
flows. Current work at Lewis using turbulence models to capture the transition 
have only limited success. Bypass transition models developed at CMOTT, which 
take into account the effect of intermittency, account for the effect of free-stream 
turbulence. The development of transition models which account for both the free 
stream turbulence and surface roughness effects is a pacing item for transitional 
flow calculations for propulsion systems. Much work still needs to be done in this 
area. 

2. Heat transfer 

Currently, most turbulent heat transfer calculations at Lewis use Reynolds anal- 
ogy (i.e. constant turbulent Prandtl number) coupled with the zero-equation or 
two-equation models for the velocity field. The turbulent Prandtl number, in gen- 
eral, is not constant for complex flows. Therefore, improvement in the heat transfer 
calculations requires a more accurate approach for calculating the turbulent thermal 
diffusivity. One such approach is to use the transport equations for the temperature 
fluctuations and its dissipation rate. Some work in this direction is in progress at 
CMOTT. 

Since, as noted above, transition occurs on turbine blades quite often, part of the 
reason for the poor predictions of heat transfer is the inability of current models 
to capture this transition. There is a need to improve the modeling of this type of 
transitional flow. 

3. Turbulent chemical reaction interactions 

In this workshop, two approaches were presented for the calculation of turbulent 
flows with finite rate chemical reactions. One is the pdf method for species in which 
the chemical production rate of species can be exactly treated without modeling. 
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The other is the conventional moment closure scheme in which the averaged chemical 
production rate of species is estimated by the mean temperature of the flow field. 
Further studies are needed for clarifying and understanding the effect of chemical 
reaction on turbulence and vice versa. DNS of simple turbulent reacting flows could 
be very helpful to sort out this issue. 

Studies on realistic chemical reaction models are also needed for simulating the 
interactions between turbulence and chemical reaction. 

4. Wall bounded flows 

Some presentations at this workshop showed that for high turbulent Reynolds 
number flows, the use of wall function as a boundary condition often produces rea- 
sonable results for even complex flows in-spite of its theoretical in-validity. In addi- 
tion, this approach can significantly reduce the number of numerical grid points 
which, sometimes, is critical for complex engineering calculations. It was also 
pointed out that sometimes a key ingredient of the flow, such as a separation bub- 
ble, is under y + = 40. Such flows, therefore, cannot be well simulated with the 
wall function approach and in such cases turbulence model equations must be inte- 
grated to the wall. ( This type of turbulence modeling is usually called low Reynolds 
number turbulence modeling.) 

Most low Reynolds number turbulence models presented in this workshop are 
not tensor ially invariant, i.e. the model equations contain the detailed coordinate 
parameters which require ad hoc treatments in applications to complex flow ge- 
ometries. In addition, existing fc-e models do not respond to the effect of pressure 
gradient very well. To address these problems, two efforts have been made at 
CMOTT. For the wall function approach, a modified new wall function is proposed 
to better account for the effect of pressure gradient. For low Reynolds number fc-e 
models, a Galilean and tensorially invariant realizable k-e model is developed which 
can also respond to the pressure gradient very well. Further testing of this model 
in complex flows is needed to find its potential as well as its deficiencies. 

5. Inflow condition 

Some calculations showed that results were very much dependent on the inflow 
conditions chosen for turbulence quantities (say, k and e). To obtain a good com- 
parison with experiments, the turbulence inflow conditions had to be adjusted. This 
raises the uncertainty about correctness of calculations as well as the performance 
of turbulence models. To resolve this problem, the determination of correct in- 
flow conditions by either experiments or other methods is important. In situations 
lacking experimental inflow turbulence quantities, methods to establish reasonable 
inflow conditions are critical and must be further studied. Turbulence model test- 
ing should not be based on the flows with unknown inflow conditions except for the 
situations in which the flows are self preserving. 

6. Numerical method 

It was noticed that different turbulence models were built in different codes. It 
is difficult to compare the performance of different models based on the results 
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produced by different codes. It is agreed that the validation and comparison of dif- 
ferent turbulence models must be based on the grid independent converged solutions 
produced by one code with the same/similar initial and boundary conditions. 

Turbulence equations are much more stiff and complicated than the Navier-Stokes 
equations. The study of numerical methods to effectively obtain the accurate solu- 
tion of turbulence equations is also very critical for turbulent CFD calculations. The 
collaboration between CFD and turbulence model developers will help to resolve 
this problem. 

7. Hierarchy of turbulence models and standards 

Several different types of turbulence models have been developed over the past 
several decades. Each of them was proposed for certain type of flows. In CFD calcu- 
lations, all models are one-point closure schemes which include zero-equation models 
(or mixing length models), one- and two-equation models, algebraic Reynolds stress 
models and second order closure models. At each level of these turbulence models, 
there are several versions (as many as 10 or more) of same type of models. Each 
later version was developed to improve the performance over its parent-model for 
the flows of interest. 

Engineers were often, and are still, confused by so many levels of turbulence 
models with so many versions at each level. Do we need all levels of turbulence 
models? Is there a best model at each level? The possible answer to above questions 
is yes. 

From zero-equation to second order closure models, each level of turbulence mod- 
els contains a certain degree of turbulence physics which gives each level of models 
a different degree of potential as well as limitation. For example, the mixing length 
models have the potential for modeling of attached boundary layer flows. While the 
good prediction of pressure distribution is often obtained, the predictions of skin 
friction and heat transfer are often poor, especially for separated flows. This is not 
surprising because the mixing length models were developed for that kind of flows 
and contain little turbulence physics. In this workshop, many presentations showed 
that two-equation isotropic eddy viscosity models perform much better than mixing 
length models for variety of flows. This is also not surprise because there are two 
turbulent transport equations involved in this type of models. However, for (lows 
in which the anisotropy of turbulence plays main role, we must use anisotropic 
turbulence models. In addition, if the evolution of Reynolds stresses, streamline 
curvature, rotation, etc. are important, the second order closure schemes will be 
more appropriate. 

In principle, the higher level of models (if they were well developed), give better 
performance in a wider range of turbulent flows and, of course, at the expense of 
higher computing cost and numerical difficulty. Further studies on the potential as 
well as deficiencies of each level of models in turbulent flow calculations are needed. 
Identifying the best model version at each level seems an urgent task for engineering 
CFD applications. To achieve this goal, systematic model testing for a variety of 
bench mark flows as well as real flows in propulsion systems is needed. 
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Research Activities at the Center for 
Modeling of Turbulence and Transition 

Tsan-Hsing Shih 


1. Introduction 

2. General developments 

2.1 Turbulent constitutive relations 

2.2 Mechanical and scalar dissipation equations 

2.3 Eddy viscosity transport equation 

3. One-point closure schemes 

3.1 One-equation eddy viscosity transport equation model 

3.2 Galilean and tensorial invariant realizable k-e model 

3.3 Reynolds stress algebraic equation model 

3.4 Scalar flux algebraic equation model 

3.5 Reynolds stress transport equation model 

3.6 Non-equilibrium multiple-scale model 

3.7 Bypass transition model 

3.8 Joint scalar PDF model 

4. RNG and DIA 

5. Numerical simulation 

Abstract 

The main research activities at the Center for Modeling of Turbulence and Tran- 
sition (CMOTT) are described. The research objective of CMOTT is to improve 
and/or develop turbulence and transition models for propulsion systems. The flows 
of interest in propulsion systems can be both compressible and incompressible, 
three dimensional, bounded by complex wall geometries, chemically reacting, and 
involve “bypass” transition. The most relevant turbulence and transition models for 
the above flows are one- and two-equation eddy viscosity models, Reynolds stress 
algebraic- and transport-equation models, pdf models, and multiple-scale models. 
All these models are classified as one-point closure schemes since only one-point (in 
time and space) turbulent correlations, such as second mo ment s (Reynolds stresses 
and turbulent heat fluxes) and third moments ( UiUjUk , urf 2 ), are involved. In 
computational fluid dynamics, all turbulent quantities are one-point correlations. 
Therefore, the study of one-point turbulent closure schemes is the focus of our 
turbulence research. However, other research, such as the renormalization group 
theory, the direct interaction approximation method and numerical simulations are 
also pursued to support the development of turbulence modeling. 
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1. Introduction 

The center for modeling of turbulence and transition was established as a special 
focus group within the Institute for Computational Mechanics in Propulsion at 
NASA Lewis Research Center in 1990. Its objective is to improve and/or develop 
turbulence and transition models for computational fluid dynamics (CFD) applied 
in propulsion systems. With the advance of computer technology and algorithms, 
accurate turbulence and transition modeling becomes the pacing item for improving 
flow calculations used in propulsion system design in all its key elements. The flows 
of interest in propulsion systems are, in general, very complex since there are wall- 
bounded three-dimensional complex geometries, chemical reactions, compressibility 
and transition, etc. In order to accurately predict these flows one must correctly 
model the turbulent stresses and scalar fluxes which are one-point (in time and 
space) turbulent correlations. For flows with finite rate chemical reactions, accurate 
modeling of the production rate of species is crucial for turbulent flow calculations. 
Based on the above considerations, turbulence modeling activities at CMOTT are 
focused on one-point closure schemes, that is, using the moment closure schemes for 
the turbulent velocity field and the joint scalar pdf method for the reacting scalar 
field. 

There are various moment closure schemes which have been developed for var- 
ious engineering applications. However, in practice, one often finds that the ex- 
isting models need to be improved and/or re-developed in order to reasonably 
simulate complex flow structures appearing in propulsion systems. For this pur- 
pose, CMOTT devotes itself to improving and/or re-developing these moment clo- 
sure schemes which include eddy viscosity (one- and two-equation) models, second 
moment algebraic- and transport-equation models, non-equilibrium multiple-scale 
models, and bypass transition models. In addition, other studies supporting the 
development of one-point closure schemes have been also carried out (for example, 
studies on renormalization group theory (RNG), direct interaction approximation 
(DIA), direct numerical simulation (DNS) and large eddy simulation (LES)). 

In this report, we first describe the general development of turbulent constitutive 
relations, turbulent mechanical and thermal dissipation and a new eddy viscosity 
equation. Second, we describe the detailed developments on each moment closure 
scheme and the pdf method. Then the RNG and DIA methods and finally, the 
numerical simulation of particular turbulence phenomena, such as rotation and 
bypass transition, etc., are considered. 

Each research subject is the joint project of several CMOTT researchers and 
visitors. In describing research activities, the names of involved researchers will be 
mentioned for reference. 
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2. General Developments 

2.1 Turbulent Constitutive Relations 


Reynolds stress 

Using the invariant theory in continuum mechanics and Generalized Cayley- 
Hamilton formulas for tensor products, a turbulent constitutive relation (or a gen- 
eral turbulence model) for any turbulent correlations can be obtained, in principle. 
Therefore, this theory provides an avenue to develop better turbulence models than 
those existing. For example, a commonly used constitutive relation for Reynolds 
stresses u^j (in terms of the mean deformation rate tensor Uj j and the turbulent, 
velocity and length scales characterized by the turbulent kinetic energy k and its 
dissipation rate e) is 


£ 


UiUj C ^ ^ ( Ui,j T ^ kS^j 


The effective eddy viscosity v? defined as 

l/rp = 


UiUj k" . • • 

— = Cp— for i ^ j 


Ui,j + Uj,i 


( 2 . 1 . 1 ) 


(2.1.2) 


is isotropic since ut is a scalar quantity. However, the invariant theory enables us 
to formulate the following general model (Shih and Lumley 1 , Johansson 2 ): 

2 iF 2 2 

UiUj — /c j T 2 q , 2 ^ ( Ui,j ~t- Ujj ^ Uj jSjj ) 

+ + (/?,- ^n,^) 

K z 1 

+ 2fte — (Ui'kUj'k — -n 2 <5 i3 ) 

£ O 

+ 2 a 7 ~(U k jU k ,j-ln 2 6 ij ) 


3 

K A 9 

+ 2a s -^(U i , k Ul k + Vl k Uj. k - -Ilstfij) 

+ 2 “io-j(u k ,,u kl + v Kj ul - |n 5 Sy) 

+ 2ai 2 ^(f7? 4 f7? l .-in 4 ^) 

J r2a n ‘^-(Ul,Ul J -\ IWy) 

K 5 9 

+ 2 a lt ~(U,, k U l , k Ul j + Uj.M.tUl - -IVy) 

K Q 9 

+ 2 * u - r (u iJ .v?. k u?j + - |n 6 «y) 

+ 2a w U {Ui kUl kU 2 mV 2 m + - ?n 7 «y) 


(2.1.3) 
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where 


III — Ui'kUk.i , II 2 = Ui, k U iJt , II 3 

n 5 = Uij.Uu.vl . n 6 = Ui, t ul k ul , 

From Eq.(2.1.3), the effective eddy viscosity 


Ui.kUl k , n 4 = ul k ulk , 
n, = v^u^vlJJlrr, 

(2,1.4) 


Mij 


UjU j 

Ui,j + Uj.i 


(2.1.5) 


is no longer a scalar and, hence, is an anisotropic eddy viscosity. It is noticed that 
the first two terms on the right hand side of Eq.(2.1.3) represent the standard k-e 
eddy viscosity model (2.1.1) and that the first five terms of Eq.(2.1.3) are of the same 
form as the models derived from both the two-scale DIA approach (Yoshizawa 3 ) and 
the RNG method (Rubinstein and Barton 4 ). 

Eq.(2.1.3) is a general model for uJu]. It contains 11 undetermined coefficients 
which are, in general, scalar functions of various invariants of the tensors in ques- 
tion, such as SijSjj (strain rate) and flijflij (rotation rate) which are (II 2 + IIi)/2 
and (n 2 — III ) /2 respectively. The detailed forms of these scalar functions must 
be determined by other model constraints, for example, realizability, and by exper- 
imental data. Eq.(2.1.3) contains 12 terms; however, its quadratic tensorial form 
may be sufficient for practical applications. We will see later in section 3.3 that the 
constitutive relation (2.1.3) has a significant impact on the development of Reynolds 
stress algebraic equation models. 


Turbulent scalar flux 9u{ 

We assume the following functional form: 


9ui 


Fi(Ui,j,T'i, fc,e, 9 2 ,eq) 


( 2 . 1 . 6 ) 


where 9 2 is the variance of a fluctuating scalar and e# is its dissipation rate. 
Eq.(2.1.6) indicates that the scalar flux depends on not only the mean scalar gradi- 
ent 7’,;, but also the mean velocity gradient Uij and the scales of both velocity and 
scalar fluctuations characterized by k, e, 9 ' 2 , s 0 . 

Applying the invariant theory, we may obtain the following general constitutive 
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relation for Ouf. 


0u i = a 1 k{--) 1/2 T, i + — (--) 1/2 (a 2 t/ i . j +a 3 C/ j<i )^j 

€ €o e £ £f) 

k 0^ 1 /9 

+ -2 ( ) 1 (a 4 U Uk U kJ + a 6 U jtk U k ,i + a & Ui, k U jM + a 7 U^U kJ )T^ 

£ £ £# 

u 4 u~a 2 , 

+ -j ( ) ' \asUi, k Ul k + a 9 Ul k Uj, k + a 10 )T.; 

£° £ £$ 

+ £ (* 

+ a,uUi' k Ui, k Ui j + ai* ) Uj, k Ui' k U’i i )Tj 
+ ^(- e ^) 1/2 (0i eU^Ufj + auUj^Ul) 


C O C Q 


(2.1.7) 

A; /fl 2 


The coefficients ai - aig are, in general, functions of the time scale ratio -/ 


e ' fo 


and the other invariants formed by the tensors in question, for example, T %k T %k , 
TjJJijTj, etc.. Again, Eq.(2.1.7) implies that the effective eddy diffusivity 


(7r)i = 


-0Ui 


is not isotropic. It is noticed that the first term on the right hand side of Eq.(2.1.7) 
is the standard eddy diffusion model, and the models derived from the two-scale 
DIA (Yoshizawa 5 ) and the RNG method (Rubinstein and Barton 6 ) are similar to 
the first two terms of Eq.(2.1.7). In practice, a form containing the first two terms 
on the right hand side of Eq.(2.1.7) may suffice. Further development of this model 
for turbulent heat transfer is described in Section 3.4. 

The Researchers involved with the subject in this section are T.-H. Shih, J. Zhu, 
A. Shabbir, J.L. Lumleyfand A. Johansson.-^ 

2.2 Mechanical and Scalar Dissipation Equation 
Mechanical dissipation £ 

In turbulence modeling, we often need turbulent characteristic velocity and length 
scales. While the turbulent kinetic energy k is used to characterize the velocity scale, 
the mechanical dissipation rate £ and the scalar dissipation rate eq are used to char- 
acterize the length scales for mechanical and scalar fields, respectively. Comparing 
with the turbulent kinetic energy equation, the exact dissipation rate equation is 

f Professor, Cornell University, Ithaca, NY 
| Professor, Royal Institute of Technology, Stockholm, Sweden 
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very complicated. In this equation, all the terms which represent important tur- 
bulence physics (for example, turbulent diffusion, generation and destruction) are 
unknown and are of complex forms that are all related to small scales of turbulence. 
Therefore, in the literature, the exact dissipation equation is not considered as a 
useful equation to work with. Instead, one creates a model equation by assuming an 
analogy to the turbulent kinetic energy equation, i.e., one assumes that the model 
dissipation rate equation also has generation and destruction terms which are as- 
sumed to be proportional respectively to the production and dissipation terms in 
the turbulent kinetic energy equation over the period of large eddy turn-over time 
characterized by k/e. The resulting model dissipation rate equation is written as 

e.t + UiS'i - (eui).i 

£ £ 2 

Cel ^UiUjUi'j Ce 2 

Recently, Lumley 7 proposed a dissipation rate equation based on the concept of 
spectral energy transfer caused by interactions between eddies of different sizes. 
This model equation mimics the physics of statistical energy transfer from large 
eddies to small eddies and is of a different form than equation (2.2.1). 

In this study, we explore another rational way to obtain the model dissipation 
rate equation which contains certain important physics and hope it will work better 
than the existing one. The idea is that first, there is a relationship between the 
dissipation rate e and the mean-square vorticity fluctuation uJiUl at high Reynolds 
numbers or in homogeneous turbulence: 


(2,2,1) 


£ = UU){UJi 


and second, all the terms appearing in the uJiJJJl equation have more clear physical 
meanings than that in the e equation so that the uJiUl equation is easier to model. 
Once the UJui equation is modeled, a model dissipation rate equation will be readily 
obtained. 

The exact equation for ZJiUTi is 


\ . . UJ jCJ: v . LO-iUJi , 1 . 

(— 2~ )•* + u ){~2~h = K-y-Xjj - 2 (<UjU}iU}i ^ j + 

U T OJiUJjUi'j T U>iU)jUi y j 


( 2 . 2 . 2 ) 


where u-, and U t are the fluctuating and mean velocities, and ui; and (1, are the 
fluctuating and mean vorticity which are defined by 


Ldi — £ijl c Uk,j — ^ijkUk.j (2.2.3) 

Tennekes and Lumley 8 clearly described the physical meaning of each term in equa- 
tion (2.2.2). Order of magnitude analysis shows that the first, third, fourth and fifth 
terms on the right hand side of Eq.(2.2.2) become small compared with all other 
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terms in the equation as the turbulent Reynolds number increases. The sixth and 
seventh terms are the production due to fluctuating vortex stretching and the dissi- 
pation due to the viscosity of the fluid. As the turbulent Reynolds number increases 
these last two terms become dominant and the balance between them determines the 
evolution of vorticity fluctuations. Neglecting terms UiUijQj, —ujuiflij, u fiUjUij 
and v(^r L ),jj, the evolution of UjUi at large Reynolds number will be described by 
the following equation, 


/ U{Ui \ . j-t / U%Ui v 1 , . /o o A \ 

( 2 )•* "f" ( 2 ' — — 2 + WiWjUij (2.2.4) 

To model UiUjUij — uuijUij, let us first estimate UiUjUij. we define an 
anisotropic tensor b^: 

/.!•/.» • 1 

(2.2.5) 


LOiiO-i 1 „ 

m = ^ U 


‘j 


ut 


*3 


then uiiUjUi'j can be written as 


UiWjUi.j = bfjUlUij 


( 2 . 2 . 6 ) 


We expect that the vortex stretching tends to align vortex lines with the strain rate 
so that the anisotropy 6^ would be proportional to the strain rate Sij, i. e., 


b'ij oc — , where s = (2 SijSij) 1 ^ 2 , = (uij + Uj^)/2 


(2.2.7) 


This leads to the following model: 


uJiUjUij a uil^SijSij) 1 / 2 a u\ ^2s lj s lJ 


(2.2.8) 


where we have assumed that u\ and (2 SfySy) 1 ^ 2 are well correlated. 

Using the relation, Ui = eijkUkj, it is not difficult to show that at large turbulent 
Reynolds number, 

UiUJi ~ 2sijSij (2.2.9) 

and Eq.(2.2.8) can be also written as 


UiUjUi,j (X u\ 



( 2 . 2 . 10 ) 


Equation (2.2.10) indicates that this term is of the order (u 3 fl 3 )R 3 / 2 as it should 
be. On the other hand, from eq.(2.2.4) the term UiUjUij - uu^Ju^~ must be of the 
order (u 3 // 3 )/?, which is the order of magnitude of all the other terms in Eq.(2.2.4), 
therefore the term -i/UijUij must cancel the term (2.2.10) or (2.2.8) such that the 
difference of these two terms is smaller than the term (2.2.10) or (2.2.8) by an order 
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1 /2 

of R, . This suggests that the combination can be modeled by 

the following two terms: 

(2,2.11) 

because the ratio of k/u to yjujj and the ratio of s to S are of order R 1 / 2 , where 

k « u 2 is the turbulent kinetic energy and S is the mean strain rate (2 SijSij) 1 / 2 . 
Equation (2.2.11) does give the right order of magnitude for — i/u~JuR~. 

Therefore, the dynamical equation for fluctuating vorticity (2.2.4) at large Reynolds 
number can be modeled as 


“I 


UJiUJ; v c OiUJi v 


1 k 




( 2 . 2 . 12 ) 
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Using £ = vu )ju)j, we readily obtain the following model dissipation rate equation, 

£ ^ 

£,t + Uj£'j = —(Uj€.)'j + C^lS £ — Cu 2 j y= (2.2.13) 

k + y l/£ 

where C u i and are the model coefficients which are expected to be constant at 
large Reynolds number. 

It should be noticed that Eq.(2.2.13) is different from the standard e equation 
(2.2.1) by both the generation and destruction terms. First, the Reynolds stresses 
do not appear in the generation term and the new form of the generation term is 
similar to that proposed by Lumley 7 which is based on the concept of spectral energy 
transfer. Second, the destruction term is well behaved so that equation (2.2.13) will 
not have a singularity anywhere in the flow field. We expect that equation (2.2.13) 
will be numerically much more robust than equation (2.2.1). 

Equation (2.2.13) can be applied to any level of turbulence modeling including 
second order closure models; however the turbulent transport term (eu7 ),; needs 
to be modeled differently at different levels of turbulence modeling. In an eddy 
viscosity model, the term (at7),i will be modeled as 


(ml* = -(—e,ih ( 2 . 2 . 14 ) 

The coefficients C u \, C w 2 , and the eddy viscosity v? must be calibrated using 

experimental data (Shih et al?) 


Scalar dissipation £q 

A similar analysis leads to the following model scalar dissipation rate equation: 
£e,t + Uj£oj = — (ujSg),] + Ce\S £0 + CqiPt §y/££e — C # 3 — - °:=— ( 2 . 2 . 15 ) 

K + y/V£ 

where $ = sJT^Tj and T is the mean scalar quantity, such as, the mean tempera- 
ture. Further development of heat transfer model is described in Section 3.4. 

The Researchers involved with the subject in this section are T.-H. Shih, W. Liou, 
A. Shabbir and Z. Yang. 
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2.3 Eddy Viscosity Transport Equation 

In eddy viscosity models, one accepts the following simple constitutive relation 


and assumes that the eddy viscosity is characterized by some kind of velocity and 
length scales u' and l: 


In two-equation fc-e eddy viscosity models, for example, one specifies that 


The eddy viscosity assumption (2.3.4) is commonly adopted in two-equation models. 
Eqs.(2.3.1) and (2.3.4) together with appropriate k and e equations have been widely 
used in engineering calculations. However, for cases where the mean flow changes 
quickly or has a strong mean stream-line curvature or rotation, etc., this kind of 
model does not work very well, since the assumption (2.3.4) is too simple to account 
for the effect of the above mean flow structure on eddy viscosity. 

The main purpose of this study is to drop the assumption (2.3.4) and to derive 
an exact equation for ut based on Eq.(2.3.1) and other exact turbulence equations 
(i.e. first principles). In this way, we hope that some important turbulent physics 
can be brought into the eddy viscosity and that a physically sound turbulence eddy- 
viscosity can be calculated. 

Using Eq.(2.3.1), we may write for incompressible flows 



(2.3.1) 


z/y oc u' t 


(2.3.2) 



(2.3.3) 


and, hence, the eddy viscosity is assumed as 



£ 


(2.3.4) 



(2.3.5) 


Differentiating both sides, we obtain 


Dt S 2 Dt 


D Sij D 

— ut = - -prf- — 



(2.3.6) 


The equation for muj can be written as 


Dt U i U j D ij U Pi.j E 1 1 ij Sij -}- Coij 


(2.3.7) 
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where 


D{j [uUiUj UjlljUk^'k 

Pij 'U'jUkUjfc UjUf; U Ik 

Hij 'P'l'lLj ”|“ PjUi 

Eij — ZvUi'kUj^ 

CO{j — 2£j m fcfi m Ufc'Uj 2£jmk^mUkUi 

Inserting Eq.(2.3.7) into Eq.(2.3.6), we obtain an exact transport equation for eddy 
viscosity 


= -§J(Dy + Pij + n« - e (j + c 0ij ) -~ 2 §- ( S 2 (2.3.8) 

In this equation, all the important turbulence physics in the Reynolds stress equa- 
tion, such as Reynolds stress diffusion term Djj, production term P^, pressure- 
velocity gradient correlation term IL^ and dissipation tensor Eij, are involved. 
Comparing with the standard eddy viscosity assumption (2.3.4), this exact eddy 
viscosity equation (2.3.8) contains very rich turbulence physics. This equation also 
implies that a second order closure model will naturally lead to a corresponding 
eddy viscosity model. 

Now, as an example, we use Launder Reece and Rodi’s 10 model and a gradient 
transport model for the triple velocity correlation (— mujUk — k ) to derive 

a model equation for vt- The resulting equation is 


D r , , vt \ T . / . 3 vt s u T.kS 2 k vt .2uxSijSij kk 

Ut a 2 a a (2.3.9) 

, 4 , ^ £ , ruri O) SikSkjSji V T D 2 

J r Tb k ~ C ^T + 2{C,-2 ) V T — 


Note that the Coriolis terms do not explicitly appear in this equation; however 
the rotation effect on vt could be carried over through the mean flow field. In 
addition, we also note that there are no extra model coefficients introduced in 
Eq.(2.3.9). All model coefficients (<r, C\ and C 2 ) are brought in from the second 
order closure model. The values of these model coefficients may need adjustment in 
model applications. Note that Eq.(2.3.9) is not a self-consistent equation since the 
turbulent kinetic energy k and its dissipation rate e are also involved. Eq.(2.3.9) 
together with k-e transport equations will provide a new three-equation model which 
may better represent the effect of mean flow structure as well as mean flow history 
on the eddy viscosity. 

The Researchers involved with the subject in this section are T.-H. Shih, Z Yang, 
and W. Liou. 


3. One-Point Closure Schemes 

In this section, we describe the developments on each of the moment closure 
scheme and the pdf method which are of concern at CMOTT. The first two sections 
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3.1 and 3.2 describe the one- and two-equation isotropic eddy viscosity models. Sec- 
tions 3.3 and 3.4 describe the new developments on Reynolds stress and scalar flux 
algebraic equation models. Section 3.5 assesses Reynolds stress transport equation 
models. Section 3.6 describes a multiple-scale model for non-equilibrium turbulence. 
Section 3.7 is about transition models. Finally, in Section 3.8 the pdf method for 
turbulent chemical reaction is described. 

3.1 One-equation eddy viscosity model 

Recently developed one-equation eddy viscosity models are either based on the 
assumption (Baldwin and Barth 11 ): 


V T = C M — (3.1.1) 

£ 

or created according to computational experience (Spalart and Allmaras 12 ). Both of 
them are successful in some flow calculations. This scheme is quite attractive in CFD 
because one only needs to solve one scalar v? equation without bothering about 
other turbulence quantities. However, comparing with k-e two equation models, the 
above mentioned one-equation vj models do not contain any more turbulent physics. 
In fact, Baldwin and Barth’s model is, basically, a change of dependent variable 
based on Eq.(3.1.1) plus some extra approximations. Therefore, in principle, we 
should not expect any superior performance over two-equation models. However, if 
we do not use the assumption (3.1.1), there is the possibility to improve and extend 
the capability of one-equation eddy viscosity models. 

The objective of this study at CMOTT is to derive a physically sound eddy 
viscosity equation which contains rich turbulent physics and accounts for various 
effects from mean flow structures. 

Note that in Section 2.3 we have already derived an exact equation for the eddy 
viscosity (2.3.8) and also a model equation (2.3.9) which is based on the Reynolds 
stress transport equation model of Launder, Reece and Rodi (LRR). All turbulent 
physics contained in the Reynolds stress equation can be brought into the eddy 
viscosity equation. Therefore, in principle, the transport equation (2.3.9) should be 
better than existing one-equation models based on Eq.(3.1.1). However, Eq.(2.3.9) 
is not self-consistent because k and e are also involved. To make Eq.(2.3.9) self- 
consistent, we must model k and k/e in terms of i/y and S. In most shear flows, 
the energy-containing eddy turn-over time k/e is of the same order as the mean 
flow time scale S' -1 , so that e/k oc S is a reasonable model. In addition, a crude 
dimensional analysis gives k oc vtS and this is, of course, reasonable only for shear 
flows. After the above considerations, the resulting self-consistent one-equation 
model is: 


— v T = [{y + — ){y T ),fc],fc — (y T ) M u t) , k + C u iS ut 

ut a a 


+ 2 (C„ 2 - 2 )rr 


Ft D n2 
2S 2 Dt 


(3.1.2) 
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where the diffusion terms from the Reynolds stress equation (2.3.7) have been ma- 
nipulated and approximated. Eq.(3.1.2) clearly exhibits the various effects of the 
mean flow on the eddy viscosity. 

The model coefficients C u i, C u 2 and a can be determined by using the experi- 
mental data of homogeneous shear flows, free shear flows and boundary layer flows 
as well as the relations in the inertial sublayer. Extensive tests of this model in 
various flows are carrying out at the CMOTT. 

The Researchers involved with the subject in this section are T.-H. Shih, W. Liou, 
Z. Yang and J. Zhu. 

3.2 Galilean and tensorial invariant realizable k-e model 

The two-equation k-e eddy viscosity model is one of the most widely used tur- 
bulence models in engineering calculations. The k-e model has versions for high 
Reynolds numbers and for low Reynolds numbers. For wall bounded turbulent 
flows, the high Reynolds number k-e model (for example, Launder and Spalding 13 ) 
must be applied together with a wall function as its boundary condition, while the 
low Reynolds number k-e model (for example, Jones and Launder 14 ) can be inte- 
grated to the wall. The high Reynolds number k-e model of Launder and Spalding 
is considered as a standard k-e model. We notice that even though the model dissi- 
pation rate equation is created by assuming an analogy with the turbulent kinetic 
energy, there was not much modification until Lumley 7 and Shih et al.° For near 
wall turbulence, in addition to Jones and Launder’s model, there are many other 
versions of low Reynolds number k-e models (such as Chien 15 , Shih and Lumley 16 , 
Yang and Shih 17 ) which have made better performance over Jones and Launder’s 
model. 

There are, probably, four or five issues worth mentioning about existing low 
Reynolds number k-e models: the model constants are not consistent with those 
in the high Reynolds number k-e model; the wall correction terms and damping 
functions are related to the wall distance so that models are not tensorial invariant; 
a nonrealistic dissipation rate near the wall is introduced; they are not always 
realizable since normal stress could become negative; and finally, they do not work 
very well for boundary layer flows with various pressure gradients. 

The objective of this study at CMOTT is to overcome the above mentioned 
problems. First, we propose a vorticity dynamics based dissipation rate equation 
as a part of high Reynolds number k-e base model. 9 Second, based on the invariant 
theory, inhomogeneous terms for the dissipation rate equation are proposed which 
enable the model to better respond to the change of pressure gradients (Yang and 
Shih 18 ). Third, the wall distance parameter is removed from the damping function 
so that the model is tensorially invariant (Yang and Shih 19 ). The model constants 
are consistent with those in the high Reynolds number k-e model. Finally, the 
non-negativity of normal Reynolds stresses, the realizability condition, is imposed. 

The Researchers involved with the subject in this section are Z. Yang, T.-H. Shih 
and C. Steffen. 
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3.3 Reynolds stress algebraic equation model 

All eddy viscosity models including one- and two-equation models are isotropic. 
For the flows where anisotropy is important, for example, the secondary flows driven 
by turbulent normal stresses in a square duct or curved duct, eddy viscosity mod- 
els do not produce correct flow structures. To overcome this intrinsic deficiency of 
isotropic eddy viscosity models, one proposes a Reynolds stress algebraic equation 
model which will provide an effective anisotropic eddy viscosity. The first such a 
model was proposed by Rodi 20 and it achieved some success in the prediction of 
anisotropic related flow structure. However, Rodi’s formulation is a set of algebraic 
non-linear system equations for Reynolds stresses and it often creates numerical dif- 
ficulty in obtaining a converged solution. Recently, Taulbee 21 obtained an explicit 
algebraic expression for the Reynolds stress using Pope’s 22 tensor expansion formu- 
lation and solved this numerical difficulty. However, in general, Rodi’s formulation 
assumes that the ratio ujuj/k is constant and, of course, this is not really true for 
most turbulent flows of interest. Therefore, sometimes, this Reynolds stress alge- 
braic equation model produces even worse results than the isotropic eddy viscosity 
models for cases where eddy viscosity models are appropriate. 

Alternative ways for obtaining effective anisotropic eddy viscosity models have 
been tried by a few researchers, for example, the DIA method by Yoshizawa 3 , 
the RNG method by Rubinstein and Barton 4 and invariant theory by Shill and 
Lumley. 1 . It is interesting to point out that the RNG and DIA methods result in 
the same formulation and that this formulation is the first five terms of a general 
constitutive relation Eq.(2.1.3) except that the model coefficients are different. 

One of our goals at CMOTT is to search for an effective anisotropic eddy vis- 
cosity model for complex turbulent flows where the nonequilibrium of turbulence 
is not very severe so that the constitutive relation (2.1.3) is more or less valid. We 
have explored the potential capability of Eq.(2.1.3) and found that a truncation of 
Eq.(2.1.3) up to the quadratic terms of the mean velocity gradients is sufficient for 
various flows of interest. The model coefficients are determined such that realizabil- 
ity for the normal stresses is ensured. The detailed analysis is described by Shih et 
alP 

The quadratic version of Eq.(2.1.3) together with the standard k-s transport 
equations, successfully predicts many complex flows as well as simple flows which 
include backward-facing step flows; confined coflowing jets; confined swirling coaxial 
jets; flows in 180° curved duct; flows in a diffuser and a nozzle; boundary layer 
flows with pressure gradient and turbulent free shear flows. See references 23-25 for 
detailed results. 

The Researchers involved with the subject in this section are J. Zhu and T.-H. 
Shih. 

3.4 Scalar flux algebraic equation model 

In parallel with Reynolds stress algebraic equation model, we have also tried to 
develop an effective anisotropic scalar eddy diffusivity model for scalar (heat) fluxes 
based on the new constitutive relation (2.1.7) and the new thermal dissipation rate 
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equation (2.2.16). We have determined that it seems sufficient to truncate Eq.(2.1.7) 
up to linear terms of the mean velocity gradient, i.e., 

Wi=a l k(--) 1 / 2 T, i + -(--) 1/2 (a 2 C/ i , i +fl 3 t/ i . i )T i (3.4.1) 

£ €0 £ £ £e 

This equation indicates that the heat flux and the mean temperature gradient are 
not necessarily in alignment due to the distortion of the flow field. This means that 
the effective scalar eddy diffusivity is anisotropic. 

Eq.(3.4.1) together with the 6 2 and £q equations will be a closed set of model 
equations for turbulent heat fluxes. The model coefficients are calibrated from 
homogeneous flows. Detailed analysis and a few model tests are described in this 
research briefs by A. Shabbir. 

The Researchers involved with the subject in this section are A. Shabbir and 
T.-H. Shih. 

3.5 Reynolds stress transport equation model 

The Reynolds stress transport equation model is considered as a next generation 
of advanced turbulence modeling for engineering applications. In principle, the 
second moment equations describe various effects of the mean flow and external 
agencies on the evolution of turbulence, hence, are the most attractive way (also 
the simplest correct way) to study turbulent flows. 

Various closure models for second moment equations have been developed. The 
success of these closures are marginal and vary with each flow. To identify the 
sources of their deficiencies, one often uses simple flows where the specific model 
term in the second moment equations can be isolated, hence, the corresponding 
model can be checked against experimental data or direct numerical simulation 
(DNS). For example, using pre-distorted anisotropic homogeneous relaxation flows, 
we may check the return-to- isotropy models with experimental data or DNS. How- 
ever, for other flows, several model terms, such as, triple velocity correlations, rapid 
and slow pressure-strain correlations, etc., simultaneously exist and can not be iso- 
lated in the experiments. In these cases (for example, in a homogeneous shear flow 
or a channel flow) only DNS can provide all the information for simultaneously 
checking various models. 

We have examined various existing closure models using experimental data as 
well as DNS data (Shih et al. 20 and Shih and Lumley 27 ). Conclusive statements are 
difficult to draw at this time. However, the following remarks can be made about 
various closures for the second moment equations, i.e., the triple velocity correlation 
Tjjk, the rapid and slow pressure related correlations n*J, n*j, and the dissipation 
rate tensor 

a) Tijk. All the existing models, such as Daly and Harlow 28 , LRR 1(I , Lumley 29 , 
etc., are not very satisfactory for highly inhomogeneous flows, such as flow near the 
wall. However, for flows where the inhomogeneity is not very high, the above closure 
models become close to each other and also closer to the DNS data. In addition, 
the triple velocity correlations in these situations are usually small comparing with 
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other terms in the equation, so that modeling of this term is not as critical as other 
terms for the results of turbulent flow calculations, except for the flow near the wall, 
b) n'J. It is very clear from all the available DNS data that nonlinear models, such 
as, Shih and Lumley 30 are much better than linear models, such as SSG 31 . It seems 
also that the following constitutive relation 

n[J = F(ujUj,Ui,j) 


is quite appropriate, i.e., its dependence on turbulent Reynolds number and other 
parameters is quite weak and can be neglected. However, one deficiency of this form 
observed by Reynolds 32 is that it can not take the rotation effect into account, 
c) II?j. This term is usually modeled together with the dissipation tensor e t j and 
the combination of the two is called the return-to-isotropy term. All existing mod- 
els are unsatisfactory at the present time. They are far from “universal”, i.e., their 
performance varies from flow to flow. It is noticed that some strange behavior of 
return-to-isotropy (for example, for some pre-distorted flow relaxation, turbulence 
evolves toward anisotropy before it returns to isotropy) occurs and cannot be pos- 
sibly modeled with the following constitutive relation: 

n -j = F(ufu], k, s ) 

In addition, the behavior of return-to-isotropy was found to depend not only on the 
Reynolds stresses at the present time but also on their history according to DNS 
data (Lee 33 ). It may be also necessary to include triple velocity correlations into 
the above constitutive relations from the definition of The term seems 
highly dependent on the turbulent Reynolds number and slowly approaches to its 
asymptote as Reynolds number goes to infinity, so that, in general, one should not 
exclude its dependence on turbulent Reynolds number even for moderate Reynolds 
numbers. In addition, n?j is also noticeably affected by the mean strain rate ac- 
cording to the DNS data 34 , so that, in general, the mean strain rate should be also 
considered in the constitutive relation. In short, much more research is needed for 
developing a better model of n?j. 

The Researchers involved with the subject in this section are T.-H. Shih and A. 
Shabbir. 

3.6 Non-equilibrium multiple-scale model 

To consider the effect of the nonequilibrium of energy spectrum on turbulent 
quantities, such as k, e and uiu], etc., Hanjelic et al . 35 are the first to propose a 
partition in the turbulent energy spectrum. Because of the nonequilibrium, the rate 
that energy enters the low wave number region, e p , does not equal to the energy 
transfer rate from low wave numbers to high wave numbers, e t . Therefore it is 
reasonable to describe the evolution of the energy contained in low wave number 
region, k p , and high wave number region k t) separately. As a result, the time scale 
or the length scale defined by different energy transfer rates will be different and 
this multiple-scale concept reflects the nonequilibrium effect of turbulence. 
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We think that this concept would be more appropriate for compressible flows 
because the compressibility often creates nonequilibrium interactions between large 
and small eddies. We first modify Hanjelic et al.’s model, test it in various free 
shear flows and boundary layer flows and then extend it to compressible flows by 
consideration of the effects of compressibility on the equations for k p and e p . The 
proposed model is tested in both compressible free shear flows and boundary layer 
flows. For detailed analysis and flow calculations see the report by Duncan et a/. 36 
and Liou and Shih 37 . 

The Researchers involved with the subject in this section are W. Liou, T.-H. Shih 
and B. Duncan. 

3.7 Bypass transition model 

The onset of turbulence transition in the propulsion system is often highly influ- 
enced by the free stream turbulence. This transition process does not go through 
the linear instability but is mainly controlled by nonlinear processes. Therefore, it 
is sometimes called “bypass” transition. Because of this highly nonlinear process of 
transition, turbulence models may be used to predict it. In fact, many two-equation 
models, for example, k-e eddy viscosity models of Launder and Sharma 38 , Chien 15 , 
etc., do mimic bypass transition on a flat plate when the free stream has a certain 
amount of turbulent intensity. However, to obtain an accurate prediction of bypass 
transition, the study of the bypass transition process and physics is needed. The 
conventional turbulence models must be modified to take into account the intermit- 
tent phenomena of transitional flows. 

We have proposed transition models based on a two-equation turbulence model 
using an intermittency factor to modify either the eddy viscosity or modeled k-e 
equations. Successful results for a flat plate boundary layer under various free- 
stream turbulence intensities are obtained. For details see the report by Yang and 
Shih 39 . 

The Researchers involved with the subject in this section are Z. Yang and T.-H. 
Shih. 

3.8 Joint scalar PDF model 

One of the critical problems in turbulent combustion is how to treat the inter- 
action between the chemical reaction on the turbulence. The estimation of the 
production rate of compositions based on the mean flow temperature would be in a 
very large error for flows with finite rate chemical reactions. The reason is that the 
production rate of compositions depends not only on the mean values of tempera- 
ture T and compositions Ci, but also very much depends on the detailed fluctuations 
of temperature 6 and compositions Cj. The moment closure scheme of modeling the 
production rate of compositions in terms of the mean flow temperature, the mean 
compositions and various correlations consisting of the fluctuating temperature and 
composition, such as 0 2 , Oa, cjc], ..., has not been successful. However, the PDF 
method allows us to treat chemical reaction exactly without modeling (Pope 40 ). 
Therefore, for the study of turbulent combustion problems, we use the joint scalar 
PDF transport equation for the scalar field and the moment closure schemes for 
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the velocity field and develop a hybrid solver consisting of a Monte Carlo scheme 
and a conventional CFD method. For detailed description of this procedure and its 
applications see Hsu 41 and Hsu et a /. 42 

The Researchers involved with the subject in this section are A. Hsu, A.T. Norris 
and J.Y. ChenT 

4. RNG and DIA 

In developing one point turbulence models, conventional modeling methods can 
be supplemented by “non-conventional” methods such as renormalization group 
theory (RNG) and the direct interaction approximation (DIA). These are two point 
theories formulated in wavevector or fourier space; one point models are derived 
by integration over wavevectors. This approach provides theoretical support for 
conventionally derived models and sometimes suggests theoretically derived forms 
for the empirical elements, whether constants or functions, which appear in these 
models. 

We have applied RNG methods to both the eddy viscosity and Reynolds stress 
transport equation models. In addition to the k — e model proposed by Yakhot 
and Orszag 43 , it is possible to obtain constitutive relations for Reynolds stress 
and heat fluxes (Rubinstein and Barton 4,6 which are special cases of the general 
results Eqs.(2.1.3) and (2.1.7). By applying the perturbation theories of Yakhot 
and Orszag 43 to the relevant correlations, expansions in powers of the mean velocity 
gradient are obtained for the stresses and heat fluxes; quadratic truncation of the 
series leads to a stress model Eq.(2.1.3) with constant 04 , 06,07 and a heat flux 
model Eq.(2.1.7) with constant 02,03 in which the constants are in good agreement 
with empirically selected values. The forms derived are also consistent with the 
DIA analysis of Yoshizawa 3,5 . 

The RNG method also provides a formulation for closing the Reynolds stress 
transport equation (Rubinstein and Barton 44 ). Perturbative evaluation of the cor- 
relations n'J and nfj leads to series expansions in powers of the mean velocity 
gradient. These series can be consolidated, or “resumed” using the known pertur- 
bation series for the Reynolds stresses by methods analogous to Pade approximation. 
Systematic lowest order summation leads to a Reynolds stress transport equation 
with a form identical to the LRR model equation and with constants in reason- 
able agreement with empirically chosen values. Higher order resummation leading 
to nonlinear models of the type described in Sec. 3.5 remains an open possibility. 
The possibility of such resummation in the context of DIA has been discussed by 
Yoshizawa 45,46 . 

Recent work has focussed on nonequilibrium time dependent relations between the 
Reynolds stress and the mean flow derived from a simplification of the DIA theory 
of shear turbulence. In this theory, shear turbulence is modeled by a non-Markovian 
eddy damping acting against the mean shear. The RNG and DIA Reynolds stress 
transport models and the LRR model all assume Markovian damping; as in the 

f Professor, University of California, Berkeley, CA 
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molecular theory of transport coefficients, Markovian damping describes long time 
behavior and is incorrect at short times. The most important consequence of non- 
Markovian damping is a strong suppression of eddy damping at short times. This 
leads to closer agreement between the present theory and rapid distortion theory 
at short time, this is important in modeling oscillating shear flows: recent work of 
Mankbadi 47 shows that RDT based models best predict such flows. In transient 
homogeneous shear flow at high strain rates, the LRR model predicts rapid onset 
of eddy damping leading to excessive growth of turbulence kinetic energy at short 
times. The suppression of eddy damping at short times in the present model should 
lead to improved predictions for this flow as well. 

Another consequence of this theory is a stress model Eq.( 2 . 1 . 3 ) in which the coef- 
ficients a2> • • • are functions of the mean strain rate. This theory can be described as 
RDT with a modified total strain determined by the response function of the DIA 
theory of isotropic turbulence. The introduction of a phenomenological modified to- 
tal strain has often been advocated in the RDT literature to improve the agreement 
between RDT and shear flow data; here the modified total strain is deduced as a 
consequence of the theory. In the special but important case of simple shear flow 
in which dUi/dxj = S6n6j2, the result can be formulated in terms of Eq.( 2 . 1 . 3 ) in 
which, for example, = a2(Sk/e) and the function 03 is found exactly from RDT. 
There are analogous results for the coefficients 04,0,6,0.7; in simple shear flow, the 
remaining terms in Eq.( 2 . 1 . 3 ) identically vanish. Extension of this theory to other 
mean shear tensors depends on the tabulation of the corresponding RDT solution. 

The researchers involved with the subject in this section are R. Rubinstein and 
A. Yoshizawa.t 

5. Numerical Simulation 

To obtain a better understanding of the effect of compressibility and rotation 
on turbulence, numerical simulations of compressible homogeneous shear flows and 
rotational flows are carried out. The effects of compressibility and rotation on the 
energy spectrum and energy cascade between turbulent eddies has been analyzed 
(Hsu and Shih 48 ). These simulations support the idea of the multiple-scale model 
for nonequilibrium compressible turbulent flows (W. Liou and Shih 37 ). 

Another numerical simulation is the transition subjected to the free stream large 
disturbances. The objective of this simulation is to obtain some insight into the 
transition physics and to provide data base for bypass transition modeling. Based 
on the assumption that the transition process is mainly controlled by large scale 
motions, we use a high accuracy finite difference Navior-Stokes solver with course 
grids to simulate the large scale motions of transition. A preliminary calculation of 
bypass transition was carried out. Various statistics of the calculated flow field are 
under examination. 

The Researchers involved with the subject in this section are A. Hsu, C. Liou^, 
Z. Yang, A. Shabbir, T.-H. Shih. 

f Professor, Tokyo University, Japan 
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Research on Two Equation Models 


Presented by Z. Yang 


NASA Lewis Workshop on Computational Turbulence Modeling 

September 15-16, 1993 


Motiva tions and Objectives; 

• k - e model is the most widely used turbulence model in engineering calculations. 

• However, the following deficiencies need to be fixed: 

- Currently, most k — e models for near wall turbulence contain geometry parameters. 

- The form of the e equation lacks a solid theoretical foundation. 

- The k — e model performs rather inadequately for flows with adverse pressure gradient. 

- The capability of k — e models in predicting bypass transition due to the freestream 
turbulence needs improving. 
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Modeling of Near Wall Turbulence 


• Near wall k — e model = Standard k — e model + near wall effect. 

• The near wall effect: 

- The time scale approaches the Kolmogorov time scale near the wall. 

— The damping function is parametrized by a new parameter which is independent of the 
coordinate system. 

• The resulting model is Galilean and tensorial invariant. 

• The resulting model is robust numerically. 


Turbulent Channel Flow at Re r = 395 
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ZPG Turbulent Boundary Layert at Reg — 1410 
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On the Wall Functions 


• The advantages of the wall function approach: 

- Reduce the number of grid points by half, at least. 

- Reduce the numerical stiffness of the dissipation rate equation, by less grid stretching. 

• The limit of wall function approach: the flow is assumed to be attached to the wall. 

• Existing wall functions are based on the flat plate BL at zero pressure gradient. 

It is inadequate when the pressure gradient is not zero. 

• A new set of wall functions axe obtained: 

- They axe based on the asymptotic behavior of the governing equations in the log layer. 

- They contain the effect of the pressure gradient. 


PPG Boundary Layer APG Boundary Layer 
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A Vorticity Dynamics Based Model for the Dissipation Rate Equation 


• The dynamic equation for the fluctuating vorticity is analyzed. 

(The terms in the fluctuating vorticity equation have clearer physical meanings than terms 
in the dissipation rate equation.) 

• For large Reynolds numbers, e = vUJiUJi. 

• The resulting model equation has a better foundation than the standard e equation. 

• The resulting model equation always gives a positive production in dissipation rate. 

- The model calculation is expected to be more robust for complex flow calculations. 


Flow Inhomogeneity and the Dissipation Rate Equation 

• The exact dissipation rate equation contains source terms due to the flow inhomogeneity. 

• However, the existing e equations are homogeneous. 

(The source terms are the same for both the homogeneous and inhomogeneous flows.) 

• A new model equation for e is proposed, which accounts for the inhomogeneity effect: 

- Flow inhomogeneity is represented by V5 and Vfc. 

- Invariant theory is used to derive such a model equation. 

• The resulting model equation accurately account for the effect of the pressure gradient. 
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Rotating homogeneous shear flows 


Q/S= 0.0 


Cl/S= 0.5 


a/5=- 0.5 



Mixing Layer 


Plane Jet 


Round Jet 
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FPG Boundary Layer 


FPG Boundary Layer 




Mild APG Boundary Layer 


Mild APG Boundary Layer 
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Moderate APG Boundary Layer 


Moderate APG Boundary Layer 




Strong APG Boundary Layer Strong APG Boundary Layer 
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Mode ling of Bypass Transition 


• Low Reynolds number k — e models could mimic transition. 

- The predictions are not very good. 

- Among these models, the Launder-Sharma model gives the best prediction. 

• New model for bypass transition is proposed by introducing the effect of the intermittency 

• The transition model recovers to turbulence model at the end of the transition zone. 

• Calculations of the benchmark flows show that the present model gives better predictions 
compared with the k — e models without the intermittency effect. 


Transitional Boundory Layer Transitional Boundary Layer 

Tu=3% Tu=6% 




Re.X Re_X 


35 


Su mmaries: 


The capabilities of k — e model are enhanced in the following areas: 

• A Galilean and tensorial invariant k — e model for near wall turbulence. 

• A new set of wall functions for attached flows. 

• A new model equation for the dissipation rate: 

- It has a better theoretical basis. 

- It contains the contribution of flow inhomogeneity. 

- It captures the effect of the pressure gradient accurately. 

• A better model for bypass transition due to freestream turbulence. 
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A REYNOLDS STRESS 
ALGEBRAIC EQUATION MODEL 


Presented by 

J. Zhu 


Motivation 

• Two-equation turbulence models (K-e, K-cv) 

- Simple, robust and 
computationally inexpensive 

- Isotropic 

- Not realizable, ujui may become negative 

• Reynolds stress transport equation models 

- Anisotropic 

- Capable of simulating iipj-transport 

- Difficult to model higher-order correlations 

- Complex and computationally expensive 
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Objective 


Develop a Reynolds stress algebraic equation model 
which combines the respective advantages of 
two-equation models and 
Reynolds stress transport equation models: 

• Simple, robust and efficient 

• Anisotropic 

• Realizable 


Derivation of Constitutive Relation 

• First step: assume some relationship between 
quantities of interest u^u] = Fij(Uij , K, e) 

• Second step: use invariant theory to determine 
the function 

1. An invariant can only be a function of other 
invariants 

2. If the LHS of an equation is bilinear in arbi- 
trary vectors A* and Bj, the RHS must also 
be bilinear in A* and Bj 
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• A general constitutive model of ~u t u 3 


O r/2 O 2 

u,u } — -fc<5y + ^ a 2 —^-(Uij + Uj t i — -Ui,i6ij) + 2 a 4 — —II' (!) t j ) 

+ 2 a 6 -^-(Ui'icUj t k — -n 2 <5ij) + 2a 7 — (Uk,iUk,j - glRtfy) 

+ 2ag (Ui % k Uj t k + - | n 3^ii) + 2a 10 ^-(U kyi Uli + U ^ U L ~ f^i) 

+ 2a l2 ^-(Ul k Ul k - -n 4 <5ij) + 2a 13 — (f/^C/fcj - 
+ 2a ki -^j-(Ui i kUi i kUij + U j<k U ltk U? ti — -n 5 5ij) 

Zf6 2 

+ 2axg~pr(Ui tk Ui k Ulj + U j<k Ul k Uli — -Ilefiij) 

+ 2.a l g~{U x<k Ui^ k U^ m U^ m + Uj tk Ui tk Ui m Ui iJn - -H-iSij) 

rii = Ui, k u k ,i , n 2 = Ui, k Ui , k , n 3 = Ui <k uf k , n 4 = u? k uf k , 
n s = Ui^Ui^ufi , rie = UijPfjPt* i iR = Ui tk Ui <k u l m u i m 


The first five terms are of the same form as that of Yoshizawa (1984) 
and Rubinstein and Barton (1990). 


Realizable Algebraic Equation Model 


U{Uj — — K6ij - v t (Uij 4- Uj t i) 


C T 1 


K ; 


A 2 + V J + £ 3 e 2 


( u i,kUk,j + Uj,kUk,i - o n<5 t?) 


A 2 - ■ - in««) 


A 2 + r, 3 + C 3 e 2 
C T 3 K 3 

'a 2 + V 3 + ^ 3 e 2 


(u kii u kJ - ±r\ 6 ij) 


K,t + [UjK - O + = -uiU]U ii:j - e 

e,t + [Uj€ - (v + -^r) e ,j],j — -Ci—uiUjUij - C 2 — 
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where 


2/3 K 2 
A 1 + V + «£ € 


I<S 


S = (2 SijSij) 1 ' 2 , 



e 


_ Kn 



e 


Cri = -4, C t2 = 13, C r3 = -2, Ai = 5.5, a = 0.1, A 2 = 1000. 
Ci = 1.44, C 2 = 1.92, a K = 1, cr e = 1.3 

Solution Procedure 

• Conservative finite-volume method 

• Non-staggered grids 

• Momentum interpolation for velocity-pressure 
coupling 

• Second-order accurate differencing schemes 

• Quadratic terms in u{i zj treated as source 

• Strongly implicit solution algorithm of Stone 
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Performance in Complex Turbulent Flows 


Diffuser flows 

Backward-facing step flows 
Confined jets 

Confined swirling coaxial jets 

Diffuser flows 

• Fraser’s case, a = 10° 

• Trupp, Azad & Kassab’s case, a = 8° 
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Wall pressure 



SKE : Standard K-e model 
RNG : RNG-based K-e model 
RRSAE : Present model 
KO : K-u model 


Skin friction 
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Uc/Uref 


Centerline velocity 



I 


Shape factor 
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Wall pressure 



x/Ro 


Skin friction 



x/Ro 
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Centerline velocity 


Trupp et ol.’s cose 



x/Ro 


Displacement thickness 


Trupp et ol.'s cose 
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Backward-facing step flows 


• Driver & Seegmiller’s case, H s /H d =1/8 

• Kim, Kline & Johnston’s case, H s /H d = 1/2 



Comparison of reattachment points 


case 

experiment 

K-e 

present model 

RSM 

ASM 

DS 

6.1 

4.99 

5.82 

- 

5.66 

KKJ 

7 ±0.5 

6.35 

7.35 

6.44 
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Wall pressure 




Reynolds stresses (KKJ-case) 
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Reynolds stresses 
(DS-case) 
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Confined jets 

• Barchilon & Curtet’s case, 6 = 0 ° 

• Binder Kian’s case, 9 = 2.5° 
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Wall pressure 
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Axial mean velocity profiles 




BK-case 


Turbulent shear stress (BK-case) 
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Qr/Q 


Separation and reattachment points (BC-case) 
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Confined swirling coaxial jets 


• Roback &. Johnson’s case 
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r/Ro r/Ro 


Mean velocity profiles (U, V, W) 
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Turbulent intensity profiles (V, v', w') 
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Conclusions 


• A general constitutive relation for u^u] has been 
derived 

1. Its quadratic form seems sufficient 

2. It is realizable 

• It significantly improves the performance of K-e 
based models 

• Easy to implement 

• No significant increase in computing time 

• Very robust 
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ASSESSMENT AND DEVELOPM ENT OK SECOND ORDER 
TURBULENCE MODELS 


Presented by 
Aamir Shabbir 


MOTIVATION 

• these models describe the effect of mean flow and external 
agencies (such as buoyancy) on the evolution of turbulence 

• therefore, in principle, these models give a more accurate de- 
scription of complicated flow fields than the two equation mod- 
els 

- e.g flows with large anisotropy in turbulence (such as near 
the leading edge of a turbine blade) 
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OBJECTIVE 


• assess the performance of the various second order turbulence 
models in benchmark flows 

• seek improvements where necessary 

- model for the pressure correlation term in the scalar flux 
equation 

- model for the scalar dissipation equation 


Transport Equations for Second Moments 


DuiUj 

Dt 


— p ij + + Tij — Dij 


Dutf 

Dt 

d!p 

lot 


—Pi + IIj + Ti 


=P + T- D 


These equations have to be closed by providing models for: 

• Pressure correlation terms (Iky, IT) 

• Transport (Diffusion) terms (Ty , T) 

• Dissipation terms (Ay, D) 
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HOW TO ASSESS M ODELS ? 


Global computation 

• Mean and turbulence equations are numerically solved 


Du{Uj 

Dt 


—Pij + 






• Results (e.g. Reynolds stresses) are then compared with ex- 
periments or DNS data 


Direct comparison 

Individual terms in the turbulence equations (such as pressure cor- 
relation terms) are directly compared with experiment or DNS data 
Note that: 

- In experiments pressure correlation terms can not be mea- 
sured but can only be obtained indirectly through balance 
of second moment equations 

- DNS allows direct computation of these correlations but 
is limited to low Reynolds number 

Most of the results to be shown in this presentation are direct com- 
parisons 
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Models for pressure correlation term in the scalar flux equation 


• a. 

• b. 


• c. 

• d. 


Launder model (1975) 

Zeman and Lumley model (1976) 

- linear in scalar flux 

- do not satisfy realizability 

Shih, Shabbir and Lumley model (1985,1991) 
Craft, Fu, Launder, Tselepidakis model (1989) 

- linear in scalar flux and Reynolds stress 

- satisfy realizability 


Application to Homogeneous Shear Flow 
(Experiment as well as DNS data) 
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Tavoularis & Corrsin Expt (1981) DNS Data of Rogers et al. (1986) 
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Shabbir and George Expt. (1990) 


Application to Round Buoyant Plume Flow Experiment 
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CONCLUSION 


Models for pressure correlation term in scalar Hux equation 

• Models involving both scalar flux and Reynolds stress give bet- 
ter performance than the models which involve only scalar flux. 


Models for pressure correlation term in the Reynolds stress equation 

• a. Launder, Reece, Rodi model (1975) 

• b. Speziale, Sarkar and Gatski model (1991) 

- linear or quasi-linear in Reynolds stress 

- do not satisfy realizability 

• c. Shih and Lumley model (1985) 

• d. Craft, Fu, Launder, Tselepidakis model (1989) 

- nonlinear in Reynolds stress 

- satisfy realizability 
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Tavoularis & Corrsin Expt (1981) Tavoularis & Corrsin Expt (1981) 
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DNS Data of Rogers et al. (1986) DNS Data of Rogers et al. (1986) 
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DNS Data of Rogers et al. ( 1 986) DNS Data ° f R og ers et al - ( ' 1 986 ) 
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CONCL USION 

Models for pressure correlation term in Reynolds stress equation 

• For the DNS data non-linear models give better performance 
than linear models. However, for the experiment no single 
model performs better for all the components 

• For the rapid part of the pressure correlation the relation 
ng = F(iiiUj , Uij) is found to be adequate 


CONCLUSION fcontd.) 

• Performance of all the slow pressure correlation models varies 
from one flow to another. 

• Furthermore, the relation = F(uiUj,k , e), is inadequate in 
certain situations 

- DNS shows that Ilg is dependent not only on the present 
time value of Reynolds stress but also on its past history 

- Definition of Ft?- implies that it is also a function of triple 
velocity moment U{UjUk 

• Therefore, more research is needed before any model for 
can be recommended for use. 
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A New Model Equation for Scalar Dissipation 

• Traditional scalar dissipation rate equation is modeled in an 
analogue fashion to the mechanical dissipation equation 

• Equation proposed here is modeled after the exact equation 
for scalar dissipation 

• Its production/destruction mechanisms are different than the 
traditional model equation 


Application to Homogeneous Benchmark Flows 

1. Homogeneous turbulence subjected to constant scalar gradient 

2. Homogeneous turbulence subjected to constant scalar gradient 
and constant shear 

Global computation of the following two equations 


Uj 

Uj 


de 2 

dxj 

dee 

dxj 



— 2e0 


—Coit-eS + Cq 2 


y/e e e$ 

V Pr 



Mechanical field (i.e. fc, e, etc.) and scalar flu x,Ui6, are taken as 
known. This way performance of the scalar dissipation equation is 
isolated. 
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Tavoularis & Corrsin Expt(1981) 
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CONCLUSION 


• The transport equation for thermal dissipation rate proposed 
here gives improvement over the standard equation in at-least 
all the simpler benchmark flows. Its performance in the wall 
bounded flows is being assessed 
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A Multiple-Scale Model for 
Compressible Turbulent Flows 


presented by 


William W. Liou 


Workshop on Computational September 15-16, 1993 

Turbulence Modeling 


Outline 

□ Motivation 

□ Model Development 

• formulation and determination of model coefficients 

□ Model Predictions: 

• compressible shear layer 

• compressible boundary layer 

□ Concluding Remarks 


69 


Motivation 


□ Incorporate the effects of non-equilibrium energy spectrum 
of compressible turbulence in computational models 


DNS of Compressible Turbulence 

□ formation of high gradient regions or eddy shocklets 

□ enhanced vortex stretching and spectral energy transfer 

□ non-equilibrium turbulent kinetic energy spectrum due to 
the flow compressibility 
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Multiple-Scale Model 



Model Equations 

□ Large Scale 



d r , pt \ dk p . ,du.~ 

[(P+r-)^rl + wfc) 2 


dy 


dy 


dy 


pe v + fci 



• fci - exchanges between the turbulent kinetic energy and internal energy 

• fc 2 - increased spectral energy transfer due to eddy-shocklets 
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Model Equations (cont.) 
□ Small Scale 




d r / — . ^T,dk t 

[(/* + — )-z-\+ P 6 p ~ P e t 
dy a~ dy 


hv+—)^]+ct 1 p^- ct2 A 


dy 
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dy 


□ Eddy Viscosity 


fj, T « pul ~ p{k v + k t \ 


i ( kp kt'j 


Model Coefficients 

□ Decaying turbulence and homogeneous turbulence 


Cpi 

Cti 


(1 - -) + -Cp 2 , Cp 2 

a a 


n + 1 


n 


x — — + Ph 
0 0 ’ 


Ct 2 = 


0 -i + Cptfh 

hp 

0 + 0 h-i 


C „ 

a k 

„ K p 

(Jr 

kt 

a 7 

G~ 

€ 

n 

q: 

P 


UP 

1.0 

1.3 

1T3 


T2 

TU5 


72 


Model Equations (cont.) 

□ Large Scale 

• fci - Sarkar et al. (1992) 

. du.o , 

fci = -0.2 + ocsM t pe p 

• fie 2 - dimensional analysis 

fc 2 = fc 2 (p,P,k,e,Jl) 


n _ P_ _ J_ TT ifl J_ 
1 ~~ pfc “ 7M t 2 ’ 2 

fc 2 = b\ + 6 2 -^ 2 + 0(.M^) 


Model Predictions 
□ Compressible Free Shear Layers 



• mean velocity profiles. M c =0.51 

• peak Reynolds shear stresses. M c = 0~1.6 

• vorticity thickness growth rates 
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Multiple-Scale Model 
□ Wall Function 
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Model Predictions 


□ Compressible Boundary Layers 



• mean velocity profiles. Moo=2.831, i?e#=420,700 

• mean temperature profiles. 

• skin friction coefficients. Moo=0~5, Reo= 10 4 
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Concluding Remarks 

□ The present compressible multiple-scale model predicts 
correctly both the spreading rate of compressible shear 
layers and the skin friction coefficient of compressible 
boundary layers. 

□ Need to implement the model into the calculation of more 
complex compressible turbulent flows. 
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PDF TURBULENCE MODELS 
FOR REACTIVE FLOWS 


Presented by 
A.T. Hsu 

OUTLINE 

1. Motivation: why PDF 

2. Works acomplished: 

Model development. 
Numerical results. 

3. Future studies 
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CLOSURE PROBLEM: 


ti,- = Wi + u' i , 

Yi = Yi + Y; t 

u"u" — Turbulence Modeling 

Y"u" — Analogy of shear stress: diffusion model. 

pWi — ??? 


pwi = pw { Y u ...,Y n ,T,p) 

But in general: 

pw] # pw { Y u ...,Y n T,p) 


EXAMPLE: REACTION RATE FROM MEAN TEMPERATURE 

CH4-tir Film* B: x/D-20 4-»t«p 


Co + Hs 0 ^ co,+H«. 
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Final Solution 


Fig. 4 Velocity distribution 
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CONCLUSIONS 
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ANALYTICAL THEORIES OF TURBULENCE 
APPLIED TO TURBULENCE MODELING 


presentation by 
Robert Rubinstein 
CMOTT 


MOTIVATION 

Turbulence models contain undetermined elements 

constants 
functions 
model forms 


Apply analytical theories to determine these elements: 

direct interaction approximation (DIA) (Kraichnan, 1959) 
renormalization group (RNG) (Yakhot and Orszag 1986) 
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EXAMPLE 1 


Determine constants in nonlinear eddy viscosity (NLEV) (Yoshizawa 1984) 
and nonlinear eddy diffusivity (Yoshizawa 1987) models 


< uu > -\kl = C„— (V U + VU T ) 

O 6 

+ C Tl ^(VUVU T - \ VU : VU T I) 

6 <5 

+ C r2 ^r(VUVU + VU T VU T - - VU : VU I) 

e z 3 

+ C t3 ^(VU T VU - i VU T : <5U I) 

e z 3 

< 9n > = C„ — V0 + C^VU • V0 + C k2 ^VU t • V0 

e e z 


Universal inertial range constants for isotropic turbulence: 
E(k) = C K e 2 ' 3 k ~ b ' 3 

i /(k) = C D e 1 / ' 3 k ~ 4/ ' 3 

Ck — 1.61 Cd — 0.49 (Yakhot and Orszag 1986) 

Assume shear is a weak perturbation of isotropic 
background state. Then (Rubinstein and Barton 1990): 

exact value experiment 

C v ( l) 2 C D /C 2 K 0.085 ~ 0.09 

C Tl 4(f) 2 /{C d C k ) 2 0.034 ~ 0.04 


c r3 toK \) 2 I{C d C k ) 2 -.014 


-.01 


94 


Passive scalar results 


Pr-p — f(Pr,Rc) (Yakliot and Orszag 1986) 

/ is a known function: for Re = oo, P-rp ~ 0.7 
= CJPr T 

Constants of the nonlinear theory are known functions of Cd^CkiP^t 
(Rubinstein and Barton 1991) 

DifFusivity ratio «i 2/«22 = 1-9 (Tavoularis and Corrsin), 2.3 (theory) 


In principle, higher order nonlinearities in NLEV 
can be computed theoretically 


EXAMPLE 2 

Determine model form for rapid pressure strain correlation II 

RNG gives a series for Id in powers of V77. 

The LRR model 


II = Cfi < uu > VU + C+ 2 < uu > VU T 

arises from systematic lowest order consolidation (resummation) of this 
series. (Rubinstein and Barton 1992): 



LRR 


RNG 

Crl 

0.7636 

16 

21 

~ 0.7619 

Cr 2 

0.1091 

2 

21 

~ 0.0952 
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Second order resummations are possible: 


n = F(< uu >)VU + G{< uu >)VU T 


n =< uu > [Cj VU + C 2 VUVU + • • •] 
These are speculative at this point. 


EXAMPLE 3 

Determine unknown functions in generalized nonlinear eddy vis- 
cosity representation for simple shear flow 

Nonlinear eddy viscosity does not apply at large strain rates: it predicts 
< uu > /k — ► oo. Introduce the generalization 

< uu > -\kl = CM — (VU + VU T ) + ••• 

o 6 

where in simple shear flow dUi/dxj = S6u8j2 , tj = Sk/e 


96 


DIA analysis of shear turbulence leads to 


< UiUj >= kFij(r)/CR ) 

where Fij is known from rapid distortion theory (RDT) (Maxey 1982) and 
Cr is a universal constant {Cr ~ 1.59, Yakhot and Drszag 1986). 

General time dependent theory 

< UiUj > ( t ) = kFij(a*(t)) 

where a;* is a modified total strain determined by DIA: 

a*(t) = [ ds G{t — s)lS(s) 

Jo 

G is the response function of DIA. 


Homogeneous shear flow at high strain rate 
r) = Sk/e — * oo 



6x2 

bn 

622 

633 

NLEV 

0 ( 77 ) 

Oirfi) 

0 (v 2 ) 

Oijfi) 

LRR 

o(tr l ) 

/11(c) 

/22(c) 

f 33(c) 

present 

Oiv- 1 ) 

2/3 

-1/3 

-1/3 

RDT 

o(v - 1 ) 

2/3 

-1/3 

-1/3 


In LRR, fij(c) are model dependent constants. The limiting normal stresses 
are not necessarily realizable. 

DNS (Lee, Kim and Moin) shows that RDT is a good description of these 
flows even for moderate times. 
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CURRENT WORK 


Time dependent model for nonequilibrium turbulence based on DIA 
New evaluation of NLEV constants 

Nonmarkovian damping: short time suppression of eddy damping. 

Required for oscillating and highly strained flows. 

Two scale theory of shear turbulence: 
large scales RDT, 

small scales ~ Kolmogorov inertial range dynamics 
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DNS OF TURBULENCE, TRANSITION 
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DNS for Bypass Transition: Numerical Aspects 

• Fully implicit, high-order finite difference scheme 
(fourth order in space and second order in time.) 

• Acceleration Techniques: 

l 

— Semi-coarsening multigrid 
— Line-box distributive relaxation 

• A new buffer outflow boundary condition 

• A fine-coarse-find grid mapping 
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ACCOMPLISHED 


DNS for Bypass Transition: Cases Calculated 


• 2D/3D linear instability stage (temporal simulation) 

• 2D/3D linear instability stage (spatial simulation) 

• Secondary instability stage (spatial simulation) 

• Simulation of the whole transition zone (spatial simulation) 


2 
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DNS: Effect of Rotation 
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EQUATIONS FOR TRIPLE CORELATION EQUATIONS 


o 

+ 



Em 

c 

4 - 




Em 

1 

in 


O 

• f — 1 

'C 


+ 

>~-o 

CP 

1-1 



• 

3 





co 




d 

d 


0 

O 


•+j 

• r— « 
■4-5 

J 

ji 

Ip 

E 

cT 

i~< 

i-t 

C 

'co 

co 

• r— < 



CP 

T5 

V4J 

CP 

T3 

CQ 


d 


1 

13, 

cO 



a 


1 

-40 


Jo 
fa c 

l-J 

/ v 


CP 


•fO 

<-+-! 

c 

d 

<p 



CO 

CP 


;5 

-P= 

O 

-<-5 

r->. 


»o 
• CO 

CP 

fed 

J3 

• 0 



w 



O 

z 

w 

p 

PQ 

p 

H 

o 

z 

I— I 

% 

o 


-54 

£ 

+ 


• ro 

-f 


co 

r~< 

c 

• r— < 

CO 

r-t 
r— ' 

C7 1 

CP 

co 

co 

CP 

l-l 

-t-5 

co 

CO 

13 

d 

>> 

CP 

K 

cp 

H 


o, 

5p G 

+ J 

o E 

55, CM 


G 


-54 CM 

3 I 


-34_ 

-34 


+ 


JO 

bO 

CP 

d 

<p 

CP 

• f— < 

-4-5 

<p 

d 


d 

CP 


-O 

d 

H 


=3 

bi | CO 


'<~2 

*«m 


a 

i I CM 


co 

co 


CP 


co 

CO 

<p 

(-1 


-34. 

-34' 


C3 

CM 


C'l 


CM ICO 
*54* I 


-34 


.s* + 


-54 


CP CM 


O 

<P 


C'l 

p 

CM 

I 


103 


0.040 


*S 



ssb* is spiouAay 


II 

C 


8 

to 

8 

u 



O 

>> 

U 

o 


V 


e 

H 



03 

o 


to 

d 


o 


\ 


6 

P 



104 


.11 One- dimensional energy spectra for the flow field at 
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o o o.2 o * o c o a i o 
Fie." The ratio of vorlicity square, representing incom- 




DNS: Reduced H2 Chemistry 
Temperature [k] 




! 


I 
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DNS: 


12 Chemistry 
Taction 



Yno 


0.000137 
0.000129 
0.000120 
0.0001 1 1 
0.000102 
9.358E-5 
• 471E-5 
7.584E-5 
8.697E-5 
5.810E 5 
4.923E-5 
4.037E-5 
3.1S0E-S 
2.263E-5 
1.376E-5 


7.5 
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National Aeronautics and 
Space Administration 

Lewis Research Center 


MODELING OF BYPASS TRANSITION 
Fred Simon 



• Presentation of the Problem 

Variables 

Meaning of Bypass Transition 
Physics/Mechanisms to be Modeled 

• Modeling of Transition Onset 

Use of DNS 
Use of k-e Model 
Use of PSE 
Summary of Results 

• Mo d eli ng o f t he Trans iti on Re gion 

Comment on Linear-Combination Method 
Use and Limits of k-e Methods 
Schmidt/Patankar Method 
Multi-Time-Scale Method 
Intermittency Method 
Reynolds Analogy 
Summary of Results 

» Recomme nd ations 
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Effect of Streamwise Curvature on Bypass. 
(Wang, 1984; Kim and Simon, 1991) 


Stanton 

number 


Transition Wall Radii of 



CD-93-66129 
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Linear Versus Bypass Path to Transition, From Suder, 
O'Brien, Reshotko (1988) 

Uoo = 100 ft/sec 



TRANSITION 


TU = 0.3% 


Re^lOOO 


I 1402 



1693 

[ . 


1898 


2096 

2294 

0 10 20 30 40 50 

Tima, ms 


TU = 1% Re * /100 ° 

209.1 
308.3 

407.5 
506.9 

606.5 

Junto} 7 ° 6 -° 

805-2 

804.7 

0 10 20 30 40 50 

Time, ms co**4«us 




Computed Skin Friction Along the Length of the 
Flat Plate (Grids A and B), (Rai and Moin, 1991) 


8x1 0 -3 □ Computation (Grid A) 
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Comparison of Computed Disturbance Energy from 
Computations Based on a k-e Model and on Direct 
Numerical Simulation (DNS) for Bypass Transition 
(Simoneau and Simon, 1993) 



CD-93-66128 


Onset Prediction Using PSE Approach 
(Stuckert & Herbert, 1992) 

Data of Sohn & Reshotko(1991) 


GrldO 
Grid 0.5 


© Grid 1 

* Grid 2 

<> Grid 3 
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Computed and Experimental Momentum Thickness 
Reynolds Number for transition Onset 



Tu, percent 

CD-93-66127 


/ 

s 





S3 
i * 
£8 
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The Numerical Prediction of Transition Flow Using 
MTS Model Compared with Experimental Data of 
T3A (Crawford, 1993) 
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Intermittency (Volino and Simon, 1991) 
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Kim Tu = 0.3 (3D) 

Kim 1.5 (3D) 

Kuan 0.9 (2D) 

Sohn 1.1 (ID) 

Sohn 2.6 (1 D) 

Blair 0.93 (3D), k = 0.2x10-6 
Blair 2.0 (3D), k = 0:2x1 0-6 
Blair 1.9 (3D), k = 0.75x10-6 
Blair 5.3 (3D), k = 0.75x1 0-6 
Narasimha 


1.5 
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Transition Length Determination 
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Variation of Nondimensional Spot Formation Rate 
with Free-Stream Turbulence (Narasimha(1 985)method) 


/✓- 
/J.- 
/•o 
0.8 - 
o 6 - 
O'/ - 
- 


VoW 


□ Schbauer and Skramstad ( 1948 ) 

O Abu-Ghannan and Shaw (1980) 

A Sohn and Reshotko(1991) 

D Suder, O' Brian and Reshotko(1988) 
A Kim ( 1990 ) 
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Use of Intermittency to Model Transition Region 




Without intermittency, Tu = 2.8% 


Use of an Intermittency Weighting Factor in Computations 
(Yang, 1992; Yang and Shih, 1992) 
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National Aeronautics and 

Space Administration MODELING OF BYPASS TRANSITION 

Lewis Research Center Fred SimOh 


REYNOLDS ANALOGY EFFECTS 


Problem: Reynolds analogy breakdown in flows with pressure 

gradient 


Possible Solution: Use of additional modeled transport equations for the 

thermal turbulence to predict the turbulent heat flux/ 
turbulent Prandtl number 



Computed and Experimental Momentum Thickness 
Reynolds Number for Transition End 
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(1) Volino, 
Simon 
(1991) 



A Kim (1990) 

O Sohn (1 991 ) 

□ Blair and Anderson 
(1 987) k = 0.2x1 O' 6 
C> Blair and Anderson 
(1 987) k = 0.75x1 Or® 


B MTS (2) Crawford 1 

■ TXM (2) Crawford 

• Jones-Launder (1) Simon, Stephens 

(1991) 

A Launder-Sharma'l (2) Sieger, Schulz, 

■ Nagano-Tagawa r Crawford, Wittig 

T K.Y. Chlen J (1992) J 


Experiments 


Computations 


(1) Based on intermittency 

(2) Based on minimum Stanton number 

Abu-Ghannam and Shaw 
(1980) 
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National Aeronautics and 

Space Administration MODELING OF BYPASS TRANSITION 

Lewis Research Center Fred Simon 

RECOMMENDATIONS 

• Develop a DNS data base with/without heat transfer for turbjne blade 
geometry for (1) understanding physics (2) supporting/guiding expe- 
riments (3) developing/testing models 

• Continue development of models more faithful to physics 

• Continue development of Multi-Time-Scale Model 

• Develop Reynolds Stress Turbulence Models 

• Initiate an effort in the use of Large Eddy Simulation(LES) 

• Develop the PSE approach as design tool 

• Develop models for unsteady flows (e.g. wakes) 

• Develop capability to predict turbulent Prandtl number 

• Test models for turbulence levels found in combustors (6-20%) 
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Title 


1 


ALGEBRAIC MODELS FOR TURBINE 
BLADE HEAT TRANSFER 


R. J. Boyle 

NASA Lewis Research Center 


Workshop on Computational Turbulence Modeling 


LeRC 


9/93 


Aim and Scope 


• VALIDATE PREDICTIONS OF TURBINE 
BLADE HEAT TRANSFER & PRESSURE LOSS 

SMOOTH AND ROUGH BLADES 

• ALGEBRAIC - BALDWIN-LOMAX TYPE 
MODELS 

• IMPLEMENTED IN N A VIER -STOKES CODES - 
RVC3D & RVCQ3D 

• ENGINEERING APPLICATIONS 


Workshop on Computational Turbulence Modeling 
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LeRC 


9/93 SMOOTH SURFACE PREDICTION 


• FREESTREAM TURBULENCE - ex LEADING 
EDGE 

• TRANSITION MODEL - ex MAYLE'S 

• ALGEBRAIC - BALDWIN-LOMAX TYPE 
MODELS 


Workshop on Computational Turbulence Modeling 



120 


LeRC 


9/93 


BLADE HEAT TRANSFER 




Workshop on Computational Turbulence Modeling 


LeRC 

9/93 


6 


APPROACHES FOR ROUGH SURFACES 


• TAYLOR, COLEMAN, and HODGE 

EXPLICITLY ACCOUNT FOR BLOCKAGE, 
DRAG AND HEAT TRANSFER 

SOURCE TERMS IN NAVIER-STOKES 
EQUATIONS 

• CEBEC1 and CHANG - USED IN ANALYSIS 

MODIFY MIXING LENGTH 

ONLY EDDY VISCOSITY AFFECTED 

Workshop on Computational Turbulence Modeling 
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LeRC 


9/93 


REQUIRED ROUGHNESS PARAMETERS 


• HEIGHT 


• SPATIAL DENSITY 


• SHAPE 


Workshop on Computational Turbulence Modeling 


LeRC 


9/93 


8 


CORRELATION FOR EQUIVALENT HEIGHT RATIO 
FROM SIGAL and DANBERG 



1.0 10.0 100.0 1000.0 


ROUGHNESS DENSITY PARAMETER, A 


Workshop on Computational Turbulence Modeling 
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LeRC 


9/93 


EXPERIMENTAL DATA USED FOR COMPARISONS 


• HEAT TRANSFER 

HOSNI - FLAT PLATE 

BLAIR & ANDERSON - LARGE LOW SPEED 
ROTOR 

DUNN & KIM - SSME FUEL TURBINE 


• TURBINE EFFICIENCY 

BOYNTON - TWO-STAGE SSME HIGH PRESSURE 
FUEL TURBINE 


Workshop on Computational Turbulence Modeling 


LeRC 


9/93 JlO 

ROUGH FLAT PLATE HEAT TRANSFER 



Workshop on Computational Turbulence Modeling 
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LeRC 


9/93 


1 1 


RATIO OF PREDICTED TO MEASURED HEAT TRANSFER 

FOR A FLAT PLATE. Re x =400,000. 

(DATA OF HOSNI et al.) 



Equivalent Height Correlation 


Sigal &3Danberg 

Dvorak 

Dirling 

L/D 

h e q/h 

stp/st m 

heq/h 

Stp/St m 

heq/h 

Stp/St m 

2 

1.58 

1.64 

2.43 

1.85 

1.69 

1.67 

4 

0.20 

1.30 

0.50 

1.57 

0.45 

1.52 

10 

0.023 

1.01 

0.001 

1.15 

0.079 

1.22 

Smooth 


1.10 


1.10 


1.10 


Workshop on Computational Turbulence Modeling 


LeRC 


9/93 


12 


PREDICTED AND MEASURED BLADE HEAT TRANSFER 



SURFACE DISTANCE, s/b x 


Workshop on Computational Turbulence Modeling 
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LeRC 

9/93 


13 


EFFICIENCY FOR SMOOTH AND ROUGH BLADES 


pr 



NORMALIZED WHEEL SPEED, 


Workshop on Computational Turbulence Modeling 


LeRC 


9/93 


GRID & MODEL ASSUMPTION EFFECTS 


14 


Grid 

y + 

K 

A + 

e 

Arj 

291x54 

0.3 

0.0 

26.0 

0.0321 




0.0 

1.0 

0.0483 

0.011 



0.3 

26.0 

0.0348 

0.002 



0.3 

1.0 

0.0550 

0.015 



1.0 

26.0 

0.0455 

0.009 



1.0 

1.0 

0.0703 

0.026 

291x54 

0.6 

0.0 

26.0 

0.0423 




0.0 

1.0 

0.0570 

0.010 



0.3 

26.0 

0.0449 

0.002 



0.3 

1.0 

0.0654 

0.016 



1.0 

26.0 

0.0546 

0.008 



1.0 

1.0 

0.0767 

0.023 

291x54 

1.2 

0.0 

26.0 

0.0375 




0.0 

1.0 

0.0531 

0.011 



0.3 

26.0 

0.0402 

0.002 



0.3 

1.0 

0.0586 

0.014 



1.0 

26.0 

0.0502 

0.009 



1.0 

1.0 

0.0730 

0.024 

Experimental 

0.021 


Workshop on Computational Turbulence Modeling 


125 


LeRC 
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ERROR DUE TO SURFACE RESISTANCE 


15 


• TEMPERATURE MEASUREMENTS BENEATH ROUGHNESS 
AFFECT HEAT TRANSFER COEFFICIENT 


^ACTUAL / ^MEASURED _ 1 + ht/k 

h - HEAT TRANSFER COEFFICIENT 

t - ROUGHNESS THICKNESS 

k - THERMAL CONDUCTIVITY OF ROUGHNESS 

• POSSIBLE REASON FOR HIGHEST AUGMENTATION 
CORRESPONDING TO LARGEST SCALE 


Workshop on Computational Turbulence Modeling 


LeRC 


9/93 


16 


REQUIREMENTS FOR BENCHMARK ROUGHNESS 
HEAT TRANSFER DATA SETS 

• COMPLETE DESCRIPTION OF THE ROUGHNESS 
CHARACTERISTICS 

• SURFACE NOT DISTURBED BY MEASUREMENT 
TECHNIQUE 

• KNOWN TEMPERATURE OF ROUGHNESS ELEMENTS 

• VERIFICATION OF ACCURACY BY MEASUREMENTS ON 
SMOOTH SURFACES 


Workshop on Computational Turbulence Modeling 
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LeRC 

9/93 

CONCLUSIONS 

17 


• SMOOTH SURFACES 

REASONABLE ACCURACY 
TRANSITION BEHAVIOR GOVERNS 
ACCURACY 

• ROUGH SURFACES 

INCONSISTENT ASSUMPTIONS 
NEEDED FOR LOSS AND HEAT 
TRANSFER 


Workshop on Computational Turbulence Modeling 
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NAVIER-STOKES TURBINE HEAT TRANSFER 
PREDICTIONS USING 

TWO-EQUATION TURBULENCE CLOSURES 
AL1 A. AMERI 

CENTER FOR RESEARCH, INC. 
UNIVERSITY OF KANSAS 
NASA LEWIS RESEARCH CENTER 



TURBI NE HEAT TRANSFER 

OUTLINE 

Brief description of the method of solution 

2- d examples using various models 

3- d calculations using algebraic model and q - co model 
Conclusions 




129 




TURBINE HEAT TRANSFER 


TURBULENCE MODELS 


LOW REYNOLDS NUMBER TWO-EQUATION MODELS 

• C OAKLEY’S q - co MODEL 

• CHIEN'S K-e MODEL 

ALGEBRAIC MODEL 

• BALDWIN -LOMAX - WITHOUT TRANSITION MODEL 



TURBINE HEAT TRANSFER 


FORMULATION 


• MASS AVERAGED COMPRESSIBLE N-S EQUATIONS 

• + MODEL EQUATIONS 

• NON-PERIODIC C GRID. 
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TURBINE HEAT TRANSFER 


SOLUTION OF THE NAV1ER-ST0KES EQUATIONS 


• FINITE VOLUME DISCRETIZATION 

• 4 STAGE RUNGE-KUTTA SCHEME 

• EIGENVALUE SCALING OF ARTIFICIAL DISSIPATION 

• VARIABLE COEFFICIENT RESIDUAL SMOOTHING 

• MULTIGRIDING 



TURBINE HEAT TRANSFER 


SOLUTION OF THE MODEL EQUATIONS 


• FINITE VOLUME DISCRETIZATION 

• 4 STAGE RUNGE-KUTTA SCHEME 

• VARIABLE COEFFICIENT RESIDUAL SMOOTHING 

• MULTIGRIDING 
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2-D CASES 

• SSME HIGH-PRESSURE FUEL TURBINE 1ST STAGE VANE 

• ALLISON’S C3X VANE DATA OF HYLTON ET AL. 

| 3-D CASES 

• LANGSTON’S CASCADE, THICK & THIN INLET BOUNDARY 
LAYERS. 


TURBINE HEAT TRANSFER 




Surface Pressure Distribution- First Stage Stator 
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TURBINE HEAT TRANSFER 



- jr 


Suction Surface 


1 — 

0 . 


n r" 


50 . 


100 - 


S/C 


Grid Spacing from the wall 
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TURBINE HEAT TRANSFER 



First Stage Vane, Low Reynolds Number 



TURBINE HEAT TRANSFER 



S/C 

First Stator, Variation with length scale 
< 7-00 Model 
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TURBINE HEAT TRANSFER 



S/C 


First Stator, Variation with length scale 
Chien's Model 


C3x 149 
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Thick Inlet Boundary Condition 




q-omega model, 193x49x65 grid 


Surface Stanton Number x 1000. , thick inlet boundary layer 


' MM«CI -J riuiuml HlMMl' 
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Langston’s Cascade 
Endwall Stanton Number xIOOO. 




Experiment, Graziani et al. Computation, TRAF3D 

Thick inlet boundary layer. 


THICK BOUNDARY LA YER INLET CONDITION 




-0.4 -0.1 0 0.2 0.4 0.6 08 1.0 1.2 1.4 |. 6 


q-omega Model, 193x49x65 Experiment Graziani 
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Surface Stanton Number x 1000. , thin inlet boundary layer 



Experiment, Graziani et al. 



Computation, TRAF3D 


SURFACE STANTON NUMIJEltx 1000, THIN INLET ROUNDARY LAYER 


I** ■ — WCUJI »«■». I .1. Miluai Iin» VI 




(j-d) model, IWxdOxftS urid 
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Langston’s Cascade 
Endwall Stanton Number xIOOO. 



Experiment, Graziani et al. Computation, TRAF3D 

Thin inlet boundary layer. 


LANCS' TON'S CASCADE 
ENDWALL. STANTON Nl IMIlERx l()(M) 




Experiment, Graziani cl. al. 


<1 in model 193x49x65 grid 
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Internal Fluid Mechanics Division 


?■ twit fitGoffipanlst $ 


TURBINE HEAT TRANSFER 


CONCLUSIONS 


• Two-equation models perform well in the fully turbulent regime. 

• The uncertainty re. length scale ofturb. needs to be resolved. 

• Transition modeling is crucial to the success ofH-T analysis. 

• Need further improvements in convergence speed. 
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9/93 


Objective and Rationale 


Objective: Improved prediction of turbine blade surface 
temperature and heat-transfer distributions 


Rationale: Existing heat-flux models lack the fundamental 
generality required for complex flows 
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Kasagi, Tomita, and Kuroda DNS 

S 


Chennai Flo* Re, . 1 to 


Channel Flow Re, • 1 (0 





— i rnTTm — i i 1 1 ii in i nruin — i i min 
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LeRC 9/93 Nagano-Tagawa with Pr t = 0.86 vs. Kasagi DNS 


7 




good overall agreement and proper asymptote for <uv> h 


LeRC 

9/93 

Nagano-Tagawa with Pr t = 

= 0.86 vs. Kasagi DNS 

8 


Chann«( Flow A« ( * ISO 


Clwnnal Flow R« f • I (0 




ini nn i mi Mill i i iii in 




o • 


- 

10' 

I / 1 - 


e • 

* 

t 


/• \ 

10-* 

1 / ' 


0 4 


: \ 

10* 

* / - 


03 



10'* 



1 

■' to' io* io* • or 1 io' io' to* t 

r »' 

mediocre overall agreement, but proper asymptote for<v0> + 

o’ 
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Two-equation thermal models 


°i 


Hattori, Nagano, and Tagawa (HNT) 

complex damping functions (evaluation in progress) 


Sommer, So, and Lai (SSL) 

Shima coincidence condition (evaluation in progress) 


Thermal time scale model (development in progess) 
simple wall boundary condition: t 0 =o 
simple near-wall asymptote: t 0 ~ y 2 / 2 a 
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A COUPLED IMPLICIT SOLUTION METHOD FOR 
TURBULENT SPRAY COMBUSTION IN PROPULSION SYSTEMS 


K.-H. CHEN 

The University of Toledo, Ohio 
and 


J.-S. SHUEN 

Sverdrup Technology, Inc., Ohio 


Aerothermochemistry Branch 
Internal Fluid Mechanics Division 
NASA Lewis Research Center, Cleveland, Ohio 


OBJECTIVES 


• Develop an efficient and robust algorithm for multi-phase 
chemically reacting flows at all speeds, with emphasis on 
low Mach number flows. 

• Calculate turbulent spray combustion flow in a gas tur- 
bine combustor. 
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MOTIVATION 


• Many reacting flows in propulsion devices cannot be efficiently calculated by 
modern compressible flow CFD algorithms, e.g., 

rocket motor — wide range of Mach numbers, from near zero 
velocity at closed end to supersonic at nozzle exit, 
gas turbine combustor — low subsonic velocity, but large density 
variation precludes incompressible approach. 

• Most low-speed reacting flow codes based on TEACH-type technologies - 
inefficient and lack of robustness for complex flows. 

• Tremendous progress made in high-speed compressible flow CFD in past 
two decades. Extending application range to low-speed flow regime highly 
desirable. 


OUTLINES 


• GOVERNING EQUATIONS 

- Gas-Phase Equations 

- Liquid-Phase Equations 

• NUMERICAL ALGORITHM 

. NUMERICAL TEST RESULTS 

• CONCLUSION 

• FUTURE PLAN FOR ALLSPD CODE 
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GOVERNING EQUATIONS 


• Gas-Phase Equations 
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ALLSPD MATN FEATURES 


Present Capabilities 

• 2-D and Axisymmetric geometries 

• Second-order central difference for both inviscid and viscous terms 

• Fully coupled, fully implicit algorithm 

• Efficient convergence for wide range of Mach numbers (from M < 10 -10 to 
supersonic) 

• Finite-rate chemistry, realistic thermophysical properties 

• Multi-block and body-fitted curvilinear coordinates for complex geometries 

• ac — e turbulence model 

• Stochastic liquid spray model (dilute spray), vortex model for droplet internal 
circulation and diffusion 

• Still a research code, require experience and knowledge of CFD and flow 
physics to use. 

Future Plans 

• Development of a more efficient solver 

• Extension to 3-D 1 

• PDF model for turbulence/chemical reaction closure 

• Thermal radiation model 

• Detailed soot and NOx kinetic models 

• Multi-grid and unstructured grid capabilities 

• Dense spray and high pressure (near-/super-critical) spray models 


Difficulties with Compressible Flow Algorithms at Low Mach 
Numbers 


• Disparities among system’s eigenvalues (stiffness), u, u + c, u — c, re- 
sulting in significant slowdown in convergence rate. 


• Singular behavior of pressure gradient term in momentum equations 
as Mach number approaches zero, 


p 4 u * 2 + 


P 4 

7 M? 


As Mach number is decreased, pressure variation (Ap* oc M 2 ) becomes 
of similar magnitude as roundoff error of the large pressure gradient 
term (p* / 7 M 2 ). 
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METHOD OF APPROACH 


Pressure Singularity Problem 

• Pressure decomposed into two parts: 


P = Po + Pg 

p g replaces p in momentum equations and retains p 9 as one of the 
unknowns. 

• Employs conservative form of governing equations, but uses primitive 
variables 

( P g t u,v,h,Yi) 

as unknowns. Conservation property preserved and pressure field ac- 
curately resolved for all Mach numbers. 

Eigenvalue Stiffness Problem 

• Pressure rescaled so that all eigenvalues have the same order of mag- 
nitude. Physical acoustic waves removed and replaced with pseudo- 
acoustic waves which travel at speed comparable to fluid convective 
velocity. 


EIGENVALUES RESCALING 


\ = u,u, i[[/(i + £) ± sju'{ 1 - ly + 4P( Q ] + *!)], u, u ■ ■ ■ 

— £x> a 2 = { v , U = a t u + ajv. 

For well-conditioned eigenvalues, scaling factor P taken to be 

P = u J + v 7 . 
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LIQUID-PHASE EQUATIONS 


• Droplet Motion Equations 


du p 

Hi 

dv p 

Hi 


dx p 

Hi 




3 Hfip 

16 p r r\ 

3 C D l^g 
16 p r rl 


K “ u p). 
(”j - v p)> 


( 2 ) 

(3) 

(4) 

(5) 


He p = 2^[(u,-u p ) 2 + (o t -v p )Y /2 , 


24 Re 2 / 3 

-ii(I + ^-) for Rt p < 1000, 
Rt p 6 

C D = 

0.44 for Re p > 1000. 


• Droplet Heat and Mass Transfer Equations 


^ = 2N,ln{\ + B), (6) 

P Dj 

hd p 2N p ln(\ + B) u ~' 

k [(1 + B)^-' - 1] K ’ 


N. = \ + 


0 MGRe'J'Pr'l 3 
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l 1 + RcpPr'l* 1 
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0 .216Re l J*Sc'/ 3 
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1 — 1^/jp 


Y„ p 
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X hp W, + (1 - X, gp )W a 
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• Droplet Internal Temperature Equations 
(Vortex Model) 




(8) 


r(t) _ L(9m\ r 

c[t > ~ 17 k, ’ ” dt 


t = t 


tnj > 


T v = T in 


a = 0 , 


yl dT P 

da 17 v k, ,rp dt 


dT p 1 , C p ipu a ~-i- 

o — i»V l. > r P zu ’ 


o = l. 


3a 16 r ” 3r 


NUMERICAL ALGORITHM 


• Gas-Phase - ALLSPD code 

• Liquid-Phase 

— Droplet motion equations (ODE) - Runge-Kutta method. 

- Droplet internal equations (PDE) - implicit method 
(Thomas algorithm). 

— Determination of spray time step for integration. 

— Stochastic separate flow model. 

• Interaction Between Two Phases 
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SPRAY TIME STEP 


• Droplet Velocity Relaxation time (i r ) 

tr = T ( f )( 7 )( Cd 

3 P 9 v 

• Droplet Life Time (<;) 


</ = 


3 rri'pi 


• Droplet Surface Temperature Constraint Time ( t ,) 


U = fn ( 1 1 &t p )/ A ' 

1 A 




• Local Grid Time Scale (t s ) 

• Turbulent Eddy-Droplet Interaction Time (t,) 

t{ = t e , if L e > t |u" -Up| 
i; = min(t e ,t t ), if L c < t\u - u p | 

where 

L. = C 3 J\ 3 '*/e 
t e = L c /{2k/ 3Y'\ 
t, = — r/n[l - LJ(t\u - u p '|)] 

• Spray Time step - A T apr 

Ai, P r = a min (t r ,t,,t„t g ,ti) 
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INTERACTION BETWEEN TWO PHASES 


1. Initialize gas and liquid phase variables. 

2. Solve liquid-phase equations. 

3. Evaluate spray source term, H/. 

4. Solve gas-phase equations and update gas-phase 
variables. 

5. Update spray source term, H j ? 

No, go to step 4. 

Yes, go to step 2. 


NUMERICAL TEST RESULTS 


• Turbulent Backward-Facing Step Flow - non-reacting. 

• Evaporating Turbulent Spray Flow. 

• Gas Turbine Spray Combustion Flow. 
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OPERATING MODE 1 

NO ATOMIZING AIR 

LIQUID FLOW RATH : 1-26 grams /sec 
AIR FLOW RATH : 0.00 grams / sec 
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Gas phase u velocity Gos phase u velocity 




I 
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Fuel vapor mass /lux fuel vapor mass flux 

1 0.30 , 


o 



lujou z~\ 



Sjtu/6»i ’ *ui s t ui/6>i *Jui 
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Axial position. 
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Fig. 9 Center-piane cut for GE’s EEE gas turbine enr.ne and .be grid for the corabus ^ Velocity vectors (colored by temperature) for (a) non-combustion and (b) spray com- 

bustion cases for the gas turbine spray combustion Bow. 


luromc spray combustion Qo 


1.2 norm 



CONCLUSION 


• ALLSPD code efficiently coupled with the SSF spray 
model. 

• Satisfactory convergence property for flows with and 
without spray. 


162 


Numerical Investigation of Complex, Transitional, and 
Chemically Reacting Flows 


S.-W. Kim 

Resident Research Associate 
NASA Lewis Research Center 
Cleveland, Ohio 44135 


Part A. On Turbulent Transport of chemical Species in 
Compressible Reacting Flows 

Contents 

Mixing of chemical species in reacting flows 
Analysis 

Chemically reacting, turbulent flow equations 

Density-weighted, time-averaged Navier-Stokes equations 
Multiple-time-scale turbulence equations 
Species conservation equations for reacting flows 

Numerical Method 

Mixing and combustion of hydrogen in vitiated supersonic airstream. 

Comparison with measured data and other numerical results obtained using 
k-e turbulence models and a PDF. 

Conclusions and Discussion 
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Mixing of chemical species in reacting flows 


Chemically Reacting Laminar Flows 

Numerical (or semi-analytical) methods 

O.D.E. solvers 

1-D Numerical methods 

Numerical uncertainty is minimal 

Chemical kinetics 

Chemical kinetics (finite rate chemical kinetics & Reduced finite rate chemical 
kinetics) for certain fuels (hydrogen and some of hydro-carbons) have been 
tested and validated repeatedly. 

Numerical and theoretical analyses yield accurate results. 


Chemically Reacting Turbulent Rows 


Numerical methods 

Boundary layer or Navier-Stokes equations solvers 
(Except for semi-analytical analyses of PSR cases) 

Uncertainty is caused partly by numerical methods. 

Chemical kinetics 

Finite rate chemical kinetics & Reduced finite rate chemical are used. 

Uncertainty is caused partly by over-simplified chemical kinetics. 

Turbulence Equations 

Numerical analysis yields not so accurate results. Unceratinlty caused by turbulence 
equations is by far greater than that caused by chemical kinetics.* 

The probability density functions are used to improve the predictive capability of 
turbulent mixing of chemical species. However, some pdf methods yield physically 
incorrect numerical results. 

It is certainly important to better understand the turbulent mixing of chemical 
species in reacting flows. 


•Westbrook and Dryer, Prog. Energy Combustion Science, vol. 10, pp. 1-57, 1984. 
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Turbulent Mixing of Chemical Species in Reacting Flows 

1. Large eddy mixing 

2. Turbulent Mixing 

3. Molecular Mixing 

The use of accurate chemical kinetics alone can not yield accurate results for 
numerical analysis of turbulent reacting flows unless turbulent mixing is resolved 
correctly. 


Large Eddy Mixing 

Caused by separated flows, recirculating flows, and organized structures. 

The extent of calculated recirculation zone or the organized structures depends 
on the accuracy of both the numerical method and the turbulence equations since 
they are nonlinearly coupled with each other. 



Therefore, the capability to correctly resolve the large eddy mixing depends on 
the accuracy of the numerical method and the turbulence equations. 

Numerical investigations carried out during last decades show that k-e, ARSM, 
and RSM do not yield accurate results for complex turbulent flows. 

On the other hand, the multiple-time-scale turbulence equations* yield highly 
improved numerical results for various complex turbulent flows. 

*Kim and Benson, Int. J. Heat Mass Transfer, vol.35, pp. 2357-2365, 1992. 
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Turbulent Mixing 


Caused mostly by energy-containing eddies and partly by fine-scale eddies. 

The capability to resolve the turbulent mixing of chemical species depends on the 
capability of the turbulence equations to correctly describe the turbulent 
transport of scalar variables (i.e., heat transfer, turbulent kinetic energy, 
concentrations, and convection-diffusion of species, etc) 

Need to be able to describe the chemical reaction-turbulence interaction (i.e., 
turbulent mixing is enhanced and shear layer thickness is widened by chemical 
reaction)* 

The multiple-time-scale turbulence equations can resolve the cascade of 
turbulent kinetic and the nonequilibrium turbulence phenomena. 

(Present numerical results show that the M-S turbulence equations can resolve 
the chemical reaction-turbulence interaction.) 


* C.T. Chang et al., "Comparison of reacting and non-reacting shear layers at a high 
subsonic mach number," AIAA paper 93-2381, 1993. 


Molecular Mixing 
Caused by molecular diffusivity. 

Theoretical and numerical analyses of chemically reacting laminar flows yield 
accurate results. 

Various molecular diffusion equations can accurately describe the molecular 
diffusion of species. 

The Lennard-Jones 12-6 potential law is used in the present study. 

Various numerical methods that yield accurate results for laminar flows can 
resolve the molecular mixing. 
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Multiple-time-scale turbulence equations 
Energy containing eddies 
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Turbulent Eddy Viscosity 

Pt=pCj I k 2 /e p where c p = c^ f et/e p 

Remark: Single-time-scale turbulence models can not resolve nonequilibrium 
turbulence phenomena. 


Species conservation equations for reacting flows 
Chemical species concentration equation 

^(pYij + V- {pYi(v + Vj)} = w, 

where the diffusion velocity, V,, is approximated using the Fick's law given as 
YjV, = - (Dj_j + D i t ) VYi 

The production rate of the i-th species, Wj, is given as 
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Chemical reactions for the combustion of H 2 in a vitiated supersonic airstream 
are described using 9 chemical species (H 2 , 0 2 , H 2 0, OH, O, H, H0 2 , H 2 0 2 , and 
N 2 ) and 24 pairs of reaction-steps (Burks and Oran, 1981; Kumar, 1989). 


A fast chemistry can not be used to describe the fine details of chemically 
reacting flows. 

A reduced chemical kinetics can not be used confidently due to the uncertainty 
contained in the reaction mechanisms. 

The use of a detailed finite rate chemistry may make it difficult to obtain a fully 
converged solution due to the coupling between the large number of flow, 
turbulence, and chemical equations. The numerical method needs to be 
strongly convergent. Accuracy also depends on the capability of turbulence 
equations used. 


The numerical method is a finite volume method that incorporates a 
pressure-staggered mesh and an incremental pressure equation for the 
conservation of mass. 

Predictor Step: Solve momentum equation. 


Numerical Method 


(pCi+A* )u**=I a£u*% S*- |f + PC 2 UP' 1 - pC 3 uf- 2 (1) 
nb ox i 


Corrector Step: Correct the velocity field to be divergence free 
Incremental pressure equation 



Incremental velocity equation 



(pCi+Ai) 9xj 


(3) 
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Velocity and pressure corrections 


(4) 

(5) 


U*“ = u** + u'i 
p“ = p* + p’ 

Solve eqs. (1-5) iteratively until all flow variables are converged. 



Combustion of H 2 in vitiated supersonic airstream (Burrows and Kurkov, 1973) 
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Total temperature 




Calculated OH contour. 


The numerical result obtained using the M-S turbulence equations indicates 
that the combustion occurs at approximately middle of the channel. The 
calculated flame location is in excellent agreement with that observed in the 
experiment. 
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Temperature contours obtained using the M-S turbulence equations 


It can be seen in the temperature contours that the calculated flame location 
obtained using the M-S turbulence equations are in correct agreement with the 
measured data. The contours also indicate that the temperature increases 
within in a short distance. The trend is in correct agreement with experimental 
observations that temperature increase occurs within a finite flame thickness. 



Temperature contours obtained using the pdf* 


The pdf fails to predict the correct flame location (i.e., ignition delay). The 
slowly increasing temperature field indicates that the pdf may not be able to 
predict a correct flame front. The numerical results obtained using the pdf are 
not in correct agreement with the physics of combustion. 


* A. T. Hsu, Y.-L.P. Tsai, and M. S. Raju, "A PDF approach for compressible 
turbulent reacting flows," AIAA Paper 93-0087, 1993. 
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Conclusions and Discussion 


The calculated species concentration profiles are in as good agreement with 
the measured data as those obtained using the pdf. 

The flame location (ignition delay) obtained using the multiple-time-scale 
turb ulence equations is in excellent agreement with the experimentally 
observed onset of ultraviolet radiation. 

Cascade of the turbulence filed is influenced by the extra strains caused by 
chemical reaction. 

Both the numerical results and the measured data exhibit enhanced mixing 
of the hydrogen and vitiated airstream for the reacting case. 

The M-S turbulence equations can resolve the chemical reaction-turbulence 
interaction. 

The pdf produces slowly increasing temperature field and it fails to predict 
the ignition delay. Thus the numerical results obtained using the pdf are not in 
correct agreement with physics of combustion. 


Part B. Unsteady Transitional Flows over Forced Oscillatory Surfaces* 

Contents 

Nomenclature 

Unsteady turbulent flow equations for flows with moving boundaries 

Navier-Stokes equations defined on Lagrangian-Eulerian coordinates 
Multiple-time-scale turbulence equations 

Numerical results 

Unsteady transitional flow field and comparison with measured data. 


* S.-W. Kim, K.B.M.Q. Zaman, and J. Panda, "Calculation of unsteady transitional 
flow over oscillating airfoil" in Separated Flows, eds. J.C. Dutton and L.P. Purtell, 
Proceeding of ASME Fluid Engineering Conference, Washington D.C., June 20-24, 
1993. 
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Nomenclature 



Lagrangian-Eulerian coordinates 

x: fixed reference coordinates 
moving coordinates 
J = la^j/aajl 


Unsteady transitional flow equations with moving boundaries 


Conservation of mass equation 



Conservation of linear momentum equation 



Convection-diffusion equation for scalar variables (i.e., <|> = (kp, £p, kp £p etc. 



+ J 


fie 


9<t> 

3x; 


- Jpf(4») 
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Numerical Method 


The unsteady transitional flow equations are solved using the same finite 
volume method. Time-integration is made using an iterative-time-advancing 
scheme. 

Comparison of Unsteady Flow Solution Techniques* 

Iterative Time-Advancing Scheme (ITA). 

Simplified Marker and Cell (SMAC). 

Pressure-Implicit Splitting of Operators (PISO). 

The ITA that can best resolve the nonlinearity of the Navier-Stokes equations, 
and yields the most accurate results. 

The SMAC is the most efficient computationally and yields accurate numerical 
results for laminar flows. 

The PISO is the most unstable numericaly and yields less accurate results. 


* Kim & Benson, Computers and Fluids, vol. 21, pp. 435-454, 1992 



Oscillating airfoil and moving mesh 
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Ensenble-averaged velocity profiles at x/ c=1.0 
The deteriorated comparison at a ~19°d 

(i) The hot wire can not accurately measure the velocity components when the 
flow is misaligned more than approximately 30° from the hot wire axis. 

(ii) The interaction between the DSV and the TEV occurs in a relatively coarse 
mesh region and the numerical method yields somewhat deteriorated results 
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Turbulent viscosity contour ( A(x/p. t = 10) 



CONCLUSIONS AND DISCUSSION 

Numerical method successfully predicts the Dynamic Stall Vortex and the 
Trailing Edge Vortex. 

The calculated ensemble-averaged velocity profiles are in good agreement 
with the measured data. 

Both the numerical results and the measured data show that the transition 
from laminar to turbulent state and relaminarization occur widely in space and 
in time. 

The good comparison between the numerical results and the experimental 
data are attributed to the capability of 

(i) the ITA that can best resolve the nonlinearity of the Navier-Stokes 
equations, 

(ii) the new pressure correction algorithm that can strongly enforce the 
conservation of mass, and 

(iii) the Multiple-time-scale turbulence equations that can resolve the 
transitional nonequilibrium turbulence field. 
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2c aj for u — u and c = Ju 2 + v 7 . 2d a, for u = « and c = \/u 2 + v 



Fig. 1 Dissipation characteristics of the CD6-RK4 scheme for CFL=0.6. 
la |g» . 2 (p, <7)| for u = v and c = \/u 7 + v 7 . 
lb lga(p» ^)1 for u = v and c = \/u 2 + u 2 . 
lc |p 4 (p, 9)| for u = v and c = \Ju 7 + v 2 . 
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SHOCKED SOUND WAVE 





185 


A LAMB VORTEX IN AN UNIFORM FLOW 
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00 l 
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Fig. 1 1 The vorticity magnitude contours for the vortex pairing. 
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, S The spatial growth rate as a function of the angular frequency from linear Fig. ;/ The vorticity magnitude contours for the vortex pairing, 

theory for the compressible free shear layer 
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waves in a compressible free shear 
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Fig. 12 Reconstructed eigenfunction components u' and ])' from the CFD results fo: 
the free shear layer perturbed by the most unstable mode. 





Fig. 14 Contours of constant vorticity of the simulated free shear layer perturbed by 
the most unstable mode at a magnitude of one hundredth. 
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Turbulent Back-Facing Step 
Flow and the k-e model: 

A Critical Comparison 



ref: /UAA-?3-l?2-? 


Christopher J. Steffen, Jr. 

Center for Modeling of Turbulence and Transition (CMOTT) 
NASA Lewis Research Center Cleveland, Ohio 


-■ — 

Lewis Research Center CFD Branch 


INTRODUCTION 

CMOTT at NASA Lewis: 

to address propulsion related turbulence issues 

□ assessment of state-of-the-art models 

□ development of new models 

□ validation of new techniques 



Objective of this work: 

to assess the performance of several low Reynolds number k-e 
formulations compared to the standard high Reynolds number 
form for separated flow over a step 
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Need to clarify this issue for code developers 
Critical comparisons between k-e models for: 

□ Channel flow: Lang & Shih (1991) 

□ Flow past a hill: Michelassi & Shih (1991) 

Avva (1990) compared the high Re t and Chien k-e models for 
several flows: noted the deficiency of Chien for step flow skin 
friction results . . . 

Contradicts idea of “more work<~>better result” 

Is this a property of all Low Reynolds Number k-e models? 


where the production of k is: V= -(uiUj)|(Uij + Uj,i). 
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k-e TURBULENCE MODEL 


Specifically, v t — C ^ 



Model transport of k and e: 



Primary difference between formulations: 
definition of / M , /i , / 2 , D, E 


Lewis Research Center 


CFD Branch 


4S& 

Low Re t forms: 

□ numerical domain extends to wall 

□ near wall profiles of U, k, and e are resolved not assumed 

□ typically first cell located y + « 1.0 

□ damping functions (fj, { 2 , and f^) and additional terms (D 
and E) tuned to specific flows 

Generality suffers somewhat . . . 

Formulation of above terms can affect the generality as well: 

□ JL and LS are functions of dependent variables alone 

□ CH and SL are functions of both dependent and independent 
variables 
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BACK-FACING STEP (BFS) FLOW 

Experiment of Driver and Seegmiller (AIAA J., ’85) 

\7////////////////////////////////////////////77A 
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streamline 
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Configuration: 

□ Re st ep of 33420 

□ inlet Mach Number of 0.128 

□ inlet tunnel 80H 

□ exit tunnel 60H 
Chosen for: 

□ LDV data for velocity, turbulent stresses 

□ wall static C p and wall C f 

□ reattachment length (time averaged) 


Lewis Research Center CFD Branch 


NUMERICAL EXPERIMENT 

Code: modified version of DTNS2D (Gorski ’88) 

Method: 

□ pseudo-compressibility (Chorin ’67) 

□ upwind differenced, TVD scheme for convective terms 
(Chakravarthy & Osher ’85) 

□ approximate factorization for time integration 

□ multiblock configuration 

□ decoupled (lagged) treatment of turbulence equations 

□ implicit treatment of turbulent source terms 
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Velocity profiles: 


o Driver and Seegmiller 

Jones and Launder 

Launder and Sharma 

Chlen 

Shih and Lumley 

Standard High Re, 
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Wall static pressure coefficient: 


o Driver and Seegmiller, stepside wall 
a Driver and Seegmiller, opposite wall 

Jones and Launder 

Launder and Sharma 

Chlen 

Shih and Lumley 

Standard High Re, 
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Turbulent kinetic energy profiles: 


o Driver and Seegmiller 

Jones and Launder 

Launder and Sharma 

- - - Chien 

Shih and Lumley 

Standard High Re, 
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Turbulent shear stress profiles: 


o Driver and Seegmlller 

Jones and Launder 

Launder and Sharma 

Chien 

Shlh and Lumley 

Standard High Re, 
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Skin friction coefficient: 
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Eddy viscosity field, uc 
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Primary reattachment length , \ T - 


model 

JL 

LS 

CH 

SL 

HR 

ja: r 

4.9 

5.4 

5.4 

5.1 

5.5 


Driver & Seegmiller recorded: 6.26H 


Secondary reattachment lenpth. Y r - 


[model 

JL 

LS 

CH 

SL 

1*1-2 

0.0 

0.0 

6.8 

0.1 


Driver & Seegmiller recorded: 0.5 < x r2 <1.8 
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CLOSING COMMENTS 

□ JL, LS, CH, SL and HR k-e models tested for back-step flow 

□ In general, results very similar between 5 different forms 

□ All do reasonably well for U and C p , albeit LR slightly 
better 

□ Turbulence profiles possess misplaced peak values as well 
as overpredictions in recirculation zone, common to all 

□ Reattachment length and Cf are difficult to properly resolve 
(14% error) especially for LR forms 

□ Definition of secondary recirculation zone inconclusive on 
this grid 

□ Perhaps a sensitivity to pressure gradient would help the 
LR damping functions . . . 
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Analysis of Supersonic Flows using k-e Model and the 
RPLUS code; Progress towards High Speed Combustor 

Analysis. 


by J. Lee 

Sverdrup Technology lnc./CFD Branch 
for Workshop on Computational Turbulence Modeling 
Sept. 1993 
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Outline 


• Problem of Interest - High Speed Combustor Flow Fields 

Parameters need to be Resolved 
Key Problems of Interest 

• k-s Model and RPLUS code 

Numerical Technique 
Models being Tested 
Some Results ' 

• Summary 



Problem of Interest 


• Analysis of Chemically Reacting flow inside of Supersonic RAM jet 
Combustors-Two Key Parameters need to be determined. 

Mixing/Combustion Efficiency 
Kinetic Energy Efficiency (Flow Losses) 

Inlet, Diffuser, etc.. 

• In order to do get some ideas on those parameter following (Potential 
Loss Mechanisms) must be modeled/determined correctly. 

Mixing, Shear, 

Turbulence, Vorticity, 

Shock-waves, Heat Transfer, 

Fuel Injector Drag, Poor Wall Pressure Integral, 

Chemical Dissociation. 

from 2nd JANNAF workshop on SCRAMjet Combustor performance workshop 
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Mixing and Injector Design 

• At High Mach Number(M ~ 5.0 +). 


Doesn’t mix well! 

The Natural diffusion mechanism very INEFFECTIVE. 

Fuel Residence Time Extremely Small- Even with Fast Fuel Such as H 2 
• Geometrical Complexities \ 

To induce Favorable mixing and Flame holding features 

Back-Step/Stream Wise Vorticity/Shock-Wave Interactions 
Unsteady Mechanism also being Envisioned as mixing enhancement 
Kumar, Bushnell and Hussani(1987) 



Introduction of Externally Generated Mixing Enhancements 


• Some External helping hand needed => Modeling Difficulties. 

• Externally Generated Vorticity Through Sweep angle of the Ramp 
injector. 

Davis(1990), Riggins and McClinton(1990), Drummond(1991). 

• Multiple Transverse Injection. 

Hartfield et. al. (1991) 

• Flame holding tricks/ Back-step with Recirculation. 

Hartfield et. al.(1991) 

• Simplified analysis of these features very difficult because of limited 
database/understanding (Attempts are being made using CFD 
solutions- JANNAF Combustor Subcommittee). 


s 
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Numerical Modeling(CFD) of Combustor Flow Field 


• CFD Analysis. 

• Numerical Modeling=> Overall Analysis of performance => Difficult 

• Overall Laminar Flow Fields with Complex Geometry /Finite Rate 
Chemistry has been demonstrated. 

• Finite Rate Chemistry Model^Yoon and Shuen(1989) 

• Multiple Grid Blocks- Moon (1991) 

• Analysis of a typical Injector Configuration with Zero Equation 
Turbulence Model using LU Scheme(RPLUS) code- Lee(1993) 



Simple Zero Equation IXirbulence model with multiple wall scaling 
Buleev-Inverse square rule can be used to extend model in to three- 
dimensional form. (Lee (1993)) 


• Good News/Bad News 

• Typical velocity profiles can be reasonably predicted. 

• Over all combustor flow features can be reasonably predicted. 

• Near-wall temperature characteristics near non-equilibrium region 
around the injector and separated flow were poorly predicted. 

• Overall spreading behavior of shear region poorly predicted. 

• Two Equation Transport Turbulence Model has the potential to ease 
some of these difficulties. 
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Sverdrup 

Technology Inc. 




Internal Fluid Mechanics Division 


Lewis Research Center 


Station I 



Normohzcd V velocity 



Normalized U velocity 



y/D 

Normalized Pressure 




Sverdrup 

Technology Inc. 


Internal Fluid Mechanics Division 


Lewis Research Center 


HARTFIELD ET. AL. (1990) 



THREE-BLOCK GRID SYSTEM 


TEST SECTION GEOMETRY FL0W SCHEMATIC 
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Sverdrup 

Technology Inc. 


IB , r, fSliRSA 

• Internal Fluid Mechanics Division 





X/H = 1 .06 
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CM 



?0 .1 .2 .3 .4 .5 .6 .7 

CONVECTIVE AM CH NUMBER 

GROWTH RATE VS. COMPRESSIBILITY EFFECT 

ye. n. lal 



Two Equations Transport Turbulence model are being Analyzed 


High speed turbulence models are some what Deficient (The 
deficiencies are well documented(Marvin(1986), Wilcox(1993)). 

Effect of Compressibility 
An-isotropy (Low/High Speed). 

Non-Equilibrium Flow Features (Low/High Speed). 

Near- Wall Flow(Low-ReynoIds Number Features (Low/High Speed)). 
Inflexibility of handling Complex Geometry- Invariance Principle 
(Low/High Speed) 

Large Dependence in the Numerical Methods Used 
(especially elliptic Solvers). 

Appropriate Initial/Boundary Conditions 
Etc ... 


209 





K-£ Model-RPLUS Development 


LU Based k-£ Model Solver-De-coupled Approach. 
Mean-Turbulence Transport Equations 

LU-SSOR- Yoon and Shuen- Explicit Terms Centrally Differenced 
LU-SW -Steger and Warming- Explicit Terms Upwind Differenced 

k-£ Models 

Convective Terms + Diffusive Terms + Source Terms = 0.0 
Model Only differ in Low-Reynolds Number Character. 

Models performance are being Evaluated. 

Implicit Source Term Handling Strategy also Being Studied 



k-£ Tbrbulence Models being studied for potential used in Three 
Dimensional RPLUS Code. 


• Low-Reynolds Number Model plus Dilatational Terms 

Chien (1976) 

Launder-Shima( 1976) 

Shih(1990) 

Various CMOTT derivatives of k-E Model 
Realizability 
Invariance 

Simplified Boundary-Conditions 

• Performance of the Low-Reynolds number K-e model in low-Mach 
number flows have been demonstrated (Patel, Rodi and 
Scheuerer(1985), Steffen(1993), Launder(1992)). 

• Some of the Potential Difficulties in high speed turbulence model are 
well documented (Marvin(1993), Coakley and Huang(1992). 
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Internal Fluid Mechanics Division 


Evaluation and Development of the RPLUS/k-e Model Solver 

Various 2D-3D problems are being studied to optimize the numerical 
method and to Evaluate model performance in supersonic flows in 
context to the LU based numerical Technique. 

Simple 2D k-e models are also being used to study various 
components of the flowfield generated by the complex combustor 
geometry previously shown. 

Studying the Numerical method/Model Behavior/Model Performance. 
2D Supersonic Turbulent Boundary-Layer- Skin Fraction/Heat transfer 
(NASA Ames Database). 

2D Supersonic Shock- Wave Boundary-Layer Interaction- Skin fraction/Heat 
-Transfer/ Shock-wave(A. Smits (1990’s)) 

2D Shear-Layer - Mixing (H. Lai(1993)) 

3D Fin/Flat Plate Interaction- 3D Comer Flows-Interaction Developed through 
a Fin generated Shock- Waves. (D. Davis(1992)) 
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(y-yO)/b 


1.5 


Supersonic Mixing Layer 

RPLUS vs. Dutton (case 2) 


- RPLUS 
C Dutton, case 2 


w 

o 

B 

o 






Boundary Layer 
Mach 2.87 
Re/m = 6.3x1 0 7 /m 


Law of Wall Profile 
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Sverdrup 

Technology Inc. Internal Fluid Mechanics Division 


Turbulent Shock- Wave/Boundary Interactions 
Mach 2.87 

Ramp Angle = 8.0 degrees 



Skin Friction Coefficient 
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Other Factors 

Optimum Numerical Strategy with in LU frame work. 

Effects of Initial condition. 

Modeling of Compressibility terms/Dilatational terms. 

Modeling of Turbulent terms in the Finite Rate Chemistry Model. 
Anisotropy of Turbulence 

Effects Upstream and Down stream Influences (Inlet(K. Kapoor) and 
Diffuser(?)). 

Chemistry-Turbulence Model Interactions (A. Hsu-PDF). 

Numerical Robustness(A. Suresh). 
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Validation of a k-e Model in RPLUS2D Code 
for Non-reacting/reacting Subsonic Shear layers 


H.T.Lai 

Workshop on Computational Turbulence Modeling 
September 15-16,1993 


OVERVIEW 

- To simulate an expirement at NASA Lewis by 
Marek and co-workers for subsonic shear layers 

- Computations for nonreacting/reacting flows 

- Grid refinement study for nonreacting case 

- Compressibility effects in k-e models (nonreacting) 

- Torch effects to sustain combustion (reacting) 
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The RPLUS Code 


- Navier-Stokes and species conservation equations 

- Finite-volume scheme 

- Options for central or upwind differences 

- Jameson-type dissipation 

- LU-ADI algorithm 

- Vectorization on oblique grid lines/planes 

- Hydrogen-air chemistry 

- Implicit chemical source terms 

- High Reynolds number k-e model 

- k-e equations uncoupled from other equations - 


NONREACTING SHEAR LAYER 


U=140 m/s, T=304.3°K, P=1.013xl0 5 N/m 5 

y =5 .08cm 

ti=2.25%, AIR 

symmetiy 

U=390 m/s, T=822° K, P=l.l013xl0 5 N/rn 
6=2.25%, AIR 

y=-5.08cm 


x=55.9cm 


GEOMETRY AND FLOW CONDITIONS 
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NONREACTING SHEAR LAYERS 

o UEASURED 

RPLUS 

WEC 



PROFILES OF AXIAL MEAN VELOCITY 



PROFILES OF TOTAL RMS FLUCTUATION 


NONREACTING SHEAR LAYER 



M=.679 


MACH NUMBER CONTOURS 




DISSIPATION RATE CONTOURS 
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Y, mm 



PROFILES OF DISSIPATION RATE 



SHEAR LAYER THICKESS VARIATION 
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CONVECTIVE hi A CH NUMBER 


GROWTH RATE VS. COMPRESSIBILITY EFFECT 
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REACTING SHEAR LAYER 


U=135m/s, T=366.5°K, P=l.013xl0 s N/m 2 
ti=5%, 3.94% H 2 and 96.06% N 2 

U=385m/s, T=810.9°K, P=1.0l3xl0 5 N/m 2 
ti=6%, AIR 

FLOW CONDITIONS 


dimensions lcm 

hydrogen torch t i=5%,T=L250° K, 1.3% H 2 + 21% 0 2 + 77.7% N 2 

V=0 (parallel inflow) 

PARALLEL INFLOW 


dimension=2. 12cm 

hydrogen torch ti=l 1%, T=1250°K, 0.1% H 2 + 21% 0 2 + 78.9% N 2 

V=50m/s (oblique inflow) 


OBLIQUE INFLOW 


REACTING SHEAR LAYERS 

o MEASURED 

RPLUS(OBUqUE) 

RPLUS( PARALLEL) 

FEC-PDf 



PROFILES OF AXIAL MEAN VELOCITY 



PROFILES OF TOTAL RMS FLUCTUATION 
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REACTING SHEAR LAYERS 
TEMPERATURE CONTOURS 




REACTING SHEAR LAYERS 
H 2 0 MASS FRACTION CONTOURS 



PARALLEL INFLOW 



OBLIQUE INFLOW 
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CONCLUSIONS 

- Growth rate underpredicted, especially reacting 

- Large decay of k near shear layer edges 

- Small production of k or large dissipation rate e 

- But center peaks in k are constant and overpredicted 

- Small improvement with grid refinement 

- Inhibition of growth with compressibility effect 

- Significant torch effects for reacting flows 
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NASA Lewis Research Center 
Cleveland, Ohio 
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• Modelled time-dependent sound source 


■n«p |«iuamuad*3 (446l) i,»JOoy»< 

’ " "Olinwrl"" *"»■“» I'”-’*"! ■»> 


Calculated Spectra of Sound Intensity in Decibels Referred to 
10' 12 W m' 2 Due to Coherent Structures 
at Various Emission Angles for U Q = 1 95 m s _1 



Sf 

CD-92-60388 
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The Study Indicates: 

II. Large-Eddy Simulations For Prediction of The Sound Source in a Sup ^ Thj _ | arge . sca | e structure seems to be the dominant sound source. 

• Results are sensitive to approximation in the sound source. 
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Harmonic excitation 
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(I) Extend Navier-Stokes computational domain to the far field: 
(a) Prohibitive storage requirements 
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(3) Finite-difTerence approach — Linearized Euler equations 
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(4) KirclihofPs method 
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JET NOISE PREDICTION USING A k-e TURBULENCE 

MODEL 


A. Khavaran 

Sverdrup Technology, Inc. 
NASA Lewis Research Center 


Lewis Internal Workshop 
on 

Computational Turbulence Modeling 
September 15-16, 1993 


MODELING APPROACH 


• Source Spectrum Calculations 

- Acoustic Analogy 

- Ribner and Batchelor Assumptions 

- Calculation of the Source Spectrum and its Characteristic Frequency Based 
on Time-Averaged Flow Calculation ( with a k-e Turbulence Model) 


• Sound/Flow Interaction 

- High Frequency Asymptotic Solution to Lilley's Eq. for Multipole Sources 
Convecting in an Axisymmetric Parallel Flow ( Balsa & Mani) 
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GOVERNING EQUATIONS 


Lighthill’s Equation 


dt* 


c'V'p = 


cPTjj 

dx-idxj 


Tij — pV\ Vj + Sfj (p — C p) €ij 


l dV L ,dV i _l f , SVu, 

e,J ^dxj + dx , 3 6,J ax t ) 

_iil 0(— ), fie = — 

pVtVj 'fie' p 

idp-dp = (^Ms 


» The effects of source convection and refraction are included in the source term 


Lilley’s Equation 


£) £> 2 ff 9 ,3cr dV, 9 2 

Dt Dt 2 9x( C 9x,' ^ dxidxj C 9x,- 


^314^ C 

9r, ax, c»x t DvDVCj, Dt 


+ viscous terms 


<r = -In— 
7 Vo 


_d _ _a a 

Dt ~ dt + **9ifc 


The effects of source convection and refraction are included in the operator term of 
Lilley’s eq. 


Source Spectrum Calculations 


- Mean-square sound pressure autocorrelation in the far field due to a finite 
volume of turbulence (in absence of convection and fluid shielding) 


P 2 (R,9,<t>) 


RiRjRuRt 

167 r 2 C*R G 


J J~(pVU){ P >VlVl)dtdy, 

Observer 




iu,T dr 


Fourth-order velocity correlation tensor 


/ +oo 

(v i vj)(v' k v' l )dt 

•OO 


Source strength (Quasi-incompressible turbulence) 

d 4 


, f d* 

•ijkt — P Sijudt 


Reduction in order of correlation tensor 


Sijkt — StkSji + StfSjk SijSkt 
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- Separable second-order tensors 


Sij(T,() = Rh(()G{t) 


- Isotropic turbulence model of Batchelor 

Rij(0 = Te-*«' L 'f x {[1 - tr(e/t t ) 2 ]«u + 

T = ^m, (’-tf + g+g 


- Gaussian correlation time delay 


G(r) = e -< r ' r ->’ 


- Source spectrum component 

Am (n)~p , **(nr.) 4 e =C ^ 
r * 3/2 Lr . 1 

T -~ 7 r fc = 2 ^ 


- Characteristic time delay of correlation 


To 


1 

(dU/dr) 


or T 0 


k 

t 


- Doppler shifted frequency 

n = 2 ) r/C, C = i/(l - M c cos6) 2 + (a c fc - 5 /^) 2 

M c = ,5W + p c Mj 
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Sound /Flow Interaction 


Mean square pressure in the far field 

p*{R, 0, Q) = f A(a*x + 4a IV + 2a vs + 2 a yt )dy 

Jv 


- Source term 


(AttRCooC) 2 {1 - Mcose)\l - M c cosef 


- Shielding function 


g\r) = 


(1 - Mcosefi^j 2 - cos'd 

(1 - MrCOsO)" 


M(r) = U(r)/C x M c — Uc/Coo 




Design Pnrmnters for G-D No/./.le 


Throat Diameter 0.1 in. 

Exit Diameter 5.395 in. 

Distance from Throat to Exit 5.525 in. 
Exit Velocity 2409 fps 

Ainhienf Velocity 400 ips 

Pressure Ratio 3.121 

Stagnation Temp. 171G"R 






“0/a 


Turb. Intensity 


Comparison of the Turbulent Intensity Profiles With Data 


Lip-llne 
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Strouhal Number 



OASPL DIRECTIVITY (ARC=40 ft) 




Design Parameters for C-D Nozzle 



Exit Diameter 5.395 In. 

Throat Diameter 5.1 In. 

Exit Velocity 2409 fps 

Pressure Ratio 3.121 

Stagnation Temperature 1716 R 


NOISE SPECTRA ( CD Nozzle) 




Freq. Band No. 



Band 151a 125 Ha 
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CHID 

COMlOUMin PLUG NOl/ll 



CHID 




MUCH NUrtDCR 


CPU |PHHC-k»I 


C OH TOUR LCUCLS 
0.20000 
0.2SOOO 
0.30000 
0.3SOOO 
0. 40000 
0 . 45000 
0.50000 
o.ssooo 
0 . 60000 
0. 6SOOO 
0.70000 
0. 75000 
n . 80000 

II . U50UU 
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I . 45000 
I . 50000 
I . 55000 
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TURD. INTENSITY ( Pf1RC-k«l 


CONTOUR LCUCLS 
0.00300 


o.ooo nrcn 
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0 00/00 
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0 03300 
0 . 03500 
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PARC-ke Data 





60 75 90 105 

Angle From Jet Axis 


120 140 160 


Design Parameters for CPN 

Exit Inner Diameter .597 in. 

Exit Outer Diameter 1.77 in. 

Pressure Ratio 3.61 

Stagnation Temperature 532 R 
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SPL, dB SPL, dB 


NOISE SPECTRA (CPN) 



15 19 23 27 31 35 39 43 




Freq. Band No. 


Freq. Band No. 


Band 15 is 125 Hz 


SUMMARY 


Source strength was evaluated using the PARC code with a k-e turbulence model 

The characteristic Strouhal no. was obtained from k and e ( with an empirical 
proportionality constant) 

Time-averaged velocity and temperature predictions were used to assess the 
sound/flow interaction 

Constants used in the supersonic convection factor are determined empirically 

The SPL directivity and spectra demonstrate favorable agreement with data ( 
specially at angles close to the jet axis) 

The empirical constants used in these predictions need to be investigated for other 
geometries 
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OUTLINE 


NUMERICAL SIMULATION OF UNSTEADY SUPERSONIC FLOW USING 
AN IMPLICIT ALGORITHM FOR THE STRONGLY COUPLED 


NAVIER— STOKES AND K-E EQUATIONS 
by 


S.H. Shiht A. Hamed^and J.J. Yeuan# 

* ICOMP, NASA Lewis Research Center 

# Dept, of Aero. Eng. & Eng. Mechanics 
University of Cincinnati, Cincinnati, 
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Fig. 41 Computational Grid for Supersonic Flow Over an Open Cavity. 
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Fig. 40. Schematic of the Solution Domain for Supersonic Flow Over an Open Cavity. 
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Instantaneous Density Contours 


a 
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Instantaneous Velocity Vectors 






Instantaneous Static Pressure Contours 


251 


o 


t 

„s _ 




a S 

S> u 

• V) KJ 
CU<U N 

X 

wclc: 


. \ 


1 1 pi rpuT| im |rnT-^ 


hJ 

\ 

X 



252 



■s- 

u « 

I- ’S 

So |*;g 

b, &<«<;■<! 
_+J 

•4-T «J '*“J +J «J 

a S ' 0 * as 

1) <u . 0) <u 
£ N Q.?i N 

£.Sm£.B 

CL.QS WOuPS 

I t ! S f 


*B I 


1 b _ 
I o 
<1 □ 


b 

I 

1 

I 

\3b 

I 

I 

I 

H 


P 

I 

II 
4 i 


\ 

>- 


r 1- 

o 


• 9 

T 


| 

s 

I 


1 

5 


(< u >) 'ias 


253 


CONCLUSIONS 
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IMPLEMENTATION OF A k-E MODEL 
IN SPECTRAL ELEMENT CODE 


Wai-Ming To 


Sverdrup Technology, Inc. 


OBJECTIVES 


• Study the flow physcis in problems that are generic to 
turbomachinery flows. 

• Assess accuracy of turbulence models. 
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EQUATIONS OF MOTION 


dui 

~dt 


duj 

dxi 


Mi- 


dp. 

dx; 



where v e - v / + v f 


TIME DISCRETIZATION 


• 2nd order predictor-corrector 

• Capacitance matrix method to solve for pressure 
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TIME DISCRETIZATION (cont.) 


.(l) - 


n . t a dp^ . d 
u n i +Atf i -At-^—+At-~ 

uX i OX j 




f du n i du n j 
+ — 


dx; dxi 


a4 1} 

dxj 


= 0 


u- 2) = u?+ At fi~ 


At dp^ + At d 


2 dxi 2 dx; 


+ 


duP 

dxi 


A t_ _d_ 
2 dx; 


= 0 


V e ^ 




r du n i 
+ — 


3Xy 


dxj 


3 U !'> 3^ 

5 + -/ 


dx,- 


ax,- 


TIME DISCRETIZATION (cont.) 
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dx. 


+ 


-n n + 1 

UU j 
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2 dxj 


v e"> 


aw? a*" 
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ax, ox,- 


where p = + p n 


257 


TURBULENCE MODEL 


• Eddy viscosity concept 

subgrid model — > 2-equation model 


Zero equation model 


// = k 2 y 2 


1 -e~ y+l26 


1 3 ? I 


v, = 0.018 Uoo 5* 


TURBULENCE MODEL (cont.) 


» Wall boundary conditions 
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(xA 2 /‘n)fKH 


NUMERICAL RESULTS (laminar flow) 
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SUBSONIC INLET FLOWS WITH TRANSITION 


BY 


Danny Hwang and Kyung Ahn 


September 16, 1993 


INTRODUCTION 

Real flow simulation is needed 

(laminar region + turbulent region) 


for: 

(1) cruise drag prediction 

(2) separation angle of attack prediction 
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Difficulties: 


• Transition is a function of many variables 

- free stream turbulence 
-surface roughness 
-streamwise surface curvature 

- pressure gradient 
-etc. 

• Skin friction coefficient which is the best 

transition indicator is difficult to calculate 


METHODS OF TRANSITION PREDICTION 

• Euler/boundary layer approach 

• Euler/boundary layer/stability analysis 

0 Navier-Stokes calculation with turbulence models 

- Baldwin-Lomax with transition 

- RNG turbulence model 
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TRANSITION PREDICTION 


• Flat plate 

• NACA0012 

• ADP inlet 
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pnickploler2.q. RNG model 



Boundary Layer Thickness of NACA001 2 Airfoil 
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Local Skin Friction Coetticient lor NACA0012 Airfoil 



ADP CONVENTIONAL NACELLE AT CRUISE MACH NUMBER 
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LA M 1 N A R/T l .1 R H U IJS N T TRA N SI TI 0 N 


bkui tilClloo to«llK*«**I.Uy»l-0 /a 



V "X_ 

I 1 « 4 • ' 



XJC 


Laininar/Turbulent Transition 


Laminar region Turbulent region 



Laminar region prediction: 

Lewis’ s boundary layer analysis : 8.17 % of chord length 
Langley’s stability analysis : 8.00 % of chord length 
On going projects: 

UNO turbulence model in I’AItt'dl) 
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CONCLUDING REMARKS 


Several methods have been tested to predict the 
transition of an ADP inlet. The prediction of the 
transition at cruise conditions agrees very well 
among the methods investigated except for the 
Baldwin-Lomax/Transition model which predicts 
early transition. The accuracy of prediction 
methods will be judged by comparing with 
upcoming laminar flow experiment. 
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Workshop on Computational Turbulence Modeling 
( Sept. 15-16, 1993 ) 


A Comparative Study of Turbulence Models in Predicting 
Hypersonic Inlet Flows 

by 

Kamlesh Kapoor 
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TEST CASE AND TURBULENCE MODELS CONSIDERED 

• THE NASA P8 INLET. WHICH REPRESENTS CRUISE CONDITION OF 
OF A TYPICAL HYPERSONIC AIR-BREATHING VEHICLE. WAS 
SELECTED AS A TEST CASE FOR PRESENT STUDY. 

• PARC2D CODE, WHICH SOLVES THE FULL TWO-DIMENSIONAL 
REYNOLDS-AVERAGED NAVIER-STOKES EQUATIONS. WAS USED 
FOR THIS STUDY. 

• THE RESULTS ARE PRESENTED FOR A TOTAL OF SIX VERSIONS OF 
ZERO- AND TWO-EQUATION TURBULENCE MODELS. 

• ZERO-EQUATION MODELS 

• THE BALDWIN-LOMAX MODEL 

• THE THOMAS MODEL 

• A COMBINATION OF THE B.L./THOMAS MODEL 

• TWO-EQUATION MODELS 

• LOW-REYNOLDS NUMBER MODELS 

• THE CHIEN MODEL 
•THE SPEZIALE MODEL 

• HIGH-REYNOLDS NUMBER MODEL 

* THE LAUNDER AND SPALDING MODEL 


EXPERIMENTAL BACKGROUND 

• THE EXPERIMENTAL INVESTIGATION OF THE P8 INLET WAS CONDUCTED 
AT NASA AMES' 3.5-FOOT HYPERSONIC WIND TUNNEL. 

. THE INLET WAS A MACH 7.4 RECTANGULAR MIXED COMPRESSION (WITH INTERNAL 
COMPRESSION RATIO OF 8) DESIGN WITH EXITING SUPERSONIC FLOW. 

• INLET COWL HEIGHT - 18.33 CM. 

•OVERALL LENGTH - 136.2 CM. 

• TEST CONDITIONS: 

•MACH NO -7.4 

•TOTAL PRESSURE - 4.14 x 10® N/rf 
•TOTAL TEMPERATURE - 811°K 

• REYNOLDS NO - 8.86 x 10 6 /m 

• MODEL WAS WATER COOLED AND ISOTHERMAL WALL CONDITIONS 
WERE MAINTAINED; THE WALLS TEMPERATURE - 302° K 

• THE TRANSITION POINTS: 

•CENTERBODY - 40 PERCENT FROM WEDGE L.E. EDGE TO INLET ENTRANCE. 
•COWL - HALFWAY BETWEEN INLET ENTRANCE AND THROAT. 
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X/XREF 



THE COMPUTATIONAL GRID 



• GRID SIZE WAS 221 x91. 

• GRID WAS NONUNIFORM IN X DIRECTION: 

• PACKED ON BOTH ENDS FROM THE WEDGE L.E. TO THE COWL L.E. 

• GEOMETRICALLY STRETCHED FROM THE COWL L.E. TO THE EXIT OF 
THE INLET. 

• IN Y DIRECTION. THE GRID WAS PACKED USING HYPERBOLIC TANGENT 
FUNCTION . YPLUS WAS APPROXIMATELY 1 AWAY FROM BOTH WALLS. 

• A SEPARATE GRID WAS MADE FOR THE LAUNDER AND SPALDING MODEL 
AND YPLUS OF APPROXIMATELY 30 WAS USED AWAY FROM THE WALLS. 
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BOUNDARY CONDITIONS 


FIXED 


SLIP 


CONTOUR LEVEL 



NON-REFLECTIVE NO-SUP ISOTHERMAL 


NO-SLIP ISOTHERMAL 


EXTRAPOLATION 


DENSITY CONTOURS FOR P8 INLET 
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PRESSURE CONTOURS FOR P8 INLET 


CONTOUR LEVEL 



MACH NUMBER CONTOURS FOR P8 INLET 

CONTOUR LEVEL 
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PITOT PRESSURE AND TOTAL TEMPERATURE 

DISTRIBUTIONS AT X/XREF = 6.09 PITOT PRESSURE DISTRIBUTITIO 






J38X/A 
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CONCLUSIONS 

• A COMPUTATIONAL STUDY HAS BEEN CONDUCTED TO EVALUATE 
THE PERFORMANCE OF VARIOUS TURBULENCE MODELS. 

• THE THOMAS MODEL COMPARES VERY WELL WITH THE EXPERIMENTAL 
DATA. AND IT PERFORMS BEST AMONG THE ZERO-EQUATION MODELS. 

• THE BALDWIN-LOMAX MODEL AND ITS COMBINATION WITH THOMAS 
MODEL ARE NOT ABLE TO RESOLVE THE PROBLEM OF SHOCK 
WAVE AND BOUNDARY-LAYER INTERACTION ACCURATELY. THE 
BALDWIN-LOMAX MODEL PREDICTS SEPARATION NEAR THE 
INTERACTION OF THE COWL SHOCK WITH THE WEDGE BOUNDARY 
LAYER, WHERE NONE IS KNOWN TO EXIST IN EXPERIMENTS. 

• THE CHIEN AND SPEZIALE MODEL COMPARE VERY WELL WITH THE 
EXPERIMENTAL DATA, AND PERFORMS BETTER THAN THE THOMAS 
MODEL, PARTICULARLY NEAR THE WALLS. THE LAUNDER AND 
SPALDING MODEL DOES NOT PERFORM AS GOOD AS THE CHIEN AND 
SPEZIALE MODELS. 

• AS THE CPU TIME REQUIRED FOR THE THOMAS MODEL IS FAR LESS 
THAN THE TWO-EQUATION MODELS. IT IS CONCLUDED THAT THE THOMAS 
MODEL IS BEST SUITED FOR THE PREDICTIONS OF PRESSURE 
DISTRIBUTIONS. AND THE CHIEN AND SPEZIALE MODELS ARE BEST TO 
CALCULATE FLOW QUANTITIES NEAR THE WALLS. 
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' 1 TOT PRESSURE AND TOTAL TEMPERATURE 
DISTRIBUTIONS AT THE THROAT 
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AN ALGEBRAIC TURBULENCE MODEL 
FOR TURBOMACHINERY 


by 


Rodrick V. Chima 
NASA Lewis Research Center 
Cleveland, Ohio 


279 


OVERVIEW 


• MOTIVATION - TURBINE ENDWALL HEAT TRANSFER 
. DESCRIPTION OF NEW MODEL 
. RESULTS 

1. FLAT PLATE 

2. ANNULAR TURBINE CASCADE 

3. TURBINE ENDWALL HEAT TRANSFER 

4. SUPERSONIC COMPRESSOR BLADE 
. SUMMARY 



EXPERIMENTAL ENDWALL STANTON NUMBER CONTOURS 
AS A FUNCTION OF 8 inlet AND Re chord 
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RVC3D (ROTOR VISCOUS CODE -3-DI 

BY R. V. CHIMA 
DESCRIPTION 

. EULER OR NAVIER-STOKES ANALYSIS 
FOR STEADY 3-D FLOWS IN TURBOMACHINERY 

FEAT UR ES 

. CARTESIAN FORMULATION, ROTATION ABOUT X-AXIS 
RECTANGULAR OR ANNULAR GEOMETRIES 

. SOLVES NAVIER-STOKES EQUATIONS 
THIN-LAYER FORMULATION, (NO STREAMWISE VISCOUS TERMS) 
RETAINS HUB-TO-TIP U BLADE-TO-BLADE VISCOUS TERMS 
BALDWIN-LOMAX OR CEBECI-SMITH TURBULENCE MODEL 
SIMPLE TIP CLEARANCE MODEL 

. NODE-CENTERED FINITE-DIFFERENCE FORMULATION 
EXPLICIT 4-STAGE RUNGE-KUTTA TIME-MARCHING SCHEME 
2ND + 4TH ORDER ARTIFICIAL VISCOSITY, EIGENVALUE SCALING 
VARIABLE Atjj & IMPLICIT RESIDUAL SMOOTHING 
HIGHLY VECTORIZED & AUTOTASKED FOR CRAY Y-MP 

. STACKED C-TYPE GRIDS 
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TURBULENT VISCOSITY PROFILE 



INTERMITTENT 

TURBULENCE 

OUTER LAYER 
INNER LAYER 


CEBECI-SMITH & BALDWIN-LQM AX MODELS 

INNER LAYER: PRANDTL-VAN DRIEST FORMULATION 
CEBECI-SMITH 
fii = pl*\dv/dy\ 

I = nyD 

D = 1 - exp(-y +/A*) VAN DRIEST DAMPING 
BALDWIN-LOMAX 

Hi = p/ J M 


OUTER LAYER: CLAUSER FORMULATION 
CEBECI-SMITH 
p, = KfryS'u, 

7 = J^l + 5.5 (f )°j KLEBANOFF INTERMITTENCY FUNCTION 

BALDWIN-LOMAX 

r/ /i • I y max J mar 

H 0 = K fn C cr mm| wake opt . on 

/(y) = yM# 
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BALDWIN -LOMAX MODEL ANALYSIS 

(SEE PAPER FOR DETAILS) 

1. ASSUME SUBLAYER-WALL-WAKE VELOCITY PROFILE 

2. CALCULATE BALDWIN-LOMAX FUNCTION f(y) 

MAX. OCCURS AT y mo , = ,646<5 
INDEPENDENT OF PRESURE GRADIENT 

NO MAX. FOR INFINITELY FAVORABLE dp/dx 

3. SPURIOUS MAX. CAN OCCUR AT EDGE OF VISCOUS SUBLAYER 
MOST LIKELY AT LOW Re & FAVORABLE dp/dx 




1 . VELOCITY PROFILES, Re, = 105,000 j. B-L FUNCTION /(y) 


3. SPURIOUS MAXIMUM IN /(») 


PJBflEQSED TURBULE NC E _M OREL 

INNER LAYER (SIMILAR TO BALDWIN-LOMAX) 

Pi = pl’M 
I = icy D 

D = 1 - exp(— y + /A*) 
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PROP O SED TURBULENCE MODEL 

PRESSURE GRADIENT EFFECTS 


. ACCELERATING FLOWS TEND TO REL AMIN ARIZE 
. MODELLED BY INCREASING A + IN FAVORABLE dp/ds 
. CEBECI’S EXPRESSION FOR A* USED: 


A* 


26 

v/1 + 11 -V 


v dp 
• 3 ds 


• PRESSURE GRADIENT EVALUATED USING: 


dp a K 


■Vp 


• “EDGE VELOCITY” EVALUATED AT A GRID LINE FAR ENOUGH 
FROM THE WALL TO GIVE THE GENERAL FLOW DIRECTION 

. KAYS-MOFFATT EXPRESSION WAS TESTED, EFFECTS TOO STRONG 


P_R.QE.QSm TU RBU L E NC E MOD E L 

LOC A L SHE A R MODEL 

. IN STRONGLY ACCELERATING FLOWS r + DECREASES WITH y + 
. MODELLED BY REPLACING t„« WITH r(y) IN D 


D = 1 - exp(-y + /A + ) 


V = V\ 


P (Mi + Me) ■ 
Pi Pi 


. ERROR IN ORIGINAL PAPER - USED p t \u\ ONLY 

. USED BY KAYS, PATANKAR-SPALDING, OTHERS 

. ALSO USED TO AVOID PROBLEMS AT SEPARATION WHEN t^i, — 0 
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BEQR QSED TURBULENCE modf.t. 

OUTER LAYER 


7 

Cut 

Ckm 


0.825 ,, 

0.55 


PROPOSED TURBULENCE MODEL 

OUTER LAYER - FUNCTION F 
. DEFINE F = ! fdy 

• INTEGRATE BY PARTS ASSUMING |w| — » 0 AS y -> 6 


F = /g^vlwldy 



= U V to - jJ ud V 

= / 0 ( u < - u ) rf y 
F = 6'u, 

. USE F DIRECTLY IN CEBECI-SMITH OUTER FORMULATION 
. ELIMINATES CONSTANT C ^ 

• DOES NOT REQUIRE KNOWLEDGE OF S OR u, 

• DISCOVERED INDEPENDENTLY BY D. A. JOHNSON, AIAA 92-0020 
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PROPOSED TURRTTT.WNCFI MODEL 



OUTER LAYER - LENGTH SCALE g 
. y IS THE CENTROID OF THE /(y) CURVE 

J% S(v)dy = j‘ fiv)dy 

• EVALUATE USING COLE’S VELOCITY PROFILES 

n y IS 

0 .5 

.5 .55 
oo .606 

• USE EQUILIBRIUM VALUE C*, lt = y/S = .55 


PROPOSED TURBULENCE MODEL 

OUTER LAYER - WAKE MODEL 

. LOWER OPTION IS A CONVENTIONAL WAKE MODEL 
. EVALUATE BY EQUATING TWO OPTIONS, ASSUMING 

y„ p = .606 S 
F,„ = u.i/2 
AV/u, as 1 

. GIVES = 0.825 
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P&QBQ SED TURBULENCE MODEL 

3-D IMPLEMENTATION 

• GRANVILLE BLENDING FUNCTION 

Mi 

Me// = Mo tanh — 

Mo 

• MODEL APPLIED INDEPENDENTLY IN BLADE-TO-BLADE (ij) 
AND SPANWISE (() DIRECTIONS 

. INNER LAYER - USE BULEEV LENGTH SCALE 
2 

y ' ~~ a, + *< +7*3 + 4 

. OUTER LAYER - USE ACTUAL DISTANCE ACROSS PROFILE 
y 0 = s, OR s< 

. BLEND 17 AND ( PROFILES VECTORALLY 

Miuri = 


FUT PLATE VELOCITY PROFILES 



COMPARISON OF FLAT PLATE VELOCITY PROFILES 
TO SPALDING’S COMPOSITE LAW OF THE WALL 
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STATIC PRESSURE DISTRIBUTION 



COMPUTED & MEASURED PRESSURE DISTRIBUTIONS 
FOR THE ANNULAR TURBINE CASCADE 
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COMPUTED & MEASURED LOSS COEFFICIENT PROFILES 
FOR THE ANNULAR TURBINE CASCADE 



COMPUTED & MEASURED EFFICIENCY CONTOURS 
IN THE WAKE OF THE ANNULAR TURBINE CASCADE 
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BOYLE’S LINEAR CASCADE 
STANTON NUMBER x 1000 
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BOYLE'S LINEAR CASCADE 
STANTON NUMBER x 1000 
R^. = 490,000 




flaBh— jam in miiiiMiini 


CONTOUR LCUELS 

n . nsruio 

U.U/bUU 
ii. f.Ttun 

«. r/r.ua 
ii. -im 
u.;*/*.i!u 

ii. 

h 
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TUROIIIENT UISCOSITV 
B-L ftOUEL 


L'UNIUUII HULLS 
: iHHitiii 
* liilmtii 
:» iHiniM 
l.ilSIiliHl 




CONTOUR LCUCLS 
i . nonon 
x.ouaua 

iKinriLi 

4.00000 

5. UOOOO 

6. U0W1I) 
/.UOUIllI 

« . uouoti 
h. tiniMin 
in. HtHtiMi 
111. IIINHMI 
XII. IIIMUMI 
M.IUUIIW 

4ti.nonnn 

•j U. UUIIOI) 
iiii.niNniii 

•iii.mwiin 

so. nonon 
mo.nnoo 

100.0000 

200.0000 

300.0000 

400.0000 

500.0000 
000.0000 

700.0000 



SUMMARY 

. SPURIOUS MAXIMUM IN B-L FUNCTION f{y) CAN GIVE INCOR- 
RECT TURBULENT LENGTH SCALE U ERRATIC St OR C , PATTERNS 

- MOST LIKELY AT LOW Re AND FAVORABLE dp/ds 
. NEW TURBULENCE MODEL PROPOSED 

- INTEGRAL RELATIONS FOR S’ u, AND S USED WITH C-S MODEL 

- EFFECTS OF dp/ds MODELED 

- WAKE MODEL PROPOSED 
. FLAT PLATE 

- B-L & NEW MODEL AGREE WITH LAW OF THE WALL 

- LOCAL SHEAR MOD. DOES NOT AGREE WITH LAW OF THE WALL 
. ANNULAR TURBINE 

- GOOD AGREEMENT WITH EXPT. PRESSURE DISTRIBUTION 

- WAKE MIXING UNDER-PREDICTED 
. TURBINE ENDWALL HEAT TRANSFER 

- VARIATIONS IN ENDWALL St WITH Re PREDICTED WELL 

- EFFECTS OF dp/ds IMPORTANT 
. TRANSONIC FAN 

- SHEAR LAYER FROM BOW SHOCK ACTS LIKE VISCOUS LAYER 

- NEW MODEL OVERPREDICTS L.E. p, 

- B-L MODEL PREDICTS REASONABLE L.E. p, 
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Lewis Research Center 


PROPULSION SYSTEMS DIVISION 


National Aeronautics and 

Space Administration COMBUSTION TECHNOLOGY BRANCH 



LOW EMISSIONS COMBUSTORS 


J. M. DEUR 

SVERDRUP TECHNOLOGY, INC. 
LEWIS RESEARCH CENTER GROUP 
BROOK PARK, OHIO 


Lewis Research Center PROPULSION SYSTEMS DIVISION 

National Aeronautics and 

Space Administration COMBUSTION TECHNOLOGY BRANCH 


APPLIED ANALYTICAL COMBUSTION/EMISSIONS RESEARCH 

• ANALYZE LOW EMISSIONS COMBUSTORS TO AID IN-HOUSE EXPERIMENTS 
AND CONTRACTOR COMBUSTOR DEVELOPMENT PROGRAMS. 

• PRESENT: 

• UTILIZE EXISTING CODES, PRINCIPALLY KIVA-II. 

• IMPROVE NUMERICS AND PHYSICAL MODELS (E.G., PDF COMBUSTION- 
TURBULENCE INTERACTION) ON LIMITED BASIS TO SATISFY 
CRITICAL NEEDS. 

• FUTURE: 

• ADOPT 3-D ALLSPD CODE WHEN AVAILABLE. 
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Lewis Research Center PROPULSION SYSTEMS DIVISION 

National Aeronautics and 

Space Administration COMBUSTION TECHNOLOGY BRANCH 


KIVA-II FEATURES 

• MULTI-DIMENSIONAL TIME ACCURATE FINITE DIFFERENCE CODE. 

• COMPRESSIBLE FLOWS. 

• k-i TURBULENCE MODEL WITH WALL FUNCTIONS OR SUB-GRID SCALE 
TURBULENCE MODEL. 

• LAMINAR KINETICS FOR ARBITRARY REACTION SET WITH QUASI- 
EQUILIBRIUM OPTION (MIXING CONTROLLED COMBUSTION MODEL ALSO). 

• STOCHASTIC SPRAY MODEL WITH VAPORIZATION, AERODYNAMIC 
BREAKUP, TURBULENT DISPERSION, AND COLLISION SUB-MODELS. 

• ADIABATIC OR CONSTANT TEMPERATURE WALL BOUNDARIES. 

• ARBITRARY MESH. 



Lewis Research Center PROPULSION SYSTEMS DIVISION 

National Aeronautics and 

Space Administration COMBUSTION TECHNOLOGY BRANCH 


APPLIED ANALYTICAL COMBUSTION/EMISSIONS RESEARCH TEAM 



Carnegie Mellon University and University ot South 
Florida leaving program at end ol FY93 
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RICH BURN - QUICK MIX - LEAN BURN (RQL) FLAME TUBE 
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ROL FLAME TUBE WALL TEMPERATURE COMPARISON 



SLANTED SLOT MIXER ( W/O SWIRL) 

TWO HOLE MIXER (W/O SWIRL) 

ONE HOLE MIXER (W/O SWIRL) 

ONE HOLE MIXER (W/ SWIRL) 


MIN 

MAX 


Lewis Research Center 
National Aeronautics and 
Space Administration 


PROPULSION SYSTEMS DIVISION 
COMBUSTION TECHNOLOGY BRANCH 



ROL FLAME TUBE TEMPERATURE COMPARISON 



SLANTED SLOT MIXER (W/O SWIRL) TWO HOLE MIXER (W/O SWIRL) 



ONE HOLE MIXER (W/O SWIRL) ONE HOLE MIXER (W/ SWIRL ) 
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Lewis Research Center 
National Aeronautics and 
Spaoo Administration 


PROPULSION SYSTEMS DIVISION 
COMBUSTION TECHNOLOGY BRANCH 



RQL FLAME TUBE NO x EMISSION INDEX COMPARISON 



ONE HOLE MIXER (W/O SWIRL) ONE HOLE MIXER (W/ SWIRL) 


MIN NOTE : CALCULATION CONSIDERS THERMAL NO x ONLY. MAX 


I Lawla R.warch C.nt.r PROPULSION SYSTEMS DIVISION 

iSS fc 

National Aeronautics and 

Space Admiration COMBUSTION TECHNOLOGY BRANCH 


RQL MIXER CONFIGURATION PARAMETRIC STUDY NO x COMPARISON 

SLANTED SLOT MIXER W/O SWIRL 

1.00 

TWO HOLE MIXER W/O SWIRL 

1.25 

ONE HOLE MIXER W/O SWIRL 

1.39 

ONE HOLE MIXER W/ SWIRL 

0.75 

SLANTED SLOT MIXER W/O SWIRL (EXPERIMENTAL) 

1.00 

NOTE: 


• VALUES ARE SPATIAL AVERAGES TAKEN AT SAMPLING LOCATION B. 

♦ VALUES ARE NORMALIZED BY SLANTED SLOT MIXER EXPERIMENTAL 

READING. 
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THE EFFECTS OF TURBULENCE MODELING ON THE 
NUMERICAL SIMULATION OF CONFINED SWIRLING FLOWS 


G. J. MICKLOW AND M. R. HARPER 
UNIVERSITY OF FLORIDA 
GAINESVILLE, FLORIDA 

J. M. DEUR 

SVERDRUP TECHNOLOGY, INC. 
BROOK PARK, OHIO 


AIAA/SAE/ASME/ASEE 29TH JOINT PROPULSION CONFERENCE 
MONTEREY, CALIFORNIA 



lewis Research Center PROPULSION SYSTEMS DIVISION 
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K - £ TURBULENCE MODEL 


DK = d_(v L 3K' 


Dt 3x; I cty 3x; 


3ui | 

f3uj 

duj 

+ *■ 

3xj 1 

( 3xj 

3xj 


- e 


K = 


c 4 K 

O, = 


= 0.09 
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K - £ TURBULENCE MODEL 




C 1 =1.44 
C 2 = 1 .92 


De _ a /u t ae 


Dt 9xi 1 a F 0xi 




£ Suj 

+ C 1»t-T 1 
K 3X: 


1 3uj 3uj ' £ 2 

| dxj’ + 9xj] “ C 2 "k" 


L.wls Restarch C.nltr PROPULSION SYSTEMS DIVISION 

National Aeronautics and 

Space Administration COMBUSTION TECHNOLOGY BRANCH 



MODIFICATIONS TO K - £ TURBULENCE MODEL 


1° 


^ I wall 


A u„ 

= 7 

= 0.0083 
= 0.09 
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C 2 MODIFICATIONS TO K - E TURBULENCE MODEL 

c 2 

= C 2 '(1-C 3 Ri) (0.1 < C 2 < 2.4) 

c 2 * 

= 1.92 

C 3 

= 0.2 

Ri 




wall 

= 1.92 
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UCI AXISYMMETRIC CAN COMBUSTOR 
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C 2 VARIATION IN UCI AXISYMMETRIC CAN COMBUSTOR 
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TURBULENT TIME SCALE IN UCI AXISYMMETRIC CAN COMBUSTOR 




Lewis Research Center 
National Aeronautics and 
Space Administration 


PROPULSION SYSTEMS DIVISION 
COMBUSTION TECHNOLOGY BRANCH 



TURBULENT VISCOSITY IN UCI AXISYMMETRIC CAN COMBUSTOR 



C M CASE 



C 2 CASE 



0.0 crr?/s 1.34 an 2 /s 
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VELOCITY VECTORS IN UCI AXISYMMETRIC CAN COMBUSTOR 


II 


PPHIIIII 

jSlllllllllllilii 

mrnmW 



.ttilllMi •III- i 


BASELINE CASE 


C u CASE 


C 2 CASE 


Lewis Research C.nler PROPULSION SYSTEMS DIVISION 

National Aeronautics and 

Soace Administration COMBUSTION TECHNOLOGY BRANCH 



COMPARISONS TO UCI AXISYMMETRIC CAN COMBUSTOR DATA (z=24 cm) 


«SOO 
*000 
m J300 

N 
€ 

w JOOO 

* 

£■ 2500 

> 2000 
< 1500 

1000 
500 

-1012345 

Rodaol Loco Iron (*. cm) 




AXIAL VELOCITY 


CIRCUMFERENTIAL VELOCITY 
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Lewis Researcn Center 
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COMPARISONS TO UCI AXISYMMETRIC CAN COMBUSTOR DATA lz=24 cm) 



VELOCITY CORRELATION 


TEMPERATURE 


Lewis Research Center PROPULSION SYSTEMS DIVISION 

National Aeronautics and 

Space Administration COMBUSTION TECHNOLOGY BRANCH 

CONCLUSIONS 

• INLET VELOCITY IS BEING MODIFIED TO MATCH FIRST INLET STATION. 

• PDF COMBUSTION-TURBULENCE MODEL OF HSU, ET AL. IS BEING ADDED. 

• MUCH MORE WORK IS REQUIRED. 
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TITLE 


DEVELOPMENT OF A 
RELIABLE ALGEBRAIC TURBULENCE 
MODEL GIVING ENGINEERING ACCURACY 
AT REASONABLE COST 


by 


B.P. Leonard and J.E. Drummond 
The University of Akron 
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OUTLINE 


I. CHOICES FOR TURBULENCE MODELLING 

A. ZERO-EQUATION (ALGEBRAIC MODELS) 

B. MULTIPLE-EQUATION MODELS 

II. THE NEED FOR A RELIABLE ALGEBRAIC MODEL 

A. RELATIVE SIMPLICITY 

B. ENGINEERING ACCURACY 

C. COST-EFFECTIVE APPLICATION 

III. DEVELOPMENT OF THE MODIFIED MIXING LENGTH (MML) 
MODEL 

IV. COMPARISON OF MODELS 

A. MML 

B. BALDWIN-LOMAX 

C. TWO-EQUATION (k-e) MODEL 

V. WHERE DO WE GO FROM HERE? 


LeRC 
CTM 
WORKSHOP 


9/93 


THEORETICAL BASIS OF THE MODEL 


EFFECTIVE VISCOSITY: p tff = p + p, (1) 

TURBULENT VISCOSITY: p, = p l 2 |co| (2) 


MIXING LENGTH: 


1 - 1 


*\c. 


C, 

fC 


1 - «p(-^ 


f = k (c L ) 3 '‘ for >c > 


for y * < C, (3) 


(4) 


where: y' - yly'\ y' = pI\[p\t w \ 


SHEAR STRESS “FILTER”: 

|t ( | = 0.1 lr M | + 0.2 | r M { + 0.4 |rj + 0.2 \r M \ * 0.1 |r w | (5) 
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REFINING MML - CAPPING LENGTH SCALE AS 
A FUNCTION OF POSITION 


4 


OUTER REGION LENGTH SCALE (CAPPING LENGTH): 

* = 'cap = y*. 

^2 


CAPPING LENGTH BASED ON BOUNDARY LAYER THICKNESS: 

^CAP = 


CAPPING LENGTH BASED ON LOCAL SHEAR STRESS: 
*cap = (3.31-10' 6 )B|C f |' 3 - 5 y’ 


LeRC 

CTM 

WORKSHOP 


VARIATION OF MIXING LENGTH 
IN ATTACHED BOUNDARY LAYER 
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9/93 

EFFECT OF VARYING C, IN EQUATION (3) 
WHILE KEEPING C 2 FIXED (=5) 

6 
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WORKSHOP 

9/93 

EFFECT OF VARYING C^ IN EQUATION (3) 
WHILE KEEPING C,/C 2 FIXED (=8) 

7 
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SHEAR LENGTH SCALE, y*, 
BASED ON | r w | AND EQUATION (5) 

8 
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THEORETICAL BASIS OF 
BALDWIN-LOMAX MODEL 

9 

INNER REGION TURBULENT VISCOSITY: 

fj. t = p l 2 |w| 


WHERE 


£ = k y [1 - exp (~y*IA *)] 


OUTER REGION TURBULENT VISCOSITY ASSUMES THE SMALLER 
VALUE OF: 

P, = pK C cp F k (y) y mu F mu 

OR 


= 0.25 P KC cp Ul„y m JF m „ 


WHERE 


3 

II 

✓ — N 
£> 
ft. 




F t (y) = [l ♦ 5.5(EZjT 

y max ' 
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FLAT PLATE GEOMETRY AND GRID 


11 


UNIFORM 
MACH 0.2 
FLOW 


/ 


DEVELOPING BOUNDARY LAYER 
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MML FLAT PLATE SKIN FRICTION 
CALCULATIONS FOR VARYING C t /<5 + 

12 
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MML FLAT-PLATE VELOCITY DEFECT 
PROFILES FOR VARYING C,/5 + 

13 
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SKIN FRICTION CALCULATIONS FOR 
VARIOUS TURBULENCE MODELS 


14 
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VELOCITY DEFECT PROFILES FOR 
VARIOUS TURBULENCE MODELS 


15 



50 
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TURBULENT VISCOSITY PROFILES 
FOR VARIOUS MODELS 


16 
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CONVERGENCE CHARACTERISTICS FOR 
FLAT PLATE CALCULATIONS 

17 


MODEL 

ITERATIONS 

CRAY Y/MP 
CPU TIME (s) 

MMLT 

4000 

500 

Baldwin-Loraax 

4000 

500 

Chien k-E 

5000 

.1100 
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GEOMETRY FOR BACKWARD 
STEP CALCULATIONS 


yyyyy//////////fy/y///y ////////// ///////y/////y 




SPECIFIED 

INFLOW 

VELOCITY 

PROFILE 




H. 


H. 

• 137 cm 

Ur.t 

> 44.2 m/s 

M r«( 

>.128 


7 777 7777/77 7777/7 7 7 7/7777/7777777 


4H. 
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GRID FOR BACKWARD STEP CALCULATIONS 

19 

WORKSHOP 






(a) Entire grid 



(b) Detail near step 
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INCORPORATING THOMAS MODEL 
FOR FREE SHEAR LAYER 


20 


LENGTH SCALE FOR UNBOUNDED REGION FOR THOMAS MODEL: 
( = M Max 0 u jl> ~ Mindiijl)] 


TURBULENT VISCOSITY IN TRANSITION REGION BETWEEN 
BOUNDARY LAYER AND UNBOUNDED REGION: 


Mtr 


Mmml( C 4 - y*) + M Th (y* - C 3 ) 


c - c 


LeRC 

CTM 

WORKSHOP 


9/93 


WALL BOUNDARY AND FREE SHEAR LAYERS 


21 
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BACKWARD-FACING STEP SKIN FRICTION 
CALCULATIONS FOR VARIOUS MODELS 



22 


L 
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REATTACHMENT POSITION FOR 
VARIOUS TURBULENCE MODELS 

23 


CASE 

REATTACHMENT 

posmoN 

(STEP HEIGHTS.Hd) 

Driver-Seegmiller Daia 

6.25 

MMLT (C|/&+ = .7, 

C 3 = IOAcap.C* = 20^'cap) 

7.42 

Thomas 

12.28 

Baldwin-Lomax 

5.41 

Chien k-E 

6.27 


316 


LeRC 

CTM 

WORKSHOP 
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VELOCITY PROFILE COMPARISON FOR 
VARIOUS TURBULENCE MODELS 

25 
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TURBULENT VISCOSITY CONTOURS FOR 
BACKWARD-FACING STEP 


27 
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CONVERGENCE CHARACTERISTICS FOR 
VARIOUS TURBULENCE MODELS 


28 


MODEL 

ITERATIONS 

CRAY Y/MP 
CPU TIME (s) 

MMLT 

45,000 

7100 

Thomas 

40,000 

6000 

Baldwin-Lomax 

45,000 

7000 

Chien k-£ 

1 10,000 

50,000 
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WHERE DO WE GO FROM HERE? 

29 


A. APPLICATION TO MORE COMPLEX FLOWS 

• NOZZLE FLOWS 

• RECIRCULATING INTERNAL FLOWS 

• EXTERNAL FLOWS 

B. COMPARISON WITH OTHER TURBULENT CALCULATIONS 


C. ESTABLISH REASONABLE ALGEBRAIC MODEL AS AN OPTION 
TO GIVE RESULTS WITH ENGINEERING ACCURACY FOR A 
VARIETY OF APPLICATIONS 
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Application of Algebraic and 
Two-Equation Turbulence Models 
to HSR Nozzle Flow Calculations 


N. J. Georgiadis 
and 

J. R. DeBonis 

Nozzle Technology Branch 
Propulsion Systems Division 


September 16, 1993 


FULL NAVIER-STOKES ANALYSES 
OF NOZZLES (WITH PARC CODE) 


• OVERVIEW OF PARC CODE AND AVAILABLE TURBULENCE 
MODELS 


• VALIDATION TEST CASES: 

1. AXISYMMETRIC PLUG NOZZLE 

2. EJECTOR NOZZLE 

• HIGH-SPEED RESEARCH (HSR) NOZZLES WITH 
POTENTIAL FOR NOISE REDUCTION: 

1. NASA/GE 2DCD MIXER/EJECTOR NOZZLE 

2. PRATT & WHITNEY 2D MIXER/EJECTOR NOZZLE 


321 


OVERVIEW OF PARC : 

• 3D AND 2D/AXISYMMETRIC VERSIONS 

• NAVEER-STOKES AND EULER MODES 

• CENTRAL DIFFERENCE DISCRETIZATION 

• BEAM AND WARMING ALGORITHM 

• SEVERAL TURBULENCE MODEL OPTIONS 

• CAPABILITY TO HANDLE MULTIPLE GRID BLOCKS 

(NONCONTIGUOUS INTERFACING) 

• GENERALIZED BOUNDARY CONDITIONS 


Full Navier-Stokes Equations in PARC 


dQ 9Fj _ 1 dGj 

3t 9xj Re 0Xj 


pu, 


0 

pu,Uj + P5 U 



(E + P)Uj 


Uktjk - q, _ 
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t,j = n 


5ij. 

V 3x i 3x„ 


2 9u. _ 


= 


-K, 


9T 


(Y-l)Pr dx i 


H e « = H + H, 


K c-eff ~ K c + 


M- t Cp 

~TT 


Thomas Model 


m = p *>l 


Wall bounded flows: 


t = Ky 


1 - e 


A* 


Free shear layers: 


*. [ Max (hl) - Min (N)] 

®c 
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Baldwin-Lomax Model 


= 


(dinner • y ^ y o 

(M-t)omer » y — y<: 


(Pl)ouler ~ P^CcpF^F, 


kJeb 


F W a)ce — (y max^max > niax^dir ^ F*max ) 


F(y) = y|co| 


f ll) 

1 - e A ‘ 

F)deb(y) - 

1 + 5.5 f C|tIeb y ) 

\ J 


V ymax ) 


Chien k-£ Model 


P, = c„f M pk 2 /e 


D(pk) _ _a_ 
Dt dx. 



+ n-pe-D 


D(pe) 

Dt 


dx. 




+ C E |f,n| - C e2 f 2 p^ - 2p^-e ( -°- 5) '* ) 
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3uj 3^ 
3xj 3Xj 

D = 2n k / y 2 
f„ = 1 - e ( ' °" 5y ) 

f, = 1.0 

f 2 = 1.0 - 0.22 e ~^ R ' /6 ^ 



Re, 



TWO - EQUATION MODEL OF SPEZIALE 


• SOLVED AS A k-e MODEL IN PARC 

• ORIGINALLY DEVELOPED AS A k-x MODEL 
IN NASA CR-182068 

• UPWIND DIFFERENCE DISCRETIZATION 
(MAIN FLOW EQUATIONS IN PARC USE 
CENTRAL DIFFERENCING) 
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BENCHMARK VALIDATION CASES 


1. Mach 0.2 flat plate 

2. Mach 0.128 backward-facing step 


SHEAR STRESS PREDICTIONS WITH THE 
DIFFERENT TURBULENCE MODELS IN PARC 
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SHEAR STRESS AFTER THE STEP WITH THE 
DIFFERENT TURBULENCE MODELS IN PARC 



x/H b 


LANGLEY SINGLE FLOW PLUG NOZZLE 

• VENTED AND NON-VENTED PLUGS 
« 15° PLUG HALF ANGLE 

• HEAVILY INSTRUMENTED TO MEASURE: 

1. PLUG SURFACE TEMPERATURES, PRESSURES, SHEAR STRESS 

2. JET PLUME QUANTITIES ( INCLUDING LDV & FLOW VISUALIZATION) 

3. FLOWFIELD ACOUSTICS 



327 



GEOMETRY FOR GRID GENERATION 
AND PARC2D CALCULATIONS 



VELOCITY PROFILES FOR FINE GRID SOLUTIONS 


X = 25 in. 



y |ii/») 


X = 50 in. 



X = 75 In. 



0 500 1000 1500 2000 2500 3000 

V (II/.) 


500 1000 1500 2000 2500 3000 
V (It/.) 
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1 


TOTAL TEMPERATURE PROFILES 
FOR FINE GRID SOLUTIONS 


X = 25 in. X = 50 In. X = 75 in. 



SHOCK FUNCTION 
(BASED ON PRESSURE GRADIENT) 





Thomas 


Baldwin-Lomax 

/Thomas 


K-Epsilon 
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STATIC PRESSURE DISTRIBUTION 


ALONG PLUG 


THOMAS 


BALDWIN-LOMAX / THOMAS 




K - EPSILON 



2D EJECTOR NOZZLE TEST CASE 


• REPRESENTATIVE OF MIXER-EJECTOR NOZZLES 
THAT ARE BEING CONSIDERED FOR SUPERSONIC 
TRANSPORT APPLICATION 


• FLOW DOMINATED BY TURBULENT MIXING OF A 
HIGH ENERGY STREAM WITH SECONDARY AIR 
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* 

2D EJECTOR NOZZLE OF 
GILBERT AND HILL 



VELOCITIES IN MIXING SECTION 


X r 7 in. 



X = 10.5 in. 
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TOTAL TEMPERATURES 
IN MIXING SECTION 


X = 3 In. 


X = 10.5 In. 




MIXING EFFECTIVENESS 
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TURBULENT VISCOSITY CONTOURS 
FOR 2D EJECTOR NOZZLE 


/ M- 


(I.IMUUHI 
IIJH.IJUlill 
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mu. iiiuifi 

Miii.imim 
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.'llll. IIMIIII 

min. mil hi 
■ mit.imun 
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i iiiiMum 

mi min 

>vi< 
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AXISYMMETRIC INVERTED VELOCITY 
PROFILE NOZZLE OF VON GLAHN, 
GOODYKOONTZ, AND WASSERBAUER 
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AXISYMMETRIC UNDEREXPANDED NOZZLE 
OF HELSTLEY AND CROSSWY 




Mixer/Ejector Nozzles 

• Entrain large amounts of secondary flow 

• Rapidly mix two flows together to lower 
jet velocity 

• Lower jet velocity results in lower noise 

• Maintain high thrust due to large mass 
augmentation 

F = mv 
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Typical Mixer/Ejector Nozzle 
Chute Geometry 



Primary Flow 
► Secondary Flow 


NASA/GE 2DCD 
Mixer/Ejector Nozzle 



335 



Shroud Static Pressures 


MAR = 1.20 

Thomas Turbulence Model 



X/L 
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amb 




Shroud Static Pressures 
MAR = 1.20 



X/L 
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Velocity Vectors in Mixing Section 


A e /A mix = 


■< i in u 4 ,i&.n 

1 1 1 |1 >C 1 1 '/>*» 1 1 1 1 1 s i«u» , , ,,/pi 1 1 1 
miiir'' miiiiui 1 1 mm' M u mu\ 

i mi;,''" ’ ' ' "«Wi Mini ii.r"" ’ ’ 1 '""mu 1 1 


i mm' 

i i l 1 1 „i v "-' » • 


•mu ii tit mw ' 




X/L = 0.0 
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Anatomy of the HSR P&W 2-D Mixer Ejector Nozzle 



PARC3D Thomas 


Flow Angle Contours 
Nozzle Exit Plane 


MW 12 


Very Large 
Flow Angles 


il : 


PARC3DK-E 


Inventory: 

• 2 Mixer Designs 

• 3 Shroud Lengths 

• 3 Shroud/Sidewall 
Acoustic 
Treatments 

• Sidewalls with 
Windows for Flow 
Visualization 
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3D COMPUTATIONAL 


GRID FOR NOZZLE FLOWFIELD 



P&W MIXER/EJECTOR NOZZLE 
Mach Number Contours 


LEVEL 
o . onooo 
o. iooon 
ii. viiiiun 
ii. million 
II. IIIIIOU 
II. Mill) 10 

n. minim 
It. VIIIHMI 
II. IJIIIII1H 
II . *11111(111 

I .1100110 

i inwin 


! -.|| 0 |||» 

1 . liOOIM) 

i . 7nooo 
i .onooo 
1 . 900U0 

2 . U00U0 
2. 10000 
2.20000 
? . noooo 



Valley Side 
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2nd GENERATION MIXER EJECTOR ANALYSIS 


Exit Plane Total Temperature 



92 - 3 - 28-6 


Conclusions 


• Algebraic models are inadequate for 
predicting mixer nozzle flowfields 

• Two-equation model available in PARC3D 
give no improvement over algebraic models 

• A reliable two-equation model is required 
in PARC3D for complex nozzle flows 
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Turbulence Modeling in Aircraft Icing 


Mark G. Potapczuk 

Icing & Cryogenic Technology Branch 


Workshop on Computational Turbulence Modeling 
NASA ICOMP 
September 15, 1993 


INTRODUCTION 


The Icing and Cryogenic Technology Branch develops computational tools 
which predict ice growth on aircraft surfaces and uses existing CFD 
technology to evaluate the aerodynamic changes associated with such 
accretions. 


Surface roughness, transition location, and laminar, transition, or turbulent 
convective heat transfer all influence the ice growth process on aircraft 
surfaces. 


Turbulence modeling is a critical element within the computational tools 
used both for ice shape prediction and for performance degradation 
evaluation. 
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CURRENT CODE DEVELOPMENT 


2D CODES 



• LEWICE/IBL - POTENTIAL FLOW / INTERACTIVE BOUNDARY LAYER 



3D CODES 

• LEWICE3D - PANEL CODE / INTEGRAL BOUNDARY LAYER 

• LEWICE3DGR - ANY GRID BASED FLOW SOLUTION 


ICE ACCRETION MODELING 

CURRENT MODEL USED FOR ICE GROWTH 


• MASS AND ENERGY BALANCE IN CONTROL VOLUMES ALONG THE 
SURFACE 


• CONVECTIVE HEAT TRANSFER IS MAJOR FACTOR IN ENERGY 
BALANCE 


• INTEGRAL BOUNDARY LAYER FORMULATION USED TO DETERMINE 
LAMINAR AND TURBULENT HEAT TRANSFER COEFFICIENTS 


• SURFACE ROUGHNESS MODELED AS SAND-GRAIN ROUGHNESS; 
ACTUAL ICE ROUGHNESS VARIES FROM SMALLER TO LARGER THAN 
BOUNDARY LAYER THICKNESS 
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ICE ACCRETION MODELING 

CONVECTIVE HEAT TRANSFER MODEL USED FOR ICE GROWTH 
SKIN FRICTION COEFFICIENT 


C J 

2 


/864.O0 


0.1681 


In 


+ 2.568 



WHERE 



0M56, S V 3.86 
y4.11 Js tr e 
e 


-| 0.8 

+ W 


ICE ACCRETION MODELING 

CONVECTIVE HEAT TRANSFER MODEL USED FOR ICE GROWTH 
LAMINAR 


\(s) = 0.296^ [Vf 88 ^' 88 *] 


- 1/2 


TURBULENT 


h t (s) = StpVc = 


c/2 


Pr, + jF/2(l/Sl k ) 


p Vc 

r e p 
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ICE ACCRETION MODELING 

CONVECTIVE HEAT TRANSFER MODEL USED FOR ICE GROWTH 
ROUGHNESS STANTON NUMBER 

Vk ~ 0 2 

S' k = 1-16 (-^) 

AND 

^ = VeF/2 


ICE ACCRETION MODELING 

ICE ROUGHNESS CHARACTERIZATION 
SAND-GRAIN ROUGHNESS 
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ICE ACCRETION MODELING 


PLANS 


• EXPERIMENTS TO CHARACTERIZE ICE ROUGHNESS GEOMETRIES AT A 
VARIETY OF ICING CONDITIONS 


• EXPERIMENTS TO CHARACTERIZE VELOCITY FIELD OVER REAL AND 
ARTIFICIAL ICE ROUGHNESS GEOMETRIES 

• EXPERIMENTS TO MEASURE HEAT TRANSFER OVER REAL AND 
ARTIFICIAL ICE ROUGHNESS GEOMETRIES 


• DEVELOPMENT OF MODIFIED COMPUTATIONAL MODEL BASED ON 
THESE EXPERIMENTS 


ICED AIRFOIL AERODYNAMICS 


NACA 0012 ICING CONDITIONS 
a = 4° V = 130 mph 


d = 20 jim 
T = 18 F 


LWC = 2.1 g/rrr 


SHEAR LAYER 


RECIRCULATION 

REGION 


5 MINUTE 
ICE GROWTH 
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ICED AIRFOIL AERODYNAMICS 

BALDWIN-LOMAX TURBULENCE MODEL 


Inner Laver 

m v~ < 2 I ( v v * ) I 

/ = mixing length 

Outer Laver 

Vr F max Y max 

F (y) = ^|coj ( 1 - exp ( (-y) / A) ) 


ICED AIRFOIL AERODYNAMICS 

BALDWIN-LOMAX TURBULENCE MODEL 


NORMAL B.L. F PROFILE 



RECIRCULATION REGION 
F PROFILE 
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ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 


\i t = pJ 2 (o 

WHERE 

1-e > 


y + Ay 


<C, l(y) = k 




1 - 




1 - 


AND WHERE 


y + Ay 

* 

y 


>q 


Hy) 



ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 



* 


y 


V 

— * 
u 
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ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 


THE CEBECI-CHANG ROUGHNESS MODEL IS ADDED TO THE 
TURBULENCE MODEL 


A/ = 


0.9 


-k + \ 

s 


JZ-K ex p(-s 


y 


5 <k + <70 


V 0.7(0 


0.58 


70 <k + < 2000 

S 


WHERE, 


a/ = (zV) (M x /v) and k* = k $ (u/v) 


ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 
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ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 





ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 
STRUCTURED GRID FOR ARTIFICIAL ICE SHAPE 
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ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 


UNSTRUCTURED GRID FOR ARTIFICIAL ICE SHAPE 



*0.10 -0.0S 0.00 0.0S 0.10 0.1S 0.20 


ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 
STRUCTURED GRID MACH NUMBER CONTOURS 
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ICED AIRFOIL AERODYNAMICS 

MML TURBULENCE MODEL 


UNSTRUCTURED GRID MACH NUMBER CONTOURS 



CONCLUDING REMARKS 


• TURBULENCE MODELING PLAYS A ROLE IN ICE GROWTH PREDICTION 
AND IN PERFORMANCE EVALUATION 


• NEW MODELING IS REQUIRED FOR THE LARGE ROUGHNESS 
ELEMENTS OF A TYPICAL ICE ACCRETION 


• AN EXPERIMENTAL PROGRAM IS CURRENTLY UNDERWAY TO 
DEVELOP A DATABASE FOR CREATION OF SUCH A MODEL 


. AN ALTERNATE ALGEBRAIC TURBULENCE MODEL HAS BEEN USED TO 
EVALUATE PERFORMANCE DEGRADATION DUE TO ICING 


. THE MML MODEL HAS BEEN USED IN AN UNSTRUCTURED GRID 
NAVIER-STOKES CODE TO CALCULATE FLOW OVER AN ARTIFICIAL ICE 
SHAPE 
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Applied RNG Algebraic Turbulence Model for 
Three-Dimensional Turbomachinery Flows 

K. R. Kirtley 


Sverdrup Technology, Inc. LeRC Group 
and 

Cambridge Hydrodynamics, Inc. 


What is RNG ? 

Start — ►Navier-Stokes 


RNG 

1) Decompose velocity into low and 
high wave number components 

2) Use perturbation theory to eliminate 

high wave number bands then renormalize 
spectrum and repeat to Infrared cutoff 

3) Correlations disappear through mode 
elimination procedure 

4) Closure Is automatic 

5) Evaluate coefficients from high Re 
limit (fixed point) of perturbation 
expansion 


Reynold’s Averaging 

1) Decompose velocity into mean 
and fluctuating components 

2) Ensemble (time) average over 
entire spectrum 

3) Correlations arise 

4) Model correlations (Bousslnesq) 

5) Evaluate coefficients from generic 
flow data, e.g., wake, Jet, etc. 


Self Consistent Problem Dependent 

Result — ►effective viscosity 


RNG-Based Algebraic Model 

1/ = I/O 


where 


2 duk\ du{ 

3 13 dxk ) dxj 

A f - Infrared cutoff = top of the Inertial range 
resolvable < A f< modelled 

a =.1186 from -5/3 spectral law 

Acc 

(J c = yr from a and the ultraviolet dissipation range cutoff 


( dv,i dUi 

F J '+ a” 1 
oxj axi 


a 


1 + F(^6A 7 4 -C c 


€ = P = V S 
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How do you go from Fourier space to physical space ? 


why 2tx of course. 

r i 2?r 

Integral scale, L f = — 

A / 

But we need a mixing length, therefore 
Assume E = Ck€*A~* 


C k = 1.6075 


3 a 

Integration from a t to infinity gives k = -CrA.^ 3 £ 3 


Assume In high Re limit r2o 1 

i/ = 2 = G 


fc 2 


.9 


i 

a< 


Then L f = ( 2 C ll Cl)‘Aj 1 =^I> 


Giving 


i' = i/q 


1+H (^ L ^ _Cc 


Attributes of the RNG-Based Model 


1) The Heaviside function mimics: 
a) Near wall damping 

V = Vo 



9.2 


b) Intermlttency 

v = v 0 
e -* 

c) Transition 


when 


e < 


C c v 


— £4 

0 In the outer flow ^ 


3 

0 


See Next Slide 


2) Energy-Based ►-Non-equIllbrlum effects can be Included 
through clever manipulation of the dissipation rate 
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=r t 0 
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t n> ~ 
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3il 

1 <3 


% i 
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Cubic vs. Quartic 



Good: Analytic Solution 

Bad: Multivalued 


v^-yy/- Vo 1 ) - o 

Single Valued 
Highly Non-linear 



40 80 
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Couette Flow 


Near Wall Behavior 


Friction Factor 
vs. 

Reynolds Number 




Turbulence Length Scales 


In boundary layers 


C-4- 


l i = °»* tanh ( 5 ) 


K = J^»T2 and (7 = 0.0845 from RNG 


In wake region use Raj & Lakshmlnarayana correlation for 


cascade wakes 


h 


min (/cs, C w b ) 


b = 6 0 + cCl 135 

Cd = coefficient of drag ( assume 0.015) 

C = local chord 

<$0 = ave. of s.s and p.s. trailing edge boundary layer thickness 
C w = 0.169 from wake behind a circular cylinder 
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go i. x o/u 001 


Flow Development In Rotor Passage 
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Langston's Cascade 



Figure 1 


Cp 
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8 0 I/O 


«ro 










I 


Cp 



Total Pressure Loss Coefficient 


Measured (Langston et at. 76-GT-50) 



Computed w/ A.lg. RG Mixing Length 



mid-span 
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Hub Surface Static Pressure 


Measured (Langston el at. 76-GT-50) Computed w/ Algebraic RG Mixing Length 



St x 10 3 
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Conclusions 


# RNG based algebraic model is a viable alternative 
to traditional algebraic models for complex flows. 


• RNG model mimics transitions and near-wall 
damping through the Heaviside function 

• Good boundary layer profiles, loss profiles, pressure 
loadings can be achieved for viscous dominated flow: 

• Stanton number predictions are not yet adequate 
when using the "plain vanilla" form of the model 

• Work is needed to include non-equilibrium effects 
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Applied k-e and Baldwin-Lomax Turbulence 
models for S-ducts 

G.J. Harloff, Sverdrup Technology, Inc. 

LeRC Group, Brook Park, Ohio 

September 16, 1993 

Experimental contributions to paper by 

B.A. Reichert, NASA Lewis Research Center 
Cleveland, Ohio 

S.R. Wellborn, Iowa State University 
Ames, Iowa 


Objectives 

Provide CFD validation data with secondary flows 
Predict and compare with experimental flow fields 
Assess algebraic turbulence model 

Compare algebraic and K - e turbulence model results 


Introduction 

Diffusing S-ducts common on aircraft: 

727, L-1011, F-16, F-18 

S-duct fabricated and tested 

CFD: 3-D full Navier-Stokes, two turbulence models 

Research extends previous PNS and FNS studies 

Overall features qualitatively correct: total pressure 
velocity vectors, exit vortices 

Boundary layer separation details generally 
not correct 

Turbulence models and/or grid resolution usually 
cited for lack of agreement 

Code Selected: PARC3D 

Originally ARC3D - Pulliam and Steger 

Modified for internal flows - Cooper 

Full 3-D Navier-Stokes equations in Reynolds (mass) 
average form 

Beam-Warming approximate factorization 

Multi-block for computer efficiency 

Cow Reynolds number K - e turbulence model 
of Speziale employed - Nicliols 

' Algebraic Baldwin-Lomax turbulence model 



Bqpqdary Conditions 

JSo slip at wall 

Total pressure and temperature specified at inlet 
Static pressure specified at exit 
Symmetry about x - z plane 
b = 0 on wall 



Ji, e : zeroth-order extrapolation at inlet, 
outlet, and centerline 
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Duct Test Conditions and Ireomelry 

Re = 2.6xl0 6 

M 1 =0.6 

D 1? D 9 = 8.04, 9.90 in 

= 1 52 

Al 

R = 40.2 ill 
e max / 2 = 30 

measurement planes S/D ^ = -0.5, 5.73 
upstream and downstream pipes =3.75D ^ 



ll x 15 
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A Algebraic turbulencemodel 


Experiment&Computation 


Experiment , Computation 


A^ial MaPh namlw pantaws at s/0| - 5.73 


k-e turbulence model 
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ftk;£^tur,bulence^model0' 


Axial skin friction coefficient 


0.0030 




0.0030 


0.0020 


0.0010 
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Conclusions 

Computed flow fields agree reasonably well with 
experimental flow field 

K - e turbulence model better predicts pressure field 
than algebraic model 

Both models underpredict length, angular extent 
and axial location of boundary layer separation 

Possible causes include: inappropriate artificial/ 

computed viscosity, especially in separated region 

Improvements are needed to account for strong 
secondary flows with separation 
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Proteus Experience with the Modified MML 
Turbulence Model 


Julianne Conley 
NASA Lewis Research Center 


The Workshop on Computational Turbulence Modeling 
NASA Lewis Research Center 
September 15-16, 1993 
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Overview 


• Background 

• Modification of MML Turbulence Model 

• Test Cases 

• Concluding Remarks 


Background 


• Based on original MML model of Potapczuk (1989) 

• Simple model for flows where Baldwin Lomax model falls short 

• Based on Prandtl’s mixing length theory: 

H, = P ‘ 2 M 
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Modification of the MML Turbulence Model 

Step 1: Get good estimate of x w 

Step 2: Evaluate and modify MML for zero pressure gradient 
boundary layer flows 

Step 3: Modify model for an adverse pressure gradient flow 
Step 4: Combine all features into one general model (MMLPG) 
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Step 1: Get good estimate of x w 

= 

Global approach -- use momentum equation and 2 interior grid points 


dr 

3 >. 


_ dp 

~ 51 



T = ay 2 + by + c 


t, — -tj-Mj’i-yj) dp 

• ‘ = 35 

w 

c = ^{r l + r 2 -aty] + yl) - ^ (y, + y 2 ) ) 


=>T = C 

w 


T = + 


Use local average of x w to avoid large values of y* near separated regions: 

M = °- 1 l x i - d + °- 2 | x , - il + °- 4 N + °- 2 l T i + il + °- 1 | x ( + J 
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Turbulent Viscosity for Flow over a Flat Plate 


MML Baldwin and Lomax 




381 


Found optimum C, at each Re 





\ 

iLer = 

1 -e' A ' 

/<C, 

/ 

f- C, K 

/ 

leap ~ KC 1 ^ 

1 -e' 1 A+ ' 

/>C, 

Found l + cap = fen (9) 

Cap = 

1860- (6.20 x 10 5 ) 9 

l + 

= min (Caner Cap) 
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Step 3: Modifications for adverse pressure gradient flow 

Benchmark test cases : Equilibrium boundary layer flows 
(P. Bradshaw, 1966) 

g 

£/_« x?, p = = constant 

a = 0,-0.15,-0.255 p = 0,1,5 
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Combining with Blending Function: 



. c. v 



'• ■ 

1 

1 

rj|^ H 


/<c, 





our 

o 

* 

II 

♦ 

• 


y* > c 

C, = 

0.4C 4 , 

C 2 = 

5C 3 k 

c, = 

fen (p) 

c 4 = /cn(p,c / ) 


C 4 - C 5 + C 6 c f 


Strength 

of 

pressure 

gradient 

P 

C 3 

C 5 

c 6 

zero 

0 

1.00 

23,300 

-7.75x1 0 6 

mild 

1 

1.25 

30,100 

-1 .1 6x1 0 7 

strong 

5 

1.53 

33,800 

-2.09x1 0 7 
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Step 4: Combine to get one general model, MMLPG 


C 3 = 1.0 
C 5 = 23, 300. 


C 5 = -7.75 x 10 6 


(3 < 0.0 


C 3 = 1.0 + 0.307 P- 0.0391 P 2 
C 5 = 23.200 + 8560P- 1230P 2 
C 6 = -7.75 x 10 6 - 4.51 x 10 6 p + 386, 000P 2 


0.0 < P < 5.34 


C 3 = 1.52 
C 5 = 33, 900 
C 6 = -20, 900 


P > 5.34 


Velocity Defect Profiles - Mild Adverse Pressure Gradient 
MMLPG 




y lb 
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Static Pressure Distribution - Weak Shock Case. FU0.82 
Top Wall 



x/H 


Bottom Wall 
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Test Cases 

Sajben Transonic Diffuser Flows 



Computational Grid 
81x51 
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IIUIIIIIIII 



Static Pressure Distribution - No Shock Case. R=0.862 


Top Wall 




Concluding Remarks 

• Developed MMLPG for adverse pressure gradient flows 

• Shown that MMLPG successfully computes boundary layer 
flows and transonic diffuser flows 

• Future work: continue validation 

modifications for separated flow 
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National Aeronautics and 
Space Administration 

Lewis Research Center 


Turbulence Model Experiences for a 
Round-to-Rectangular Transition Duct 

Workshop on Computational 
Turbulence Modeling 

September 15-16, 1993 


Charles E. Towne 

NASA Lewis Research Center 
Cleveland, OH 



TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 
Lewis Research Center 


Outline 

• Proteus code description 

• Geometry and flow conditions 

• Numerical details 

• Turbulence model modifications 

• Convergence history 

• Computed results and 
comparison with experimental data 

• Concluding remarks 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 
Lewis Research Center 


Proteus Navier-Stokes Code 

• Reynolds-averaged, unsteady, compressible 
Navier-Stokes equations 

• Strong conservation-law form 

• Fully-coupled ADI solution procedure, 
Beam-Warming generalized time differencing 

• Second-order central differencing for 
spatial derivatives 

• Implicit steady/unsteady boundary conditions 

• Convection and diffusion terms linearized 
using second-order Taylor series expansion 

• Generalized nonorthogonal body-fitted coordinates 

• Baldwin-Lomax and Chien k-e turbulence models 



TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Proteus Navier-Stokes Code 

• 2-D, axisymmetric w/wo swirl, or 3-D flow 

• Thin-layer and Euler options 

• Wide variety of boundary conditions 

• Constant stagnation enthalpy option 

• Constant-coefficient or adaptive artificial viscosity 

• First- or second-order time differencing 

• Variety of time step selection methods 

• Output files for CONTOUR and PLOT3D 

• Highly vectorized for Cray computers 

• Extensively commented 

• Three-volume documentation set 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Circular-to-Rectangular Transition Duct 



STATION 1 


TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Flow Conditions 

• M ref = M inlel = 0. 2 to simulate 
incompressible experiment 

• Re ref = (puR/fi) in i el = 195,000 

• Constant c p and c v , with y = 1.4 

• Constant //, at 519 °R 

• Constant stagnation enthalpy 
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transition duct flow 


National Aeronautics and 

Space Administration 
Lewis Research Center 



Total Pressure Contours 




TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Boundary and Initial Conditions 


Boundary conditions 

Boundary Conditions Specified 


Inlet 

Exit 

Outer wall 
Centerline 

0 = 0 
0 = 90 


Pt -(Pt)cip- u x =0,v = w = 0 
m = rhap, u x = v x = w x = 0 
p r =0, u = v = w = 0 
(p, u, v, w) CL = ave (p, u, v, w) ct+ , 

9 

p, = u, = w r = 0, v = 0 
p, = u, = v r = 0, w = 0 


Initial conditions 

P = Poo. “ = (“apL/t,. v = w = 0 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Numerics 


• 61 (x) x 31(0) x 50(r) and 61 x 31 x 99 meshes 


• Local time step, with: 


Mesh 

Time Level 

CFL 

61x31x50 

1-1001 

1 


1001-2001 

10 


2001 - 3001 

10 


3001-4001 

20 

61x31x99 

1 - 1501 

1 


1501-3001 

10 


3001 -4501 

10 


4501 -6001 

15 


6001 -7501 

15 


• Cray Y-MP CPU time — just under 5 x 10" 5 
sec/grid point/time step, or 5.2 and 18.8 hours 
for the two cases 



TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Computational Mesh 
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Distant* along coordinate line. 



TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 
Lewis Research Center 


Baldwin-Lomax Turbulence Model 


Inner region 

• p t = pl 2 \Q.\Re re f 
Outer region 

• fl, = KC cp pF Kteb F wake^ref 

where F wake = {y^^F,^, F max = max [y n O(l - 
and {y^max is the corresponding y n 

• Model based on boundary layer flows, without accounting 
for secondary flows and streamwise vorticity 

• At each streamwise station, the search region for F^ was 
limited to the value of (yj^ in the y = 0 plane 

• Rationale — to prevent the secondary vortices, at the 
“ends" of the cross-section, from affecting the p t values 


Effect of Turbulence Model Modification 



TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 



0.0 


0.0 0.2 0.4 0.6 0.8 1.0 

Baldwin-Lomax vorticity parameter. F 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Spaca Administration 

Lewis Research Center 



Effect of Turbulence Model Modification 





TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Effect of Turbulence Model Modification 


Modified n, model 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 



Effect of Turbulence Model Modification 





TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Mass Flow Rate Convergence 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Convergence History 





TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Peripheral Static Pressure Distribution 



• l.o - 0.6 - 0.2 0.2 0.6 1.0 

Circumferential distance, s/s re( 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Centerline Static Pressure Distribution 
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Secondary Velocity Vectors 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 
Lewis Research Center 


Secondary Velocity Vectors 
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Total Pressure Contours 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Total Pressure Contours 




TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 


Total Pressure Profiles Along Semi-Minor Axis 



0 0 0.4 0.8 1.2 

Total pressure coefficient, (p T -pj)/O rc t 
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TRANSITION DUCT FLOW 


National Aeronautics and 
Space Administration 

Lewis Research Center 



Total Pressure Profiles Along Semi-Major Axis 



0.0 0.4 0.8 1 2 

Total pressure coefficient. (p r pi)/Q f# , 



TRANSITION DUCT FLOW 
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Space Administration 

Lewis Research Center 


Concluding Remarks 

• The Proteus 3-D Navier-Stokes code has been used to compute 
the flow through a round-to-rectangular transition duct 

• The search region for in the Baldwin-Lomax turbulence model 
was limited to prevent excessively large fi, values 

• The code correctly captures the basic physics — the generation 
of secondary flows and the resulting distortion 

• Agreement between computed and experimental results is good 
through the transition section 

• The computed secondary flows are overly damped in the 
downstream straight section 

• Further work is needed to investigate the effects of mesh 
resolution and turbulence model 
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Proteus experience with three different types of turbulence 

models 


Trong T. Bui 

NASA Lewis Research Center 
Cleveland, Ohio 
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Overview of turbulence models used 


1 . Algebraic model 

• Baldwin-Lomax (BL) 

2. One-equation model 

• Baldwin-Barth one-equation (BB) 

3. Two-equation models 

. Chien k-e (CH) 

• Launder-Sharma k-e (LS) 


NASA Lewis Research Center 
Cleveland. Ohio 


INTERNAL FLUID MECHANICS DIVISION 


Contents 

1 . Overview of turbulence models used 

2. Practical issues 

• CPU time requirement 

• Initialization of multi-equation models 

• y + dependency 

• Compressibility corrections 

• Where the turbulence models work well, and where they do not 
work too well 

3. Comments 


NASA Lewis Research Center 
Cleveland. Ohio 
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uhu time requirement 


Not a major issue in using multi-equation turbulence models: 


(CPU time in sec/iter/grid point) 


Turbulence 

models 

Incomp, flat 
plate 

Comp, flat 
plate 

S-duct 

3-D shock/ 
B.L. inter. 

B.- L. 

6.031E-5 

9.042E-5 

5.168E-5 

7.888E-5 

B.- B. 

6.650E-5 

9.436E-5 

N/A 

N/A 

Chien 

6.822E-5 

9.836E-5 

6.506E-5 

9.757E-5 

L.- S. 

8.667E-5 

11.77E-5 

N/A 

N/A 


• For flat plates, the Chien k-e model takes only about 10% more CPU 
time than the B.- L. model 

• For 3-D S-duct and shock/B.L. inter., the Chien k-e model takes about 
25% more CPU time than the B;- L. model 


NASA Lewis Research Center 
Cleveland. Ohio 


INTERNAL FLUID MECHANICS DIVISION 


Initialization of multi-equation models 

Initial profiles of turbulence variables are needed to start the time 
marching process. The following steps have been found to work well: 

1 . m is computed using the BL model. 

2. BB model - k 2 /e is computed using and the BB formula. 

3. k-e models — e is calculated from the local equilibrium assumption, 
and k is calculated using e, p t , and the CH damping function f^. 


NASA Lewis Research Center 
Cleveland. Oho 
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Automatic Initial k Profile for Incompressible Flat Plate TBL 




Initial k profile Converged k profile 


NASA Lewis Research Center 
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Initial k Profile for Turbulent S-Duct Flow 


Initial k profile 


Converged k profile 
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Cleveland. Ohio 
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Initial k Profile for Glancing 
Shock Wave/TBL Interaction 



NASA Lewis Research Center INTERNAL FLUID MECHANICS DIVISION 

Cleveland. Ohio 


Converged k profile 


Initial k profile 


Y + dependency 


y + is needed to feed the Van Driest style damping functions in the near 
wall region. Not desirable for two main reasons: 

1 . Geometrical considerations: 


• Need to identify the nearest solid walls. This is difficult for complex 
3-D geometries with multiple solid boundaries. 

• Need to compute distances to solid walls. This is difficult to do for 
unstructured grids. 

• What is y + in base/wake flows? 


2. Ambiguity in the definition of y + for variable-property flows: 


i Re r 

w r 


• ■ ft* 

• Local or wall values for v-, and 


However, models that use y + appear to be more stable numerically than 
the ones that do not. 
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Effects of y + Computation, BL 
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Effects of y + Computation, BB 




Re, Re, 


NASA Lewis Research Center 
Cleve'-and. Ohio 


internal fluid MECHANICS DIVISION 


408 


Effects of y + Computation, CH 
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Effects of Compressibility Corrections, CH 



Re i 



0 5.000 10' 1.00010' 1.500 10' 2.000 10* 2.500 10* 3.000 10* 

Rfli 
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Contents 


1 . Overview of turbulence models used 

2. Practical issues 

• CPU time requirement 

• Initialization of multi-equation models 

• y + dependency 

• Compressibility corrections 

— ► • Where the turbulence models work well, and where they do not 
work too well 

3. Comments 
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Where the turbulence models work well 


Attached, thin layer flows with little or no pressure gradients 
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Where they don’t work too well 


Separated/Reattached flows and flows with large pressure gradients 






NASA Lewis Research Center 
Cleveland. Ohio 


INTERNAL FLUID MECHANICS DIVISION 


Comments 

• Two-equation k-e models can be used without serious CPU time 
penalty. However, no real gain in predictive capability compared with 
algebraic formulations. 

• To improve the predictive capability, need to extend the two-equation 
formulations with algebraic Reynolds stress models or anisotropic k-e 
models. 

• y + dependency limits the generality of turbulence models, and it should 
be removed if possible (preferably without significant cost in the 
stability and accuracy of the models). 

• Models are needed for the Reynolds heat flux terms. 
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Cleveland. Ohio 


INTERNAL FLUID MECHANICS DIVISION 


411 


Page intentionally left blank 



Several Examples Where Turbulence Models 
Fail in Inlet Flow Field Analysis 

Bernhard H. Anderson 
NASA Lewis Research Center 
Cleveland, OH 
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Inlet Flow Field Analysis 

Computational Uncertainties 


• Turbulence Modeling for 3D Inlet Flow Fields 

(1) Flows Approaching Separation 

(2) Strength of Secondary Flow Field 

(3) 3D Flow Predictions of Vortex Liftoff 

(4) Influence of Vortex-Boundary Layer Interactions 

• Vortex Generator Modeling 

(1) Representation of Generator Vorticity Field 

(2) Relationship Between Generator and Vorticity Field 


Inlet Flow Field Studies 

Goals and Objectives 


To advance the understanding, the prediction, and the control 
of intake distortion, and to study the basic interactions that 
influence this design problem. 

• To develop an understanding of and predictive capability 
for the aerodynamic properties of intake distortion and 
its management. 

• To establish a set of design guidelines to maximize the 
effectiveness of vortex flow control for the management, 
intake distortion 
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Inlet Flow Field Benchmark Data Sets 

Turbulence and Vortex Generator Modeling 

• Fraser Flow A, Stanford Conference 1968 

• 727/JT8D-100 S-Duct Confirmation Experiment, 1973 

• Univ. Tennessee Diffusing S-Duct Experiments, 1986 & 1992 

• Univ. Washington TD410 Transition Duct Experiment, 1990 

• M2129 Imperial College Laser-Doppler Experiment, 1990 

• DRA-Bedford Experiments on the M2129 Intake S-Duct 

(2) DRA Surface Pressure and Engine Face. Experiment, 1989 

(3) DRA Phase IB Hot-Wire Flow Experiment, 1990 

(4) DRA Phase 2 Yawmeter Flow Experiment, 1991 

(5) DRA Phase 3 Vortex Generator Experiment, 1992 

• TD118 Bi-Furcated Transition Duct Experiment, 1994 


Reduced Navier-Stokes Analysis 

RNS3D Computer Code 

• Velocity decomposition approach, Briley and McDonald (1979 & 1984) 

(1) Conservation form of the vorticity transport equation 

(2) Mass flow conservation, m — j^pUpdA = constant 

• Non-orthogonal coordinate system, Levy, Briley and McDonald (1983) 

• Arbitrary geometry gridfile description, Anderson (1990) 

(1) Recluster existing gridfile mesh distribution 

(2) Redefine the centerline space curve 

(3) Alter cross-sectional duct shape 

• McDonald Camarata turbulence model 
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Full Navier-Stokes Analysis 

PARC3D Computer Code 


• Originally developed by Pulliam & Steger as AIR3D (1980) 

(1) Conservation form of the governing equations 

(2) Beam & Warming approximate factorization algorithm 

(3) Central differencing within a curvilinear system 

• Addition of Jameson artificial dissipation by Pulliam, ARC3D (1981) 

• Developed for internal flow by Cooper, PARC3D (1987) 

(1) Baldwin-Lomax turbulence model 

(2) Diewert approximation to turbulence model in the 
reverse flow region of flow field 


Fraser Flow A, Stanford Conference 1968 

Geometry and Mesh Definition 
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Distance from the Wall, inches 


Fraser Flow A, Stanford Conference 1968 

Comparison of Turbulence Models 
nz = 49, y + = 1.17 


— RNS Analysis, MC Model 

— FNS Analysis, BL Model 

— FNS Analysis, (k — e) Model 
O Data 



Velocity Ratio, £// Uj 


Fraser Flow A, Stanford Conference 1968 
Comparison of Turbulence Models 
nz = 49, y + = 1.17 



5 r 
4 
3 
2 

1 
0 


RNS Analysis, MC Model 
FNS Analysis, BL Model 
FNS Analysis, ( k - t) Model 


n ^ n - J 


O Data 


OO 


(a) Incompressible Shape Factor (HI) Distribution 



(b) Wall Skin Friction (C/) Distribution 


417 


Fraser Flow A, Stanford Conference 1968 

Effect of y + on Flow Field Solution 


ni = 49, y' = 0.584 

ni = 49,y* = 1.168 

ni = 49, y * = 5.739 


O Data 


O |Q 2_ 


oo 


(a) Incompressible Shape Factor (HI) Distribution 



(b) Wall Skin Friction (Q) Distribution 

RNS Analysis, McDonald-Camarata Model 


Fraser Flow A, Stanford Conference 1968 

Effect of Mesh Resolution on Flow Field Solution 


X 



S. 

a 


t n 


5 

4 

3 

2 

1 

0 


— ni = 199, y* = 0.477 

— ni = 99, y* = 0.520 

— ni = 49, y* = 0.584 
O Data 


■O-C— ' 9* — 0= 



(a) Incompressible Shape Factor (HI) Distribution 



HNS Analysis, McDonald-Camarata Model 
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Fraser Flow A, Stanford Conference 1968 

Conclusions 

(1) The current generation of turbulence models were unable to 
predict the complete state of the diffuser boundary layer 
approaching flow separation. 

(2) Both near wall and mesh resolution separately played an 
important role in accurate solutions to wall skin friction 
distribution in flows characterized as "approaching separation", 
but had little effect on the solution for the incompressible 
shape factor development. 

(3) It is important that grid independent solution be demonstrated 
before judgements about the turbulence models be stated, and 
that the complete state of the wall boundary layer be 
considered within this evaluation. 


727/JT8D-100 Inlet S-Duct 

Geometry and Mesh Definition 
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727/JT8D-100 Inlet S-Duct 

Engine Face Flow Field 
Generator Config. 12 




Analysis 

Without Engine Dome 


RNS Analysis, McDonald-Camarata Model 


727/ JT8D-100 Inlet S-Duct 

Engine Face Ring Distortion Characteristics 
Generator Config. 12 


O Experiment (With Engine Dome) 
_ Analysis (Without Engine Dome) 



Engine Face Radius, Inches (Model Scale) 



Radial Distortion 


60° Sector Circumferential Distortion 


RNS Analysis. McDonald-Cnniarnta Model 
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727/JT8D-100 Inlet S-Duct 

Conclusions 


(1) The current turbulence models predict the overall performance 
level of vortex generator installation remarkably well, although 
much of the detailed flow structure was not resolved. 

(2) Turbulence models in 3D inlet flow field analysis can also be 
evaluated on the basis of standard engine performance parameters, 
such as radial and circumferential ring distortion descriptors, 
which provide a sensitive discriminator measuring the state of 

the overall compressor face flow field. 


Univ. Tennessee Diffusing S-Duct 

Geometry and Mesh Definition 
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Univ. Tennessee Diffusing S-Duct 

Topology of Vortex Liftoff 



Experiment 


RNS Analysis. McDonald-Camarata Model 


Univ. Tennessee Diffusing S-Duct 

Topography of Vortex Liftoff 



RNS Analysis, McDonald-Camarata Model 
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Univ. Tennessee Diffusing S-Duct 
Total Pressure Coefficient Contours 
Without Vortex Generators 



Experiment Analysis 


RNS Analysis, McDonald-Camarata Model 


Univ. Tennessee Diffusing S-Duct 
Total Pressure Coefficient Contours 
With Vortex Generators 



Experiment Analysis 


RNS Analysis, McDonald-Camarata Model 
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Univ. Tennessee Diffusing S-Duct 

Effect of y + on Engine Face Flow Field 
Total Pressure Coefficient Contours 




y + = 0.5 


RNS Analvsis. McDonald-Camarata Mode! 


Univ. Tennessee Diffusing S-Duct 

Effect of y + on Circumferential Distortion 


y + — 0.5 
j) + = 8.5 



RNS Analysis, McDonald-Camarata Model 


424 


Univ. Tennessee Diffusing S-Duct 

Conclusions 

(1) Initial value space marching 3D RNS procedures adequately 
described the topological and topographical features of 3D flow 
separation associated with vortex liftoff. 

(2) The current turbulence models predicted the overall structure of 
vortex generator installation remarkably well, although much of the 
detailed flow structure was not resolved. 

(3) Adequate near wall resolution was necessary to obtain an accurate 
solution of the phenomena of vortex lift-off. 

(4) Circumferential ring distortion is a sensitive discriminator in 
measuring the state of the engine face flow field. 


Univ. Washington TD410 Transition Duct 


Geometry and Mesh Definition 



MM 


illija 


S-' 














fill 

wMM 

W0: 


sfe . v- •• 
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Univ. Washington TD410 Transition Duct 

Case Definitions 


Case 

Grid 

Total 

CPU (min) 

F+ 

td410.$ 

199 x 121 x 521 

12,545,159 

149.0 

0.565 

td410.1 

99 x 121 x 521 

6,241,059 

73.4 

0.623 

td410.3 

49 x 121 x 521 

2,851,849 

36.0 

0.568 

td41 1.1 

99 x 91 x 521 

4,693,689 

57.1 

0.623 

td412.$ 

199 x 61 x 521 

6,324,419 

72.9 

0.565 

td412.1 

99 x 61 x 521 

3,146,319 

38.2 

0.623 

td412.3 

49 x 61 x 521 

1,557,269 

18.7 

0.568 

td412.4 

49 x 61 x 401 

1,198,589 

14.7 

0.568 


Univ. Washington TD410 Transition Duct 

Effect of Radial Grid Resolution 



(a) Velocity Profiles Along the Minor Axis Surface 



(b) Velocity Profiles Along the Major Axis Surface 


RNS Analysis, McDonald-Camarata Model 
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Wall Skin Friction, C,x 10> Wall Skin Friction, C,x 10> 


Univ. Washington TD410 Transition Duct 
Effect of Radial Grid Resolution 
Wall Skin Friction 

m = 199, y* = 0.565 

.5 - 
.4 - 
.3 - 
.2 - 
.1 - 
.0 

(a) Skin Friction Q along the Minor Axis Wall 




Axial Centerline Distance, XIR , 

(b) Skin Friction Q along the Major Axis Wall 


RNS Analysis. McDonald-Camarata Model 


Univ. Washington TD410 Transition Duct 

Comparison with Experimental Data 
Velocity Contours 


uiu Mk 



RNS Analysis, McDonald-Camarata Model 
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Univ. Washington TD410 Transition Duct 

Conclusions 

(1) The current generation of turbulence models predicted the 
overall development of vortex formation reasonably well, 
although there were important discrepencies which could not 
not be explained as inadequate near wall or mesh resolution. 

(2) Radial mesh resolution had the largest impact in the region 
along the major axis where the vortex pair was formed. 

(3) It is important that grid independent solution be demonstrated 
before judgements about the turbulence models be stated. 

(4) Fully 3D grid independent solutions were achieved with a 
Reduced Navier-Stokes analysis. 


DRA M2129 Diffusing Inlet S-Duct 

Geometry and Mesh Definition 
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DRA M2129 Diffusing Inlet S-Duct 

Performance Characteristics 


1.00 


^ 0.98 
>1 

S 0.96 
o 

| 0.94 
| 0.93 
£ 0.90 



~ RNS Analysis, MC Model 
• FNS Analysis, BL Model 
O Data 



(b) Engine Face DC i0 Distortion 


DRA M2129 Diffusing Inlet S-Duct 

Separation Characteristics 

RNS Analysis, MC Model 

• FNS Analysis, BL Model 
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Static Pressure Ratio, p/Ptt 


DRA M2129 Inlet S-Duct 

Wall Static Pressure Distribution 
AGARD Test Case 3.1 

RNS Analysis, MC Model 

FNS Analysis, BL Model 

O Data 

.9 r 



.6 - 
.5 - 
A - 


(a) 6 = 0.0° Surface Element. 



Axial Centerline Distance, X/R, 

(b) 0 = 180.0° Surface Element. 


DRA M2129 Diffusing Inlet S-Duct 

Engine Face Distortion Characteristics 
AGARD Test Case 3.1 




(b) 60°-Sector Circumferential Distortion 
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DRA M2129 Inlet S-Duct 

Conclusions 


(1) Both Full Navier-Stokes (FNS) and Reduce Navier-Stokes (RNS) 
analyses adequately describe the overall flow physics of vortex 
liftoff, but consistently predict the location of liftoff further 
downstream in the duct inlet than was indicated by data. 

(2) The current generation of turbulence models were unable 
to desribed the influence of separation on the main pressure 
field for "strong" vortex liftoff interactions. 

(3) The current generation of mixing length turbulence models 
give remarkable good performance results, while the existing 
discrepancies between data and analysis can be attributed 
primarily to the over prediction of the liftoff location. 


TD118 Bi-Furcated Transition Duct 

Geometry and Mesh Definition 



.... , , 

f * ‘ : ^ 
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TD118 Bi-Furcated Transition S-Duct 

Research Objectives 


• To demonstrate diffuser duct technology advancement by using 
CFD to design a "conventionally shorter" transitioning S-duct 
configuration for application towards high speed inlet systems. 

• To develop a computational protocol whereby turbulence model 
evaluations can be made between different computer codes. 

• To develop a benchmark data set to evaluate CFD analysis and 
turbulence models, which cover fundamental flow phenomena as 
well as overall flow field physics as determined by standard 
engine performance parameters. 


TD118 Bi-Furcated Transition S-Duct 

Effect of Turbulence Model on Inlet Performance 


Analysis 

Turbulence Model 

Ptavel Pto 

DH 

DC i0 

RNS3D 

McDonald-Camarata 

0.955 

0.167 

0.116 

PARC3D 

Baldwin-Lomax 

0.959 

0.163 

0.135 

PARC3D 

P.D. Thomas 

0.970 

0.153 

0.073 

PARC3D 

Launder-Spaulding 

0.975 

0.150 

0.177 


DH — {Ptnax ~ Pt m \ n )l Pt„vc 
DCfrO — ( Ptave Ptnlillfio)! (Javc 
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TD118 Bi-Furcated Transition S-Duct 

Effect of Turbulence Model on Engine Face Flow Field 




(n) RNS Analysis, McDonald-Camarata Model 


(c) FNS Analysis, P. D. Thomas Model 




(b) FNS Analysis, Baldwin-Lomax Model 


(d) FNS Analysis, Launder-Spalding ( k - e) Model 


TD118 Bi-Furcated Transition S-Duct 
Effect of Turbulence Model on Wall Boundary Layer 
0 = 90.0° Surface Element 


X 



E 

o 

u 

c 


4 - 
3 - 
2 - 
1 - 
0 - 


Analysis 

- FNS, BL Model 

- FNS, PDT Model 

- FNS, LS (A - i) Model 
O RNS, MC Model 



(a) Incompressible Shape Factor HI Distribution. 



(b) Wall Skin Friction Q Distribution. 
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TD118 Bi-Furcated Transition Duct 

Effect of Turbulence Model on Ring Distortion 


Analysis 

— FNS, BL Model 




(b) 60°-Sector Circumferential Distortion 


Inlet Flow Field Analysis 

Concluding Remarks 


(1) Difficulties in complex 3D flow fields often arise because 
fundamental 2D aerodynamic interactions have not been 
adequately resolved. 

(2) Near wall (y + ) and radial mesh resolution (nz) play an 
important role in fundamental 2D and complex 3D flow 
field analysis. 

(3) Judgements about turbulent models should not be stated until 
grid independent solutions have been established. 

(4) Adequateness of turbulence models in inlet flow field analysis 
should also be made on the basis of fundamental performance 
parameters used to quantify the "goodness" of the flow 
entering the engine. 
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